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PREFACE 


The  production,  testing,  and  disposal  of  various  chemical  and 
biological  substances  at  RMA  has  lead  to  suspected  contamination  of 
land  areas,  water,  ecology,  buildings,  and  equipment.  Little  infor¬ 
mation  exists  as  to  the  concentration  levels,  overall  extent,  or  the 
diffusion  process  of  these  substances  within  the  ecosystem.  The  need 
to  identify  and  assess  the  potential  hazards  created  by  the  presence 
of  toxic  substances  in  the  environment  at  RMA  has  been  recognized, 
and  efforts  are  underway  for  interim  containment  of  off-post  migration. 
A  comprehensive  survey  will  be  initiated  to  gather  and  analyze  all 
necessary  information  which  would  enable  a  rational  assessment  of  the 
degree  of  contamination  on  the  post  and  provide  the  data  base  needed 
for  the  design  of  remedial  measures. 

The  purpose  of  this  study  is  to  provide  the  guidelines  for  soil 
and  water  sampling  procedures  and  field  data  gathering  and  recording 
protocols  to  be  used  in  the  comprehensive  survey  at  RMA.  The 
approach  developed  is  aimed  at  obtaining  the  maximum  likelihood  of 
intercepting  the  contaminated  horizons  in  three  dimensions  and  thus 
gaining  the  optimum  benefit  for  long-run  objectives  from  the  field 
measurements  with  minimum  field  effort  and  thus  reducing  time  and 
cost  by  minimizing  the  possibility  of  gathering  extraneous  data. 


INTRODUCTION 


A  pollutant  is  any  compound  or  substance  that  is  undesirable 
where  it  is.  A  contaminant  or  pollutant  may  become  undesirable  by 
being  introduced  into  an  environment  that  is  foreign  to  it.  The 
pollutant  is  usually  toxic  to  the  host  environment  and  will  disturb 
the  equilibrium  state  of  the  ecosystem.  The  degree  of  concern  of  the 
presence  of  contamination  depends  upon  the  toxicity  level  to  living 
organisms,  especially  man. 

Pollutants  may  be  either  man-made,  released  as  byproducts  of 
industrial  or  military  processes,  or  are  naturally  occurring  and 
become  toxic  when  concentrated,  for  example  salt. 

Transport  processes  often  create  a  pollutant  simply  by  moving  a 
compound  to  a  different  environment.  Sources  of  pollutants  are 
usually  classed  as  "point"  and  "nonpoint."  Control  measures  are 
simplified  whenever  the  pollution  plume  is  traced  to  a  point  source, 
as  compared  with  diffuse  or  nonpoint  pollutants  which  discharge  into 
waterways  by  dispersal  pathways. 


FACTORS  INFLUENCING  CONTAMINANT  BEHAVIOR 


The  success  in  locating  the  most  probable  areas  of  contamination 
and  in  tracing  the  breadth  (source  and  plume)  of  such  contamination 
within  the  ecosystem  is  dependent  upon  an  accurate  interpretation  of 
the  means  of  contaminant  emplacement.  The  dispersion  and  the  accumu¬ 
lation  of  pollutants  in  the  environment  is  dependent  on  a  number  of 
inter-related  factors.  Some  of  the  principal  factors  are: 

The  Physical  and  Chemical  Properties  of  the  Contaminants  -  The 
important  chemical  properties  of  the  pollutants  are  the  following: 
the  contaminant  phase,  concentration,  chemical  stability,  water 
solubility,  ion  exchange,  adsorption,  equ i 1 i br ium  wi th  the  host 
environment,  and  regional  geochemical  properties.  These  properties 
determine  the  chemical  behavior  of  the  contaminants,  whether  they 
are  present  as  a  solution  in  water  or  are  dispersing  as  a  solid. 

The  Form  of  Transport  -  The  form  of  transport  determines  the 
dispersion  process  and  influences  accumulation.  It  is  possible  by 
evaluating  the  manner  that  the  pollutants  were  transported  and 
knowing  their  pathways  to  predict  the  major  contaminated  areas. 
Contaminants  of  greatest  concern  at  RMA  are  transported  either: 

(1)  as  solid  particles  (particulate)  occurring  in  a  crystal ized 
form  or  combined  with  dust  or  soil  particles,  or  (2)  as  a  fluid 
in  solution  with  water.  For  example,  the  movement  of  pesticides 
on  sediment  is  affected  by  the  enrichment  process.  Pesticides  are 
adsorbed  primarily  on  organic  soil  colloids  which  remain  in  suspen¬ 
sion  longer  than  the  coarse  soil  particles.  Most  mobile  herbicides, 
on  the  other  hand,  rarely  contaminate  the  groundwater  because  the. 
adsorption  capacity  in  most  soils  above  the  water  table  is  sufficient 
to  retain  the  chemicals  until  they  decompose.  However,  in  an. area 
with  a  high  percolation  potential,  groundwater  contamination  is 
poss i bl e. 

The  partitioning  of  the  compounds  between  water  and  soil  trans¬ 
port  is  necessary  in  determining  the  mass  movement  of  chemicals  in 
run-off.  The  quantity  of  a  compound  dissolved  in  water  is  obtained 
from  water  solubilities  of  the  chemicals.  The  amount  adsorbed  at 
different  solution  concentrations  can  be  calculated  by  mixing  solu¬ 
tions  of  the  chemicals  with  known  quantities  of  soil.  (This  study 
is  currently  in  progress  for  DIMP  by  the  Aerojet  Ordnance  and 
Manufacturing  Company,  Downey,  California.) 
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The  transporting  medium  for  particulate  contaminants  is  either 
iwi nd  or  surface  water  movement.  The  contaminants,  therefore,  con¬ 
centrate  in  areas  controlled  by  these  transporting  media  —  on  or 
close  to  the  ground  surface.  The  transporting  media  for  liquid 
contaminant,  on  the  other  hand,  is  either  surface  run-off,  subsurface 
run-off,  or  groundwater  run-off.  These  contaminants  tend  to  concen¬ 
trate  much  deeper  in  the  unconsolidated  mantle. 

The  Characteristics  of  the  Physical  Site  -  The  characteristics 
'of  the  site  include:  CD  Geologic  properties  of  the  aquifier  such 
as  texture,  porosity,  permeability,  transmissivity,  heterogeneity  of 
the  subsurface  regolith,  thickness  of  the  unsaturated  zone  of  aera¬ 
tion,  thickness  of  the  saturated  groundwater  zone  boundary  conditions, 
and  (2)  The  geomorphology,  the  topography,  and  the  surficial  cover. 

This  includes  slope,  aspect,  surface  relief  (high  rel i ef  is  indicative 
of  high  sediment  delivery  ratio),  channel  density,  the  condition  of  the 
channels  (whether  clogged  or  open,  meandering,  or  straight),  stream 
gradients,  the  rel ief/1 ength  ratio.  The  geomorphology  of  a  terrain 
provides  potential  channels,  sinks,  or  traps  for  pollutant,  and  (3) 

The  stratigraphy  and  character  of  the  sediments  in  the  aquifer.  This 
would  include  the  ionic  exchange  capacity  of  the  minerals  in  clay  and 
soil,  which  has  a  direct  bearing  on  contaminant  sorption  and  leaching. 
The  type  of  both  primary  and  secondary  minerals  present  and  the  oxida¬ 
tion  reduction  reactions  which  affect  some  of  the  pollutants  and  the 
leach  ability  of  the  sediments. 

Cl  i mate  -  Humid  environments  are  those  in  which  the  annual  pre¬ 
cipitation  exceeds  evapotranspi rat  ion  rates.  They  have  a  net  water 
surplus  which  transports  the  mobile  compounds  faster  and  to  greater 
depths.  In  arid  climates  such  as  RMA,  much  of  the  water  is  removed  by 
evapotranspi  rat  ion ;  and  soil  moisture  is  brought  to  the  surface  by 
capillary  action.  Contaminated  groundwater  may  be  transported  to  the 
surface  and  the  contaminants  are  deposited  on  the  ground  surface.  The 
duration  and  intensity  of  storm  activity  influences  groundwater  recharge 
and  possible  flushing  of  the  contaminant  from  the  soil.  Flow  is  upward 
during  periods  of  high  evapotransp i rat i on ,  provided  there  is  no  ponding 
and  downward  when  water  is  available  (as  under  the  water  ponds  and  bog 
areas).  The  micro  climate  also  influences  the  viscosity  and  surface 
tension  properties  of  soil  water  solution,  which  is  temperature- 
dependent.  Soils  drain  slowly  during  cold  weather,  and  water  is 
retained  as  storage  when  frozen.  As  a  solution  freezes,  salts  con¬ 
centrate  in  the  nonfrozen  liquid,  and  it  is  not  unusual  to  find  brines 
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in  partially  frozen  soils.  These  brines  become  selectively  mobile 
because  of  their  hydraulic  properties.  During  the  spring  after  the 
frozen  ground  thaws  at  RMA,  much  of  the  groundwater  recharge  results 
from  snowmelt.  If  the  concentrations  of  pollutants  in  soil  water 
reaches  high  enough  levels,  precipitation  (or  salting  out)  may  occur. 

The  Geohydrolic  Conditions  -  Hydrologic  transport  is  the  move¬ 
ment  of  compounds  in  or  by  water  and  may  occur  along  surface  or  sub¬ 
surface  channels  at  RMA.  Hydrologic  transport  is  therefore  dependent 
upon  the  amount  of  water  and  the  pathways  of:  (1)  Surface  run-off 
traveling  over  the  ground  surface;  (2)  Subsurface  run-off  that  has 
infiltrated  the  surface  soil  and  moves  laterly  through  the  upper  soil 
horizon  within  the  zone  of  aeration;  and  (3)  Groundwater  run-off 
below  the  groundwater  table  in  the  zone  of  saturation.  Each  must  be 
known  before  the  transport  of  pollutants  and  hydrodynamic  dispersion 
can  be  described  in  detail. 

The  hydraulic  conductivity  (K  variable  in  Darcy's  Law)  of  sedi¬ 
ments  is  the  transmissive  properties  of  a  porous  medium.  Generally, 
in  hydrolic  calculations,  this  value  is  averaged  vertically  assuming 
a  homogeneous  aquifer  and  substituting  it  for  transmissivity,  expres¬ 
sed  as  T  =  Kh,  where  h  is  the  saturated  thickness.  The  transmissivity 
can  be  determined  by  pump  tests. 

An  analysis  of  the  dynamics  of  groundwater  flow  permits  a  reliable 
computation  of  the  level  of  contamination  of  the  aquifer  in  space 
through  time.  Changes  in  chemical  concentrations  in  the  groundwater 
system  are  due  to:  (1)  Convective  transport,  expressing  the  amount 
of  dissolved  chemicals  moving  with  the  flowing  groundwater;  (2)  Addi¬ 
tion  or  removal  of  contaminants  through  chemical  and  physical  reactions 
taking  place  in  the  water  or  between  water  and  soil;  (3)  Mixing  with 
recharge  water  of  different  concentrations;  and  (4)  Hydrodynamic 
dispersion  in  which  molecules  or  ionic  diffusion  and  minute  changes 
in  velocity  of  flow  in  the  porous  media  causes  the  paths  to  diverge 
from  the  normal  direction  of  groundwater  movement. 

Grouping  -  Certain  elements  tend  to  group  together  into  compound 
associations  in  a  specific  environment  due  to  similar  migration  capa¬ 
cities  (relative  mobilities).  The  tendency  for  compounds  to  group 
together  simplifies  the  mapping  task  in  that  by  testing  for  either 
indicator  or  pathfinder  elements,  it  is  possible  to  trace  the  spatial 
behavior  of  the  pollutants. 
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Contaminant  Migration  through  Soil  Erosion  -  Several  contam¬ 
inated  areas  exist  on  RMA  which  are  sources  of  ground  surface  pollu 
tion  by  soil  movement.  Field  soil  losses  have  been  estimated  to 
range  from  three  to  eight  tons  per  acre  per  year.  A  ton  of  dry  soil 
is  approximately  one  cubic  yard  in  volume.  Spread  uniformly  over  an 
'acre  of  surface,  it  would  have  a  depth  of  less  than  one-hundredth  of 
an  inch.  Ero'sion  hazards  both  by  wind  and  due  to  rainfall  and  surface 
run-off  are  local  in  nature.  This  can  be  significant  in  areas  which 
are  not  protected  by  a  vegetative  cover  (example  Basins  A,  C,  and 
'most  of  the  Arsenal  during  the  winter  months  when  vegetation  is  dor¬ 
mant).  The  annual  soil  losses  from  sheet  and  rill  erosion  by  rain¬ 
fall  can  be  determined  by  the  Universal  Soil  Loss  Equation  (USLE) . 

The  equation  is: 

A  (RKLSCP) 

Where  A  is  estimated  average  annual  soil  loss  measured  in  tons 
per  acre. 

R  is  the  rainfall  and  run-off  erosivity  index. 

K  is  the  soil  erodibility  factor.  It  is  the  average  soil  loss 
per  unit  of  R  and  depends  on  soil  properties  (see  Figure  ). 

LS  is  dimensionless  topographic  factor  that  represents  the  com¬ 
bined  effects  of  slope,  length,  and  steepness. 

C  is  the  cover  and  management  factor. 

P  is  the  support  -  practice  factor  for  land  use. 

'The  USLE  can  provide  a  measure  of  the  contaminated  sediment 
delivered  to  a  stream. 

The  object  of  the  RMA  environmental  sampling  program  is  to  be 
able  to  identify,  measure,  and  characterize  each  of  the  contaminated 
horizons. 

Generally,  immediately  underneath  a  contaminated  source,  each 
layer  in  the  unconsolidated  mantle  collects  the  contaminants  from 
above,  mixes,  absorbs,  reacts,  and  then  drains  the  soil  to  field 
capacity  to  supply  lechate  to  the  layer  below.  Thus,  the  process 
proceeds 
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FIGURE  i 


STRATIFICATION  OF  THE  UNCONSOL  I  DATES  MANTLE 

When  a  vertical  section  is  made  through  the  unconsolidated  mantle, 
a  series  of  distinct  layers  or  zones  are  observed  which  differ  in  their 
chemical  and  physical  properties.  See  Figure  The  uppermost  A  hori¬ 
zon  varies  from  a  few  inches  to  several  feet  in  thickness  and  repre¬ 
sents  the  eluvial  (washed  out)  layer.  Varying  amounts  of  organic 
material  is  added  to  this.  Materials  removed  from  the  A  horizon  are 
carried  downward  in  solution  and  are  deposited  in  the  B  horizon; 
however,  the  fine  colloidal  clays  are  mechanically  removed  by  descend¬ 
ing  soil  water  (illuvial  layer).  impervious  hardpan  and  claypan 
layers  may  develop  in  the  B  horizon  by  the  accumulation  of  downward 
moving  material  which  include  iron  and  aluminum  minerals,  including 
clays.  The  A  and  B  horizon  together  form  the  solum.  In  arid  cli¬ 
mates,  calcium  carbonate,  sodium  sulphates,  and  other  salts  may 
accumulate  just  below  the  B  horizon  in  the  C  layer.  The  C  layer  or 
parent  material  does  not  exhibit  a  high  degree  of  alteration;  however, 
it  may  display  weathering,  including  oxidation,  hydration,  and  lime 
accumul at  ions. 

A  mature  soil  profile  exhibits  well  developed  horizons;  however, 
layering  may  be  indistinct  or  lacking  in  young  soils.  Topography 
affects  both  soil  profile  development  and  contaminante  emplacement. 

On  steep  slopes,  the  soil  horizons  are  much  thinner  or  may  never 
develop  because  erosion  removes  weather  products  as  rapidly  as  they 
form.  Topography  controls  the  rate  of  surface  run-off,  the  infiltra¬ 
tion  of  water,  the  position  of  the  water  table,  the  depth  of  penetra¬ 
tion  of  chemical  weathering  processes,  and  contaminant  migration. 

Local  and  regional  relief  affect  the  hydrolic  gradient  which  governs 
the  circulation  of  water  through  the  soil  and  unconsolidated  mantle. 

Knowledge  of  the  soil  profiles  developed  within  a  suspected  con¬ 
taminated  area  is  necessary  to  enable  a  constant  relative  sampling 
procedure  to  be  maintained.  All  samples  must  be  taken  from  the  same 
horizon,  otherwise  they  will  not  be  comparable.  Also,  a  uniform 
sampling  procedure  increases  the  probability  of  locating  the  con¬ 
taminated  zones.  See  Figure 

An  analysis  of  the  vegetative  cover  can  provide  valuable  infor¬ 
mation  about  the  character  of  the  contamination  below  ground  surface. 

For  example,  deep-rooted  plants  may  penetrate  the  pollution  plume  in 
the  overburden  and  be  brought  to  the  surface  and  enriched  in  the  soil 
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THE  SOIL  PROFILE 


Figure  2 


Organic,  original  forms  recognized 
Organic,  original  forms  not  recognized 

Mineral,  mixed  with  humus,  dark  colored 

Horizon  of  maximum  eiuviation  of 
silicate  clays,  Fe,  Al  oxides  etc. 

Transition  to  B,  more  like  A  than  B 
Transition  to  A,  more  like  B  than  A 

Maximum  illuviation  of  silicate  clays, 

Fe,  Al  oxides,  some  organic  matter 

Transition  to  C,  more  like  B  than  C 

Zone  of  least  weathering,  accumulation 
of  Ca,  Mg  carbonates,  cementation, 
sometimes  high  bulk  density 
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Nutrient  recycling  is  an  important  factor  in  determining  the  relation1; 
ship  between  vegetation  and  the  soil  that  develops. 


Figure  3 


Bioaccumulation  of  DIMP  by  Plant  Parts  in 
20  ppm  Irrigation  (37  days  from  Original  Inoculation) 


Total  DIMP 

added  to  Pot 

DIMP 

Vol  of 

Wt.  of 

Cone.  In 

Bio- 

Plant 

20ppm 

DIMP 

Tis  sue 

Accumulation 

Part 

Irr. (cc) 

■  (mg) 

(ppm) 

F  actor 

Sugar  Beet  - 

9500 

190 

Root 

45.  6 

2.  28 

Stem 

37.  1 

1. 86 

Leaf 

129.  2 

6.46 

Carrot  - 

9200 

184 

Root 

12.4 

0.  62 

Stem 

6.6 

0.  33 

Leaf 

36.  9 

1. 85 

Bean  - 

9200 

184 

Root 

45.4 

2.27 

Stem 

28.  9 

1.45 

Leaf 

150.0  ' 

7.  50 

Wheat  - 

9200 

184 

Root 

31. 5 

1. 58- 

Stem 

14.  2 

0.  71  ’ 

Leaf 

105.  5 

5.  28 

Figure  4 


Bioac cumulation  Factor 


■Plant  Part 

Bioaccumulation  of  DIMP  by  plant  parts. 

Soil  culture,  37  days  exposure  to  20  ppm  DIMP  in  irrigation  water 
Figure  5 


i  layer  by  decay  of  the  vegetation  in  the  production  of  humus.  The 
| concentrat i on  of  elements  in  plants  is  dependent  on  a  variety  of 
physical,  chemical,  and  organic  factors  which  control  intake  of 
nutritional  solutions.  Geobotanic  survey  may  be  valuable  to  monitor 
contaminant  movement  in  the  soil,  especially  in  areas  where  soils 
profiles  are  young  or  poorly  developed.  See  Figure  Results 
from  biochemical  sampling  tend  to  be  hetrogeneous ,  and  care  should 
be  taken  to  assure  representative  samples  and  their  interpretation. 
Seasonal  variations  in  concentration  occur  within  plant  tissue  in 
response  to  nutritional  requirements  of  the  plants.  Different  vege¬ 
tative  species  have  selective  absorption  capacities  for  intake  for 
different  elements  which  may  vary  considerably.  Lysimeter  studies 
conducted  by  Aerojet  Ordnance  and  Manufacturing  Company,  Downey, 
California,  have  indicated  that  DIMP  will  bioaccumulate  at  much 
higher  concentrations  in  the  leaves  of  plants  than  either  plant 
roots  or  stems  (see  Figure  ).  The  leaves  account  for  the  greatest 
proportion  of  organic  matter  addition  to  the  soils  which  results  in 
a  concentration  of  pollutants  in  humus  in  the  A1  horizon.  The  infor¬ 
mation  gathered  from  geobotanical  surveys  will  be  used  in  conjunction 
with  soil  and  water  surveys. 


SOILS  OF  THE  ROCKY  MOUNTAIN  ARSENAL 

\ 

There  are  eight  soil  series  and  four  land  types  at  RMA.  A 
soil  series  is  a  subdivision  of  soil  families  --  the  most  generalized 
soil  description.  A  soil  series  is  a  "group  of  soils  developed  from 
the  same  kind  of  parental  material  by  the  same  genetic  combination  of 
processes  and  whose  horizons  are  quite  similar  in  their  arrangement 
and  general  characteristics."  Soils  series  are  divided  into  specific 
mapping  units. 

ASCALON  SERIES  -  Consists  of  well  drained  soils  occurring  on 
level  to  moderately  sloping  land  on  uplands.  They  are  formed  in 
loamy  material  containing  varying  amounts  of  sand  and  gravel.  The 
infiltration  capacity  of  Ascalon  soils  is  moderate  to  high,  and  the 
available  soil  moisture  is  high.  This  soil  is  suitable  for  plant 
root  development. 

Representative  Profile  -  The  surface  layer  is  brown  sandy  loam 
approximately  six  inches  thick  and  noncalcareous.  The  upper  part  of 
the  subsoil  is  brown  sandy  loam  and  sandy  clay  loam  about  15  inches 
thick.  It  is  noncalcareous.  The  lower  part  of  the  subsoil  is  brown 
sandy  loam  about  six  inches  thick,  highly  calcareous,  and  contains 
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much  visible  lime.  The  underlying  material  is  a  highly  calcareous, 
pale-brown  to  very  pale  brown  fine  sandy  loam.  It  is  found  at  a 
depth  of  27  to  60  inches  or  more. 

Mapping  Uni ts 

(AsB)  Ascalon  sandy  loam,  one  to  three  percent  slopes.  This 
unit  is  found  on  upland  ridges  and  benches.  Surface  run-off  is 
slow  to  medium,  and  the  hazard  of  water  erosion  is  slight  to  moder¬ 
ate.  Soil  blowing  is  a  severe  hazard  in  unprotected  areas. 

(AsC)  Ascalon  sandy  loam,  three  to  five  percent  slopes.  This 
unit  is  found  in  areas  adjacent  to  major  streams,  but  it  is  only  on 
the  higher  upland  ridges  and  side  slopes.  It  has  a  thinner  surface 
layer  and  subsoil  than  the  representative  profile.  Surface  run-off 
is  medium.  The  hazard  of  water  erosion  is  moderate.  Soil  blowing 
is  severe  in  unprotected  fields. 

(AsD)  Ascalon  sandy  loam,  five  to  nine  percent  slopes.  Located 
on  uplands  near  the  edges  of  the  breaks  adjacent  to  stream  channels. 
The  surface  layer  and  subsoil  are  thinner  than  the  representative 
prof i le. 

(ArC)  Ascalon-vona  sandy  loams,  one  to  five  percent  plagues. 
Ascalon  sandy  loam,  one  to  three  percent  slopes  (AsB)  makes  up  45 
percent  of  this  complex  and  occupies  the  smoother  and  lower  lying 
areas.  Vona  sandy  loam,  three  to  five  percent  slopes  make  up  35 
percent  and  is  normally  in  the  more  sloping  part  of  the  landscape. 

A  f ew  small  flowouts  are  on  the  higher  ridge  points. 

BLAKELAND  SERIES  -  Consists  of  somewhat  excessively  drained 
soils  occurring  on  gently  sloping  to  moderately  sloping  on  uplands. 
These  soils  formed  in  wind  and  water  worked  sandy  material.  Blake- 
land  soils  absorb  water  rapidly,  and  the  available  water  capacity 
is  low.  The  entire  profile  becomes  unusually  hard  as  it  dries. 

When  dry,  these  soils  hold  a  nearly  vertical  bank.  Permeability 
is  rapid.  The  entire  soil  is  suitable  for  plant  roots. 

Representative  Profile  -  The  surface  layer  is  grayish-brown 
loamy  sand  about  five  inches  thick  and  noncal careous.  The  next 
layer  is  brown  sandy  loam  about  four  inches  thick  and  is  noncal- 
careous.  The  underlying  material  from  a  depth  of  about  9  inches 
to  60  inches  or  more  is  noncal careous ,  brown  loamy  sand  and  pale- 
brown  sand. 
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(BoD)  Blakeland  loamy  sand  three  to  nine  percent  slopes.  Occurs 
on  uplands.  Surface  run-off  is  medium  to  slow.  The  hazard  of  soil 
blowing  is  severe. 

( B t )  Blakeland  -  truck  ton  association  -  Blakeland  loamy  sand, 
three  to  nine  percent  slopes,  constitutes  at  least  60  percent  of 
this  association,  occurring  in  areas  of  semidune  relief,  though  it 
is  also  on  long,:  narrow  ridges.  Truck  ton  loamy  sand,  three  to 
nine  percent  slopes  makes  up  20  -  kO  percent  of  mapped  areas.  It 
is  depressional  or  has  the  more  gentle  slopes.  The  hazard  of  soil 
blowing  is  moderate  but  is  very  severe  where  adequate  cover  is  not 
maintained. 

(Gr)  Gravelly  Land  -  Shale  Outcrop  Complex  -  this  complex  con¬ 
sists  of  steep  and  hilly  land  types.  Areas  are  generally  elongated 
and  roughly  parallel  river  channels.  The  complex  consists  of  very 
shallow,  clayey  soils  and  raw  shale  overlain  by  a  discontinuous 
layer  of  gravel  deposits  four  inches  to  many  feet  thick.  These 
deposits  are  old  enough  that  some  soil  is  forming  in  places  within 
the  gravelly  cap  layer.  Little  if  any  soil  is  forming  in  the  shaly 
exposed  areas,  which  is  eroding  by  moderate  wind  and  water  action. 
Some  areas  contain  steep  colluvium  of  mixed  materials.  All  of  this 
complex  is  subject  to  severe  erosion. 

Loamy  Alluvial"  0nly  one  un‘t  l°amy  alluvial  is  present  at 

RMA. 


(Lu)  Loamy  alluvial  land  is  found  in  dra inageways .  Slopes  are 
less  than  three  percent.  It  is  subject  (in  varying  degrees)  to 
flooding  from  adjacent  slopes  and  main  stream  channels. 

The  profile  consists  of  a  surface  layer,  commonly  dark-colored 
loam  or  clay  loam  about  six  to  ten  inches  thick  and  normally  non- 
calcareous.  The  underlying  material  is  stratified  loam,  silt  loam, 
and  clay  loam  containing  varying  amounts  of  fine  sand,  sand,  and 
fine  gravel.  It  is  20  to  ifiore  than  60  inches  thick,  brown  to  dark 
brown  in  color,  and  in  places  it  has  a  weak  structure.  It  is  nor¬ 
mally  calcareous.  Sand  or  sand  and  gravel,  stratified  with  their 
lenses  of  silt  and  loam,  are  in  some  places  below  a  depth  of  36 
i nches. 
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Water  is  absorbed  at  a  moderate  rate,  and  their  available  water 
capacity  is  high.  They  are  normally  well  drained.  Natural  fertility 
is  high.  The  hazard  of  water  erosion  is  severe  in  all  cultivated 
areas. 

NUNN  SERIES  -  Consists  of  wel 1  drained,  nearly  level  soils  on 
terraces.  These  soils  formed  in  loamy  alluvial  material.  Water  is 
absorbed  at  a  moderate  rate,  and  the  available  water  capacity  is 
high.  Permeability  is  slow  in  soils  with  a  clay  loam  surface  layer 
and  moderately  slow  in  those  with  a  loam  surface  layer.  The  entire 
soil  is  suitable  for  plant  roots. 

Representative  Profile  -  The  surface  layer  is  a  grayish-brown 
loam  about  six  inches  thick  and  noncal careous.  The  subsoil  is  dark 
grayish-brown  and  grayish-brown  clay  about  17  inches  thick.  It  is 
noncal careous.  The  underlying  material  below  a  depth  of  23  inches 
is  light  yellowish-brown  and  grayish-brown  loam  and  silt  loam.  It 
is  highly  calcareous  and  stratified. 

(NIA)  Nunn  loam,  zero  to  one  percent  slopes.  This  soil  is  on 
stream  terraces  along  most  of  the  major  drainageways  in  the  county. 
Surface  run-off  is  slow,  and  the  hazard  of  water  erosion  is  slight. 

The  hazard  of  soil  blowing  is  moderate  in  cultivated  fields  without 
protective  cover. 

The  profile  of  this  unit  is  similar  to  the  reprensentat i ve  pro¬ 
file,  except  for  a  thicker  surface  layer  and  subsoil  and  less  sus¬ 
ceptible  to  erosion. 

(NuB)  Nunn  clay  loam,  one  to  three  percent  slopes.  This  soil  is 
found  on  stream  terraces  along  major  streams.  Areas  are  long  and  nar¬ 
row  and  roughly  parallel  to  stream  channels.  Run-off  is  medium,  and 
the  hazard  of  water  erosion  is  moderate.  Soil  blowing  can  be  control¬ 
led  by  adequate  cover.  The  profile  is  similar  to  the  representative 
profile,  but  it  has  a  more  clayey  surface  layer  and  a  thinner  subsoil. 

PLATNER  SERIES  -  Consists  of  well  drained  nearly  level  to  gently 
sloping  soils  on  uplands.  These  soils  are  formed  in  old  alluvium. 

They  absorb  water  slowly,  and  the  available  water  capacity  is  high. 
Permeability  is  slow,  and  the  entire  soil  is  suitable  for  plant  roots. 

Representative  Profile  -  The  surface  layer  is  grayish-brown 
loam  about  nine  inches  thick  and  nonca Icareous.  The  upper  part  of 
the  subsoil  is  brown  clay,  about  nine  inches  thick.  It  is  noncal- 
careous.  The  lower  part  of  the  subsoil  is  light-gray  clay  loam 
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about  10  inches  thick.  It  is  highly  calcareous,  and  part  of  the  lime 
is  visible  as  splotches.  The  underlying  material  below  a  depth  of 
28  inches  is  very  pale  brown  loam  that  is  highly  calcareous.  At  a 
depth  of  about  49  inches,  it  is  a  white  and  very  pale  brown  sandy 
loam  that  is  highly  calcareous  and  contains  some  fine  gravel. 

(P1B)  Planter  loam,  zero  to  three  percent  slopes.  This  unit  is 
found  on  uplands.  Surface  run-off  is  slow,  and  the  hazard  of  water 
erosion  is  moderate  to  slight.  The  hazard  of  soil  blowing  is  severe 
"  in  dry  farmed  areas  unless  protected  by  plants  or  stubble. 

(Smj  Sandy  Alluvial  Land  "  Consists  of  an  unstable  accumulation 
of  gravelly  and  sandy  alluvium.  It  is  in  and  adjacent  to  beds  of 
intermittent  streams.  During  periods  of  heavy  rain,  the  streambeds 
are  subject  to  flooding  and  channels  are  relocated  and  sediment  is 
shifted  and  redeposited  at  slightly  different  locations.  During  dry 
spells,  this  land  type  is  extremely  droughty. 

The  material  of  this  soil  was  transported  by  water  from  the 
sand  and  gravel  beds  in  or  adjacent  to  the  area.  It  is  stratified 
because  of  periodic  flooding.  Thin  lenses  or  small  pockets  of  silt, 
clay,  and  sand  are  also  mixed  with  the  gravel. 

This  land  type  differs  from  wet  alluvial  land  primarily  in  that 
it  is  coarser  textured  throughout  and  is  not  affected  by  the  water 
table. 

STONEHAM  SERIES  -  Consists  of  wel 1  drained,  nearly  level  to 
moderately  sloping  soils  on  uplands.  These  soils  formed  in  loamy, 
old  alluvial  material.  Water  is  absorbed  at  a  moderate  rate,  and 
the  available  water  capacity  is  moderate.  Permeability  is  moderate. - 
The  entire  soil  is  suitable  for  roots. 

Representative  Profile  -  The  surface  layer  is  light  grayish- 
brown  loam  about  five  inches  thick  and  is  noncal careous.  The  sub¬ 
soil  is  brown  sandy  clay  loam  about  eight  inches  thich  and  noncal- 
careous.  The  underlying  material  is  light-gray  sandy  loam  and 
gravelly  loam  that  is  highly  calcareous.  At  a  depth  of  30  inches, 
it  is  pale-brown  gravelly  sandy  loam  that  is  highly  calcareous. 

It  extends  to  a  depth  of  about  60  inches. 
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( St D)  Stoneham  loam,  three  to  nine  percent  slopes,  is  found  on 
uplands  near  major  stream  channels  in  scattered  areas.  Run-off  is 
generally  medium  (but  rapid  if  the  soil  is  unprotected).  The  hazards 
of  water  erosion  and  soil  blowing  are  severe  in  cultivated  areas. 

TRUCK-TON  SERIES  -  Consists  of  well  drained  to  somewhat  exces¬ 
sively  drained,  nearly  level  to  moderately  sloping  soils  on  uplands. 
These  soils  formed  in  wind-worked  sandy  material.  They  absorb  water 
rapidly,  and  the  available  water  capacity  is  low.  Permeability  is 
rapid  and  the  entire  soil  is  suitable  for  roots. 

Representative  Profile  -  The  surface  layer  is  grayish-brown 
and  dark  grayish-brown  loamy  sand  about  nine  inches  thick.  It  is 
nonca 1 careous .  The  upper  part  of  the  subsoil  is  dark-brown  and 
brown  sandy  loam  about  12  inches  thick.  It  is  noncalcareous.  The 
lower  part  of  the  subsoil  is  yellowish-brown  loamy  sand  about  11 
inches  thick;  also  noncalcareous.  The  underlying  material  is  yellow¬ 
ish-brown  coarse  sand  that  extends  to  a  depth  of  60  inches  or  more. 

(TeD)  Truck-ton  loamy  sand,  three  to  nine  percent  slopes. 

Found  mostly  along  major  stream  channels.  Surface  run-off  is  medium, 
and  the  hazards  of  soil  blowing  and  water  erosion  are  severe. 

(TuB)  Truck-ton  sandy  loam,  one  to  three  percent  slopes.  The 
profile  is  similar  to  the  representative  profile  except  for  a  sandy 
surface  layer  and  a  thicker  surface  layer  and  subsoil  surface  run¬ 
off  is  slow  to  medium,  and  the  hazard  of  water  erosion  is  slight. 

The  hazard  of  soil  blowing  is  severe. 

(TuD)  Truck-ton  sandy  loam,  two  to  nine  percent  slopes.  The 
profile  is  similar  to  the  representative  profile  except  for  a  sandy 
loam  surface  layer.  Run-off  is  medium,  and  the  hazard  of  water  ero¬ 
sion  is  moderate.  The  hazard  of  soil  blowing  is  severe. 

VONA  SERIES  -  Consists  of  well  drained,  nearly  level  to  moder¬ 
ately  sloping  soils  on  uplands.  They  were  formed  in  wind-deposited 
material.  The  soils  absorb  water  rapidly,  and  the  available  water 
capacity  is  moderate.  Permeability  is  rapid,  and  the  entire  soil  is 
suitable  for  plant  roots. 

Representative  Profile  -  The  surface  layer  is  light  brownish- 
gray  loamy  sand  about  nine  inches  thick.  It  is  noncalcareous.  The 
subsoil  is  pale-brown  coarse  sandy  loam  about  13  inches  thick;  also 
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is  noncalcareous.  The  underlying  material  is  light  yellowish-brown 
sandy  loam  and  very  pale  brown,  highly  calcareous  loamy  sand  that 
contains  lime  visible  as  splotches. 

(VnD)  Vona  loamy  sand,  three  to  nine  percent  slopes.  Surface 

run-off  is  medium,  and  the  hazard  of  soil  blowing  is  severe. 

(VoC)  Vona  sandy  loam,  three  to  five  percent  slopes.  The  pro¬ 

file  is  similar  to  the  representative  profile  except  for  a  sandy 
loam  surface  layer  and  a  thinner  subsoil.  The  hazards  of  water 
erosion  and  soil  blowing  are  severe. 

WELD  SERIES  -  Consists  of  well  drained,  nearly  level  soil  on 
uplands.  The  soils  were  formed  in  wind-worked  loamy  materials.  The 
soils  absorb  water  at  a  moderate  rate,  and  the  available  water  capacity 
is  high.  Permeability  is  slow,  and  the  entire  soil  is  suitable  for 
plant  roots. 

Representative  Profile  -  The  surface  layer  is  brown  loam  about 
six  inches  thick.  It  is  noncalcareous.  The  upper  part  of  the  sub¬ 
soil  is  dark  brown  clay  about  six  inches  thick.  It  is  noncalcareous. 
The  lower  part  of  the  subsoil  is  pale-brown  and  very  fine  sandy  loam 
about  20  inches  thick.  It  is  highly  calcareous,  and  much  of  the  lime 
is  visible  as  splotches.  The  underlying  material  is  highly  calcar¬ 
eous,  pale-brown  silt  loam  and  light  yellowish-brown  fine  sandy  loam. 

(WmB)  Weld  loam,  one  to  three  percent  slopes.  This  soil  is 
found  on  uplands.  Surface  run-off  is  medium,  and  the  hazard  of 
water  erosion  is  moderate  to  severe.  The  hazard  of  soil  blowing  is 
severe  in  clay  farmed  areas,  if  rainfall  is  below  normal. 

(Wt)  We t  Alluvial  Land  “  This  is  found  on  the  nearly  level  bot¬ 
tom  lands  of  the  larger  streams  next  to  stream  channels.  The  land 
is  wet  most  of  the  time  and  is  flooded  by  streamflow  once  to  several 
times  a  year  during  periods  of  high  water. 

Materials  are  extremely  variable  in  texture,  consisting  of 
stratified  layers  of  dark-colored  silt,  loam,  and  clay.  The  layers 
are  generally  less  than  six  inches  thick  and  are  underlain  by  sand, 
fine  sand,  and  some  gravel  at  depths  of  one  to  three  feet.  They  are 
wet  at  a  depth  of  two  feet  most  of  the  time  and  are  commonly  wet  to 
the  surface  throughout  the  growing  season.  Water  tolerant  plants 
form  the  vegetative  cover. 
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SAMPLING  METHODS 


One  of  the  objectives  of  the  comprehensive  survey  is  to  appraise 
the  interna]  organization  of  contaminant  distribution;  the  location  of 
the  elements  of  the  distribution  with  respect  to  each  other.  The 
"spatial  structure"  represents  the  location  of  each  element  relative 
to  each  of  the  others  and  also  the  location  of  each  element  relative 
to  all  of  the  contaminants  taken  together;  that  is,  the  internal  loca¬ 
tional  organization  of  contaminant  distribution  in  space.  The  "spa¬ 
tial  processes,"  on  the  other  hand,  are  mechanisms  which  produce  the 
spatial  structure  of  the  contaminant  distribution. 

Spatial  structure  is  a  determinant  of  spatial  process  as  much 
as  process  is  a  determinant  of  structure.  For  example,  the  geo¬ 
hydrologics  characteristics  of  the  Arsenal  environment  control 
groundwater  movement,  which  also  affects  contaminant  migration  and 
the  spatial  pattern  of  contamination.  An  understanding  of  the  spa¬ 
tial  process  can  be  achieved  by  the  deductive  approach,  using  the 
system's  concept  which  requires  fewer  samples  then  would  be  required 
to  map  out  the  spatial  pattern. 

Contaminant  emplacement  and  migration  will  be  approached  from 
the  systems  viewpoint.  Systems  modeling  necessitates  a  clear  defi¬ 
nition  of  the  objectives  of  the  study,  to  identify  and  model  the 
components  of  the  system  andsubsystems ,  and  to  elucidate  their  mutual 
relationships,  especially  with  regard  to  the  spatial  patterns. 

System  modeling  has  the  ability  to  handle  a  dynamic  process. 

It  has  the  ability  to  look  at  the  total  ecosystem,  makes  efficient 
use  of  information  gathered,  has  the  ability  of  checking  consistency, 
has  the  ability  of  information  integration,  is  applicable  both  as  a 
planning  tool  and  an  analysis  tool,  and  is  an  excellent  means  of 
organizing  data  for  decision  makers.  The  process  response  system, 
incorporating  both  morphological  and  cascading  subsystems,  can  serve 
as  a  working  model  to  demonstrate  the  linkage  between  form  and  process. 
The  strength  and  direction  of  connectivity  between  the  components  is 
analyzed  by  statistical  methods.  Such  methods  can  range  from  analy¬ 
sis  of  variance,  correlation  and  multiple  regression,  to  cross- 
spectral  analysis. 

The  problem  in  sampling  and  the  selection  of  a  sampling  design 
for  meaningful  data  collection  is  predicated  by  the  research 
obj  ect i ves. 
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The  purpose  of  sampling  is  to  accurately  portray  the  nature  of 
the  environment  from  which  the  samples  were  taken.  These  are  two 
standard  categories  of  sampling  methods,  probability  sampling  and 
nonprobability,  or  judgment  sampling. 

1.  Probability  Sampling  -  Probability  sampling,  based  upon  random 
selection,  i's  used  in  modeling  and  forecasting  the  behavior  of 
naturally  occurring  environmental  systems.  It  requires  that  each 
element  have  an  equal  chance  of  being  chosen;  that  is,  independent 
of  any  event  in  the  selection  process.  Random  selection  procedures 
eliminate  conscious  or  unconscious  bias  by  the  scientist  and  enables 
the  researcher  to  use  probability  theory  as  a  method  for  estimates 
of  population  parameters  and  estimates  of  error. 

a.  Simple  Random  Sampling  -  After  a  sampling  frame  has  been  esta¬ 
blished,  numbers  are  assigned  to  each  of  the  elements  in  the  list. 
The  elements  to  be  measured  are  then  selected  randomly  from  the  list. 
If  spatial  patterns  are  to  be  sampled,  then  a  grid  overlay  is  used 
and  the  X  and  Y  coordinates  are  selected  randomly  to  identify  either 
the  coordinate  intersection  to  be  measured  as  a  point  sample  or  the 
appropriated  grid  cell  to  be  measured  as  an  area  samples. 

b.  Systematic  Sampling  -  The  first  element  is  selected  randomly, 
then  every  KtTT  element  in  the  total  list  is  chosen  for  inclusion  in 
the  sample.  Periodicity  and  cyclical  patterns  that  coincide  with 
the  sampling  interval  may  be  missed  and  may  bias  the  results. 

c.  Stratified' Sampi ing  -  Stratification  organizes  the  population 
into  homogeneous  subsets  (with  heterogeneity  between  subsets)  and 
an  appropriate  number  of  elements  are  measured  from  each.  Stra¬ 
tification  ensures  proper  representation  of  the  variables  describing 
the  parent  population.  An  appropriate  number  of  elements  are  drawn 
from  homogeneous  subsets  of  the  population.  This  minimizes  sampling 
error. 

d.  Multi-Stage  Cluster  Sampling  -  Cluster  sampling  involves  sampl¬ 
ing  groups  of  elements  (clusters)  followed  by  selection  of  elements 
within  the  clusters.  Multi-stage  cluster  sampling  involves  the 
repetition  of  listing  and  sampling.  As  the  homogeneity  increases 
within  each  cluster,  relatively  fewer  elements  may  be  needed  to 
adequately  represent  a  given  natural  cluster.  However,  a  larger 
number  of  clusters  may  be  needed  as  the  diversity  among  clusters 
increases. 
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:\  ,  e.  Hierarchied  Sampling  (Nested  Sampling)  -  A  region  is  first  sub- 

1  divided  into  major  areas  of  equal  size.  Several  major  areas  are  then 
I  chosen  at  random  and  partitioned  into  smaller  sub  regions.  The  sub 
regions  are  randomly  selected  and  further  subdivided.  Hierarchied 
sampling  ensures  that  every  part  of  the  target  region  is  represented. 

f.  Disproportionate  Sampling  and  Weighing  -  A  probability  sample 
is  one  in  which  each  population  element  has  a  known  nonzero  prob¬ 
ability  of  selection,  even  though  different  elements  may  have  dif¬ 
ferent  probabilities.  If  controlled,  probability  sampling  procedures 
can  be  used  and  may  be  representative  of  the  population  from  which 
it  was  drawn,  if  each  sample  element  is  assigned  a  weight  equal  to 
the  inverse  of  its  probability  of  selection. 

2 .  Nonprobability  Sampling;  Purposive  or  judgment  sampling  - -  a 
researcher  may  sample  with  prior  knowledge  of  the  population,  its 
elements  and  research  objectives.  This  is  useful  to  pretest  a 
sample  plan  or  to  select  a  typical  or  representative  sample. 

b,  Quota  Sampling  -  Quota  sampling  begins  with  a  matrix  describ¬ 
ing  the  characteristics  of  the  target  population.  The  ration  of 
the  population  occurring  in  various  categories  is  known.  A  pre¬ 
determined  proportion  is  selected  for  sampling  from  each  cell  of 
the  matrix. 

c.  Traverses  -  Traverse  sampling  is  the  technique  used  to  sample 
the  environment  along  a  straight  line  whose  length  and  direction 
have  been  determined.  Traverse  sampling  may  be  used  to  gather  data 
without  disturbing  the  environment.  The  assumption  that  the  sum  of 
the  intercepts  along  a  traverse  line  can  be  used  to  estimate  the 
proportion  of  the  total  area  occupied  by  the  phenomenon  being  inter¬ 
cepted.  The  intercept  lengths  reflect  the  spatial  geometry  and  the 
degree  of  fragmentation  of  the  occurrences  of  the  attribute  being 
sampled  in  the  region. 

Probability  sampling  avoids  conscious  or  unconscious  bias  in 
element  selection  which  closely  represents  the  population  of  all 
elements.  Probability  sampling  permits  estimates  of  sampling  error. 
The  selection  of  a  particular  sampling  design  depends  upon  the  accur¬ 
acy  and  precision  required;  the  efficiency  and  feasibility  of  a  par¬ 
ticular  procedure;  and  the  cost  in  time,  manpower,  and  money.  The 
design  selected  also  depends  upon  the  objectives  of  the  investigation, 
the  structure  or  functional  characteristics  of  the  population. 
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A  combination  of  the  purposive  sampling  approach  (which  allows 
intensive  study  of  structure  and  relationships  and  therefore  hypo¬ 
thesis  formulation)  with  a  probability  sampling  approach  (which 
allows  the  hypotheses  to  be  formally  tested  with  respect  to  a  given 
population)  is  advocated  for  environmental  studies. 

Estimating  the  Sample  Size 


Sampling  involves  a  compromise  between  precision  and  economy. 

A  parameter  can  only  be  defined  with  100  percent  accuracy  by  measur¬ 
ing  or  counting  the  whole  population.  The  smaller  the  sample,  the 
wider  the  margin  of  error  of  any  estimate  made  from  it.  A  homo¬ 
geneous  population  produces  samples  with  smaller  sampling  errors 
than  does  a  heterogeneous  population.  If  the  margin  of  error  that 
can  be  tolerated  for  a  particular  purpose  is  known,  it  is  possible 
to  calculate  the  minimum  sample  size  required  and  thereby  same  time 
in  the  collection  of  an  unnecessarily  large  sample.  Different  methods 
in  environmental  studies  can  be  used  depending  on  whether  the  samples 
are  to  be  measured  or  counted.  Both  involve  prior  collection  of  a 
pilot  sample  of  about  30  observations.  These  samples  must  be  collected 
by  the  same  methods  that  are  to  be  employed  in  the  full  sample;  and 
they  can  be  used  as  part  of  the  full  sample,  provided  they  do  not 
affect  the  selection  of  additional  items. 

The  formular  for  finding  the  required  sample  size  using  measured 
data  is:  n  =  (24) 2 ,  where  £  is  the  required  size  of  the  sample,  4 
is  the  standard  deviation  of  the  pilot  sample,  d_  is  the  tolerable 
margin  of  error  at  a  specified  level  of  confidence,  and  £  is  a  con¬ 
stant  taken  from  the  z-table  (table  of  probabilities  associated  with 
values  of  z  in  a  normal  distribution)  corresponding  to  the  same  level 
of  confidence.  Therefore,  if  the  pilot  sample  has  a  standard  devia¬ 
tion  of  10  and  if  the  margin  of  error  that  can  be  tolerated  is  2.5  at 
the  95  percent  level,  then  the  required  sample  size  is  n  =  (2x10/2. 5)^ 
64. 


The  formula  for  finding  the  sample  size  using  counted  data  is 
n  =  p%  x  q%  x  (z/d)^  percent,  where  n_,  z_,  and  d_  have  the  same  mean¬ 
ings  as  in  paragraph  above  and  where  £  and  £  are  the  percentages  of 
the  pilot  sample  belonging  and  not  belonging  to  the  specified  category 
respect  ?  vel y. 


17 


geographic  coordinate  system  used  at  rma 


Several  different  coordinate  systems  have  been  used  in  the  past 
on  RMA.  Each  coordinate  system  has  its  own  merits.  However,  there 
exists  a  need  to  coordinate  all  of  the  data-gather ing  tasks  into  one 
standardized  reference  system.  The  identification  code  of  each  sample 
collected  wil,.l  be  located  with  reference  to  the  State  Plane  Coordinate 
System,  and  all  existing  data  will  be  converted  to  this  system. 

1  •  State  Plane  Coordinates  represent  a  system  set  up  for  each 
state  by  the  U.  S.  Department  of  Commerce,  Coast  and  Geodetic  Survey, 
to  be  used  by  the  private  land  surveyors  of  those  states.  The  system 
consists  of  a  rectangular  coordinate  grid  that  covers  the  entire  state, 
which  is  oriented  due  north  and  is  reduced  to  sea  level.  These  coord¬ 
inates  can  be  tied  into  a  more  universal  Geographic  Coordinate,  which 
utilizes  latitude  and  longitude.  Any  point  to  be  referenced  within 
the  state  can  be  identified  on  the  State  Plane  Coordinate  Grid.  This 
has  the  advantage  of  being  more  flexible  and  universal  when  tracing  data 
that  extends  beyond  the  Arsenal  boundaries. 

2.  RMA  Coordinate  System  -  This  system  was  established  in  1 942 
on  the  original  survey  of  the  Arsenal.  The  origin  of  the  system  was 
a  Geological  Survey  monument  named  Sand  1936  (now  destroyed).  The 
grid  is  oriented  to  magnetic  north.  Out  of  about  160  coordinate  points 
established  in  the  past,  only  eight  points  have  been  recorded.  This 
system  is  completely  independent  of  all  other  coordinate  systems  used 
in  the  area.  The  origin  of  this  grid  system  is  located  in  the  southwest 
corner,  such  that  all  X  and  Y  values  on  the  Arsenal  are  positive  numbers. 
Many  of  the  original  Arsenal  design  plans  use  the  RMA  Coordinate  System. 
This  system  is  inadequate  for  !R  needs  because  negative  coordinates 
would  occur  for  locations  beyond  the  Arsenal. 

3*  The  U.  S.  Rectangular  System  is  based  on  a  grid  consisting  of 
six  mile  wide  east-west  tiers  and  north-south  ranges.  Each  grid  cell 

is  a  township  comprised  of  36  one  square  mile  sections  which  are  num¬ 

erically  numbered  in  snake  fashion  beginning  in  the  northeast  corner. 

Each  one  mile  block  of  land  or  section  has  its  own  grid  system  and  is 
therefore  difficult  to  use. 

The  Sect  ion  Grid  is  a  coordinate  system  set  up  for  each  one 
mile  block  within  a  section  of  land.  The  southwest  corner  of  the 
section  is  arbitrarily  set  at  1,000  feet  north  and  1,00  feet  east. 

This  small  grid  is  used  with  very  little  control  surveying.  To  set 

up  this  grid,  only  the  boundaries  of  the  section  need  to  be  measured. 

A  control  survey  is  then  completed  to  locate  the  section  corners  on 
another,  more  universal,  system.  The  section  coordinates  can  then 
be  converted  to  the  control  survey  grid.  This  is  a  localized  grid 
system  which  is  easy  to  use  for  specific  tasks. 
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PRINCIPAL  COMPONENT  FACTOR  ANALYSIS 


Factor  analysis  is  a  statistical  technique  which  is  used  to 
reduce  a  square  symmetric  (mxm)  matrix  to  a  related  rectangular 
(mxp)  matrix  such  that  m>p.  The  initial  matrix  is  a  correlation 
coefficient  matrix  consisting  of  m  variables.  The  output  matrix 
presents  restructuring  the  pattern  of  correlations  into  p  clusters 
of  characteristics,  accounting  for  much  of  the  variance  of  the  m 
set  of  initial  variables. 

The  p  set  of  final  vectors  then  constitutes,  a  set  of  surrogate 
variables  (called  factors),  which  represent  the  m  set  of  variables. 

Within  IR,  we  will  deal  exclusively  with  principal  component 
factor  analysis  which  selects  an  orthogenal  rather  than  an  oblique 
(correlated  factor)  solution  and  which  assigns  all  of  the  variance 
for  interpretation  rather  than  only  part  of  it. 

Principal  components  method  is  preferred  because  of  ortho¬ 
gonality  and  the  uniqueness  of  the  solution.  There  is  only  one 
principal  components  solution  to  a  correlation  matrix,  thus  the 
factor  structure  is  not  likely  to  be  an  artifact  of  procedure.  The 
factors  in  the  solution  are  uncorrelated  and  are  statistically 
independent.  Therefore,  surrogate  variables  are  independent  by 
design.  If  there  are  N  observation  drill  sites  in  each  cluster 
and  each  site  provides  information  with  respect  to  m  variables, 
then  at  each  point  in  time  a  data  matrix  (N  x  m)  is  obtained. 
Simulation  of  time  trends  may  result  in  error  because  the  signifi¬ 
cant  correlations  among  the  variables  are  likely  to  be  interdependent. 
However,  the  correlation  structure  among  the  m  variables  can  be  redu¬ 
ced  into  p  orhoganal  components,  with  which  them  variables  are  associ¬ 
ated  by  means  of  linear  equations.  Time  trends  of  the  p  components 
may  now  be  analyzed  without  fear  of  independence  problems.  The' 
results  of  the  simulations  representing  p  components  can  be  recon¬ 
structed  back  to  the  original  variables.  For  example,  the  following 
flow  diagram  illustrates  the  steps  in  factorial  simulation: 
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Principal  component  analysis  can  be  used  in  problems  which 
involve  an  analysis  of  spatial  structure  of  contamination  and  spa¬ 
tial  taxonomy.  For  example,  an  area  of  interest  such  as  the  northern 
boundary  is  subdivided  into  N  subareas  or  clusters.  For  each  subarea, 
information  on__m  descriptive  variables  is  measured.  This  constitutes 
a  data  matrix  X  (N  x  m) .  The  questions  we  now  pose  are: 

1.  How  may!  the  internal  variations  of  the  m  variables  over  the 
entire  area  be  parsimoniously  described? 

2.  How  may  the  subareas  be  arranged  to  reveal  the  similarity 
and  difference  so  as  to  assign  a  meaningful  regional  structure? 

A  data  matrix  consisting  of  variables  measured  at  strategically 
located  sites  positioned  to  describe  the  environment  can  be  subjected 
to  correlation  analysis  to  reveal  how  the  descriptive  variables  are 
intercorrel ated  and  therefore  the  degree  to  which  they  report  redun¬ 
dant  information. 

Principal  components  analysis  can  be  used.  The  oringal  variables 
may  be  correlated  and  numerous.  The  patterns  of  variation  may  be  inde¬ 
pendent  and  confusing.  These  can  be  reduced  to  p  factors  which,  on  the 
other  hand,  are  few  and  independent  and  thus  describe  the  region  more 
parsimoniously.  Furthermore,  there  is  no  way  of  knowing  which  of  the 
variables  are  the  most  powerful  descriptive  tools;  whereas,  the  princi¬ 
pal  components  method  ranks  the  factors  in  descending  order  of  ability 
to  interpret  the  correlation  matrix  variance.  If  factor -fj  interprets 
b2  percent  of  the  variance  and  factor  f£  3b  percent,  we  can  ascertain 
which  is  more  general  and  by  how  much. 

For  mapping  and  regionalizing,  it  is  desirable  to  be  able  to  esti¬ 
mate  the  strength  of  each  factor  within  each  tract.  For  this  purpose, 
the  matrix  A  (m  x  p)  which  gives  the  relationships  between  the  variables 
and  the  factors  is  combined  with  a  matrix  X  Z  (N  x  m)  of  standardized 
values  of  the  initial  variables.  The  result  is  a  matrix  of  dimension 
(N  x  p)  which  estimates  the  strength  of  each  factor  in  each  subarea  of 
observation.  This  is  called  the  factor  score  matrix  F.  The  sequence 
is  as  fol lows : 


Flow  Diagram  of  Principal  Components  Analysis 
of  a  Data  Matrix 
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Each  row  of  the  factor  score  matrix  of  a  given  sampling  bore  hole 
(or  subarea)  across  the  independent  dimensions  represented  by  the  columns. 
It  is  therefore  actually  a  point  in  the  p-d imens i onal  space  of  the  fac¬ 
tor  structure.  One  may  therefore  evaluate  the  similarity  or  difference 
between  subarea  i  and  subarea  j  by  means  of  the  pythagorean  theorem  which 
evaluates  the  distance  between  any  two  points  in  a  p-dimension  space. 


(1) 


d?j  = 


(e 


il 


-  e 


jl 


)2  +  (e 


i2 


-  e 


j  2 


)2 


+  (e 


‘P 


-  e 


JP 


Equation  (l)  may  be  evaluated  for  every  pair  of  subareas,  and  the 
set  of  { d 2 j }  will  then  form  the  B2  matrix.  Analysis  of  this  matrix  -- 
especially  by  such  techniques  as  discriminant  analysis  —  will  cluster 
those  subareas  which  most  resemble  one  another  into  families  linked  into 
a  regional  hierarchy  of  subareas.  The  grouping  of  the  {dj-.}  by  discri¬ 
minant  analysis  will  be  presented  at  a  later  date.  J 


The  Procedure  of  Principal  Components  Analysis 

If  we  consider  the  matrix  Z  (N  z  m)  of  standardized  values  of  the 
N  subareas  on  the  m  variables,  we  may  regard  each  variable  as  a  vector 
in  a  N  dimensional  space  of  orthogonal  axes.  The  columns  of  the  matrix, 
in  this  case,  form  the  vectors  which  define  each  variable's  positions 
in  this  space.  Since  we  are  dealing  with  a  standard i zed  value  matrix, 
all  vectors  are  of  equal  (and  unit)  length.  The  similarity  of  any 
variable  to  any  other  is  therefore  measured  by  the  angle  between  them. 
But,  the  cosine  of  the  angle  between  two  standardized  vectors  in  a 
space  of  N  dimensions  is  given  by  the  inner  vector  product  between 
them.  The  cosine  function  varies  between  -1  and  +1  and  is,  in  fact, 

(2)  cos  0 i j  =  zj  .  zj 

A/ 

the  correlation  coefficient  between  the  i£il  and  jtD  variable.  Thus, 
we  find  R,  the  correlation  coefficient  matrix,  by  expression  (3). 

__  _c  _ 

(3)  R  -  Z  .  Z 

R  is  a  matrix  which  measures  the  angles  between  variable  vectors. 
The  principal  components  solution  starts  by  discovering  how  orthogonal 
axes  may  be  passed  through  this  pattern  of  intersecting  vectors  so  that 
each  axis  will,  as  far  as  it  is  possible,  coincide  with  the  center  of  a 
"cone"  of  vectors  (i.e. ,  be  correlated  with  it).  This  is  analogous 
with  the  physicist's  problem  of  finding  the  resultant  of  a  set  of  force 
vectors.  In  the  F  matrix  case,  the  resultant  is  found  by  solving  the 
characteristic  equation  of  the  matrix. 
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In  ordinary  algebra,  an  expression  of  the  type 

1  W  ay  =  gy 

implies  that  a  =  g.  In  matrix  algebra,  a  common  factor  on  each  side 
of  the  equation  does  not  imply  the  equality. 

'  (5)  AY  .=*  Y 

In  expression  (5),  A  may  be  a  (m  x  m)  matrix,  and  Y  may  be  a 
(m  x  1 )  vector.  A  may  be  a  scalar,  and  it  will  then  be  true  that  the 
left  hand  term  and  right  hand  term  will  both  be  (m  x  1)  vectors;  yet 
A  ^  X..  Scalars  such  as  A  and  vectors  such  as  Y,  which  have  this 
property  with  respect  to  a  given  matrix,  A,  are  called  eigenvalues 
and  eigenvectors.  When  the  correlation  coefficient  matrix  R  is  fac¬ 
tored  into  eigenvalues  and  eigenvectors  (factor  analysi s) ,  the  eigen¬ 
values  are  orthogonal  and  have  the  properties  which  we  seek.  Further¬ 
more,  the  eigenvalues  prove  to  be  scaling  coefficients  such  that  they 
are  proportional  to  the  amount  of  variance  interpreted  by  each  factor. 
Any  factor  with  a  A>  1.0  interprets  more  variance  than  an  unmodified 
variable  alone,  while  any  factor  with  a  x< 1-0  is  of  relatively  weak 
explanatory  power.  Generally,  the  researcher  rejects  factors  with 
A<  1.0,  and  the  ZAi  of  these  rejected  factors  is  proportional  to  the 
unexplained  variance  in  the  model. 

The  factor  loading  matrix  (A)  is  synthesized  from  the  eigenvalues 
and  eigenvectors  so  that 

(5)  R  A  =  A.  A  . 

where  A  is  a  diagonal  matrix  whose  elements  on  the  principal  diagonal 
are  the  A‘s. 

A  matrix  F  can  then  be  derived  whose  columns  are  vectors  f],  f2, 
corresponding  to  the  scores  of  each  factor  in  a  certain  tract,  by  seek 
ing  some  relationship  based  upon  the  property  that 

(6)  z  =  A  f 

which  is  to  say  thaj:  z  (a  column  vector  of  the  standardized  variables) 
for  the  tract  and  f  (a  column  vector  of  the  factors)  are  linearly  rela 
ted  through  A.  Then,  premultiplying  by  A  gives 

(7)  A'z  =  AAf 

and  finally 

(8)  f  =  (A7  A)'1  A'z 
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Expression  (8)  provides  the  means  for  estimating  the  values  of 
i the  p  factors  for  all  the  tracts.  These  are  the  values  which  may  be 
(mapped. 

Rotation  of  the  factor  structure  to  any  one  of  an  infinite 
number  of  configurations  will  yield  different  resultants.  The  most 
useful  rotation  for  IR  is  that  which  uses  the  varimax  criterion. 

The  purpose  of  this  is  to  place  the  component  axes  in  the  unique 
position  where  as  much  of  each  factor  as  possible  is  interpreted  by 
•strong  loadings  by  as  few  variables  as  possible.  This  rotation  sim¬ 
plifies  factor  structure  in  that  each  factor  axis  is  placed  as  para¬ 
llel  as  possible  to  a  limited  number  of  variable  vectors. 

Rotation  serves  a  very  useful  purpose  in  clarifying  t,he  inter¬ 
pretation  of  the  model  by  a  rule  of  parsimony.  Furthermore,  the  vari¬ 
max  rotation,  oke  the  principal  components  solution,  yields  a  unique 
solution  whose  factors  are  orthogonal  --  but  with  the  added  advantage 
of  simplicity  of  interpretation. 


FACTOR  ANALYSIS  OF  CONTAMINANTS  ALONG  THE  NORTHERN  BOUNDARY 

The  following  example  will  illustrate  the  use  of  principal  com¬ 
ponent  factor  analysis  for  (l)  data  reduction,  (2)  analysis  of  spa¬ 
tial  structure,  and  (3)  spatial  mapping. 

Twenty-four  bore  hole  samples  were  drilled  250  feet  apart  in  two 
parallel  east  to  west  lines  along  the  northern  boundary  at  RMA.  Chemi¬ 
cal  analysis  was  performed  on  the  groundwater  collected  from  each  sample 
and  tested  for  eight  compounds  (eight  variables).  The  raw  data  matrix 
(2k  x  8)  is  shown  on  Table  1,  in  which  the  rows  represent  the  bore  holes 
(samples,  N) ,  and  the  columns  represent  the  variables  (m) . 

Table  3  is  an  eight  by  eight  correlation  matrix  (M  x  M) ,  which 
represents  the  product  moment  correlation  between  each  of  the  eight 
compounds  analyzed.  When  the  number  of  variables  in  the  set  is  small, 
it  is  easy  to  evaluate  those  variables  which  correlate  highly  with  each 
other.  When  the  number  of  variables  increases,  it  becomes  difficult  to 
identify  those  variables  which  are  intercorrelated  and  group  together 
as  clusters.  These  groups  can  be  isolated  as  orthogonal  factors.  Table 
6  shows  three  factors  represented  as  columns.  The  variables  in  this 
(m  x  p)  matrix  are  shown  as  the  rows  where  the  factor  loadings  and  the 
factors  represent  their  relationships. 

The  cl"  ion,  DIMP  and  the  Na  ion  load  up  on  factor  1;  therefore, 
they  must  also  be  highly  intercorrelated  and  occur  as  a  group  in  this 
data  set. 
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SAMPLES 


RAW  DATA  MATRIX 


VARIABLES 


A-l, 
A-Z 
A- 3 
A~4 
A- 5 
A- 6 


PPM7  PPM  PPM  rP PM  9PM 
Vr"  V 2  v 3  : 

[■44 54120 4 550 .206407 .  6222000031* 
i'27630 001950 . 04:3207.860940  0290 
'33530002100,043857.4807301600 
0750750  0850. 091 207. 900  0920073 
.230:26304151 . 404707.5904000000 
i22027306552, 956507.8001700000 


A- 7 

A- 3 

D-1 

D-Z 

D-3 

D-4 

D-5 

D-6 

D-7 

D-8 

UD-1 

UD-2 

UD-i} 

UD-5 

UD-6 

UD-7 

UD-8 

UD-3 


j50034307500. 127257.4015000000 
•20037309552.557057.9100820000 
'28051302590 . 043507.7712000290 
214:20801520 . 062757.9105300240 
095,07501300. 061808.0100660000 
24526805503. 006407.6702700000 
[52552505300 . 186507.5529000000 
36527002000.044757.6507501200 
24519501800 .213407. 7004800088 
225'32508552, 008557. 6701500010 
24221402000. 144207 .6005700320 
124023403250 .045007.6504200010 
44 561501 550 . 0574 0;7 . 49*190  0  0  000 
132527003950 . 044357. 7407000093 
|170345:1 0  003. 158507.6200630  000 
128529002100  •  06295|7 . 980960  0093 
^04200030 02. 004658. 2603400  012 
jl 6 92  5 30  6 1 53 . 3 0l66 5j7 . 520  2  4  0  0  0  0  0 


TABLE  1 


VARIABLE 


cl~  1 
Ha  2 
SOif  3 
NO2/NO3  4 
Hardness  5 
pH  6 
DIMP  7 
DCPD  8 


MEAN 


273.125011 

289.125000 

407.333340 

.906667 

506.249996 

7.722917 

735.541687 

181.250010 


ST. DEV. 


115.321591 
117.819796 
273.620590 
1.241141 
203.855494  . 

.200531  • 
735.724869 
394.055733 


TABLE  2 
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000005  ,947572 

000005  .889357 


SAMPLES 


ROTATED  FACTOR  MATRIX 


VARIABLE 

1 

FACTOR 

2 

3 

■■  A  A  0*7  Q 

cl'  -1 

Na  ^ 

S°  if  f 

N02/N03  4 

Hardness  5 

pH  6 

DIMP  7 

DCPD  8 

.95432 
.86626 
.07578 
-.36283  1 

.42726 
-.59616 
'.93632 
.02320 

.05771 

-.32869 

-.93095 

-.86278 

-.87938 

.32518 

.23634 

.31009 

.119/0 

-.04775 

-.09502 

-.10680 

-.00878 

-.52109 

-.12260 

.89032 

TABLE  6 


A-1 
A- 2 
A- 3 
A- 4 
A- 5 
A-6 
A-7 

A- 8 
D-l 
D-2 

D-3 

D-4 

D-5 

D-6 

D-7 

D-8 

UD-1 

UD-2 

UD-4 

UD-5 

UD-6 

UD-7 

UO-8 

UD-3 


CASE 
~~  1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


FACTOR  SCORES 

Factor  1 


Factor  2 


- 1763071 

. 04633 
.13449 
-1.46349 
-.28705 
-.63546 
1.47940 
-.50271 
.27072 
-.55679 
-1.44698 
-.46282 
2.46890 
. 16364 
-.35982 
-.22932 
-.21799 
-.17469 
2. 12827 
..15308 
-.48184 
.00252 
-1.01234 
]  -.64515 


. 14408 
.92987 
.57079 
1.25532 
-.26986 
-1. 16170 
-.57263 
-1.52047 
.78231 
1.02852 
1.19991 
-1.09225 
.09476 
.59346 
.69941 
-1.60303 
.51007 
.23713 
.10989 
.39788 
-2.18231 
1.04124 
.24606 
-1.44025 


Factor  3 
-.63347 
-.24377 
3.43871 
-.37236 
.07464 
-.26449 
.13603 
-.57213 
-.12060 
-.26528 
-.75781 
.01025 
i  -.85149 
2.25481 
-.03682 
.04838 
.64720 
-.05952 
-.45420 
-.23061 
.18710 
-.93460 
-1.39069 
.39619 


'Var 


1  *■  c  1 

2  =  Na 

7  =  D1MP 


|var  3  =  SO^  1 

4  =  NO2/NO3 

5  =  Hardness 


8  =  DCPD 


TABLE  7 


500  UeT 


Factor  2  represents  the  S0]J  ion,  NC^/NO-j,  and  hardness  variable 
group.  DCPD  behaves  independently  represented  as  Factor  3  in  this 
data  set.  It  is  possible  from  the  factor  loading  to  identify  which 
variable  (model  variable)  has  the  highest  correlation  with  the  respec¬ 
tive  factor.  This  model  variable  can  then  represent  the  group  in 
future  tests  and  thus  permit  data  reduction  with  minimum  information 
loss. 

Table  4  gives  the  eigenvalues  which  can  be  interpreted  as  the 
.proportion  of  the  variance  in  the  data  matrix  which  is  explained  by 
each  of  the  factors.  By  weighting  the  proportion  each  eigenvalue  is 
of  the  total,  we  can  determine  the  relative  explanatory  power  of  each 
of  the  factors  in  terms  of  the  total  common  variance  in  the  data  set. 
Factor  1  accounts  for  42  percent  of  the  total  variance,  while  Factor 
2  explains  33  percent.  All  three  factors  together  explain  88  percent 
of  the  total  variance  in  this  example. 

Table  7  gives  the  factor  scores.  These  numbers  represent  the 
weight  of  each  of  the  factors  at  each  of  the  24  drill  sites.  By  map¬ 
ping  and  contouring  the  factor  scores  for  each  of  the  factors  separa¬ 
tely,  we  can  construct  maps  showing  multivariate  distribution  repre¬ 
senting  the  spatial  behavior  of  groups  of  contaminants  which  behave 
spatially  together.  See  map 

This  information  is  invaluable  in  (l)  tracing  each  of  the  con¬ 
taminant  migration  paths  independently,  (2)  determining  the  source 
of  contamination,  and  (3)  determining  the  position  and  direction  for 
further  sampl i ng. 

SAMPLING  SITE  CHARACTERIZATION  &  DRILL  CORE  LOGGING 

OVERVIEW 

Soils  are  organized  natural  bodies  which  vary  markedly,  both 
horizontally  across  the  surface  of  land  and  vertically  as  seen  in 
cross  section.  The  vertical  section  provides  the  best  means  of 
examining  soil.  The  component  layers  that  form  the  basic  unit  for 
the  classification  of  soils,  in  turn,  provides  a  basic  framework 
for  the  analysis  and  mapping  of  the  behavior  of  the  pollutants  in 
the  ground. 

The  profile  form  which  expresses  the  overall  visual  impact  of 
the  physical  soil  properties  in  their  association  with  one  another 
serves  as  the  basis  for  the  following  key.  The  profile  is  regarded 
as  a  physical  system  and  those  physical  properties  capable  of 
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observation  and  record  'which  may  incorporate  other  features  or  prop¬ 
erties,  either  physical,  chemical,  or  biological  are  used  as-  the 
critical  ones  to  distinguish  subtle  changes  in  the  environment  and 
contaminant  behavior. 

An  attempt  has  been  made  to  identify  and  code  most  of  the 
relevant  parameters  which  best  describe  the  physical  characteristics 
of  the  soil.  •  The  following  list  of  variables  are  to  be  measured  at 
each  sampling  site  in  the  comprehensive  survey  at  RMA.  This  list 
includes  a  coding  key  which  has  been  successfully  used  in  previous 
studies  and  should  be  used  to  enter  the  data  gathered  into  the  com¬ 
puterized  data  bank.  The  scale  for  each  of  the  variables  describing 
topography,  site  drainage,  physical  soil  properties,  and  bedrock 
lithology  has  been  designed  so  that  the  scale  can  serve  as  input 
directly  into  numerical  analysis  and  computer  graphics  routines. 

The  key  is  flexible  in  that  if  necessary,  it  can  be  expanded  without 
detriment  to  those  portions  of  the  key  already  completed. 

PHYSICAL  VARIABLES 

Each  sampling  site  will  be  identified  by  a  ten  digit  index 
located  on  the  State  Plane  Coordinate  System.  The  first  five 
characters  of  the  index  represent  the  X  axis  of  the  east-west  State 
Plane  Coordinates  measured  to  the  nearest  foot,  and  the  last  five 
characters  represent  the  Y  axis. 

A.  Topographic  Variables 

1.  Elevation  is  the  height  of  a  given  point  above  mean 
sea  level.  Elevation  is  a  critical  variable  which  is  used  in  most 
calculations  to  provide  information  on  relative  position. 

Units  =  feet  (height  above  mean  sea  level) 

2.  Local  relief  is  the  difference  in  elevation  between  the 
highest  and  lowest  points.  It  measures  the  physical  size  of  the 
topographic  features. 

Units  =  feet  (difference  between  highest  and  lowest 

points) 

3.  Ruggedness  is  a  combination  of  gully  depth  measured 
in  feet  and  gully  count  per  100  square  yard  site.  Ruggedness  mea¬ 
sures  the  degree  of  direction  of  a  landform  and  provides  information 
about  the  degree  and  stage  of  landform  erosion.  It  also  indirectly 
measures  the  pattern  of  surface  water  flow. 
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(a)  Depth  in  feet 


>  1  foot -  1 

i  1-6  feet  - -  2 

6-15  feet  — - .  3 

I  >15  feet  -  4 

’(b)!  Count  across  100  yards  square  site 

None -  0 

1  -  5 - - — "  1 

5-10 -  2 

>  10 - 3 


4.  Slope  describes  the  angular  inclination  of  the  surface 
at  the  sampl ing  site.  It  is  measured  as  a  ratio  of  the  amount  of 
drop  (h)  in  100  feet  of  length  (£.)• 

Units  =  ratio  h  -  c  (recorded  as  a  decimal  number) 

5.  Slope  geometry  represents  the  shape  of  the  sloping  sur¬ 
face  and  is  an  important  variable  in  characterizing  soil  erosion  or 
soil  accumulation  on  the  ground  surface. 

(a)  Shape 


Convex  -  1 

Concave  -  2 

Straight  -  3 

(b)  Relative  location 

Top -  1 

Middle  -  2 

Bottom  -  3 


6.  Aspect  describes  the  directional  orientation  of  a  topo¬ 
graphic  feature  with  reference  to  north.  it  is  most  convenient  to 
divide  the  compass  direction  into  16  divisions. 


N.N.E.  -  01 

N.E.  -  02 

E.N.E.  -  03 

E.  -  04 
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E.S.E.  -  05 

S.E.  - 06 

5.5. E.  - 07 

$ #  -  08 

5.5. W . . . - . —  09 

S.W.  — - . 10 

W.S.W.  . 11 

W.  - - - - -  12 

W.N.W.  — . 13 

N.W.  - 14 

N.N.W.  -  15 

N.  - 15 


7.  Drainage  is  characterized  by  two  different  components. 
Site  drainage  describes  the  movement  or  potential  for  movement  of 
water  and  moisture  —  the  amount  of  water  contained  in  the  soil. 

(a)  Site  drainage 


Shedding -  01 

Normal  -  02 

Receiving -  03 

Regular  flooding,  Short 

Term -  04 

Regular  flooding,  Long 

Term - 05 

Occasional  flooding,  Short 

Term -  06 

Occasional  flooding,  Long 
Term -  07 


(b)  Moisture 


Dry  (becomes 
D,  (air  dry) 


D2 


darker  when  we 


t) 


°3 


01 

02 

03 


Moderately  moi st 

MM,  -  04 

MM2 -  05 

MM3  (wilting  point)  -  06 


Moi  st 


M1  - - 

M2 - 

M3  (field  capacity) 


07 

08 

09 
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Wet 


10 

(material  will 

wet 

11 

and  stick  to 

fingers 

12 

when  molded) 

B.  Soil  Physical  Properties 

1.  Depth  to  sample  in  the  vertical  distance  measured  from 
the  ground  surface  to  the  top  and  the  bottom  of  the  sample. 

units  =  feet  (measured  to  the  nearest  hundreth  of  foot) 


2.  Horizon  is  a  layer  within  the  soil  profile  haying  mor¬ 
phological  characteristics  and  properties  different  from  those  layers 
which  occur  below  and/or  above  it.  Horizons  are  usually  parallel  to 
the  land  surface,  except  where  tongues  of  material  from  one  horizon 
may  penetrate  vertically  into  the  neighboring  horizon.  The  horizons 
are  recognized  by  the  nature  of  their  organization.  A0  results  largely 
from  biologic  processes.  The  A  horizon  results  from  biologic  processes 
and  those  physico-chemical  processes  that  occur  at  surfaces.  The  B 
horizon  results  from  physio-chemical  processes  which  occur  below  sur 
faces  that  are  relatively  deep  seated.  C  horizons  are  recogn i zed  by^ 
their  lack  of  pedologic  development  and  remains  of  geologic  organiza¬ 
tion.  The  D  horizons  result  from  some  earlier  cycle  of  soil-forming 
process;  and  thus,  their  pedologic  organization  contrasts  with  that 
found  in  the  solum  of  the  soil  with  which  they  are  associated. 

Weathered  bedrock  is  an  example. 

Eight  categories  will  be  recognized  in  this  study.  They 
are  given  in  descending  order:  the  AQ  horizon  is  thesurface  soil 
layer  and  is  composed  of  organic  material;  the  A^  horizon  is  the  top¬ 
most  mineral  layer  with  a  high  content  of  humified  organic  matter, 
relatively  dark  in  color  and  has  a  maximum  biological  activity  for 
any  given  soil  profile;  the  A2  horizon  is  the  zone  of  maximum  leach¬ 
ing  (eluviation),  it  may  be  defined  by  color  difference.  It  may  be 
bleached  to  a  white  color  or  sporadic  by  bleaches  which  appear 
irregular  or  as  blocks  less  than  one-fourth  inch  in  thickness  at 
the  interface  with  the  B  layer.  The  B  horizon  is  a  transition  zone 
between  the  A  and  B  horizons  though  accumulation  of  material  (illu¬ 
viation)  starts  to  take  place.  The  B  horizon  is  characterized  by  a 
concentration  of  clay  and/or  iron  and/or  aluminum  and/or  translocated 
organic  material  and/or  having  a  structure  or  consistency  different 
from  A  or  C.  B  also  has  stronger  colors.  The  B2  horizon  is  the  zone 
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of  maximum  illuviation  of  mineral  and  organic  matter;  the  B-C  horizon 
i  ]  s  a  transition  zone  between  the  illuvial  B  horizons. and  the  uncon 
isolidated  material  of  the  C  horizon.  The  C  horizon  is  unconsolidated 
material  which  may  or  may  not  be  of  the  same  parent  material  rom  w  ic 
the  u^per  horizonl  are  forced.  It  is  not  in  the  zone  of  major  b,olog,cai 
activities;  the  D  horizon  is  weathered  bedrock. 

Horizon 


3.  Horizon  thickness  —  represents  the  size  of  the  horizon. 

units  =  feet  measured  to  nearest  .01  of  a  foot. 

k.  Horizon  boundary  is  a  zone  of  change  occurring  between 
soil  horizons  ranging  from  the  sharp  and  clearly  defined  to  diffuse. 

The  boundary  defines  the  nature  of  the  change  and  is  spec  1 f . ed  by  two 
terms  -  one  a  measure  of  the  thickness  or  width  of  the  transition 
zone  between  horizons,  the  other  a  measure  of  its  shape. 

(a)  Shape 

1  -  Even  =  boundary  is  almost  a  plane  surface  - _ 

2  -  Wavy  =  boundary  waves,  depressions  formed  are 

wider  than  their  depths  - -  .  . 

3  -  irregular  =  boundary  waves,  depressions  formed 

are  deeper  than  their  widths 

h  -  Broken  =  boundary  is  discontinuous  1 

(b)  Sharpness 

1  -  Sharp  =  <  2  cm 

2  -  Clear  =  2  -  5  cm 

3  -  Gradual  =  5  ”  10  cm 
k  -  Diffuse  =  >  10  cm 

5  Color  is  defined  in  terms  of  hue,  value,  and  chroma  us'n9 
the  Munsell  s^TTTolor  charts.  Hue  is  a  function  of  the  wave  length  of 
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the  light.  Value  is  the  intensity  of  color  and  chroma  the  relative 
purity  or  saturation  of  a  color.  Color  should  be  recorded  for  moist 

soil  only. 

Color  (moist) 


1  Hue 


2,  5  YR . — .  32 

5  YR - -  16 

7,  5  YR -  8 

10  YR  - . - .  * 

2,  5  Y -  2 


2  Value 

1-2-3-A-5-6-7-8 


3  Chroma 

l_2-3-4-5-6-7-8 

6.  pH  measures  the  hydrogen  ion  concentration  of  soil.  It 
is  a  fundamental  property  of  both  physical  and  chemical  charactei  is- 
tics  of  both  the  organic  and  inorganic  constituents  of  the  soil. 

units  =  dimensionless  number  of  nearest  .1 

7.  Soil  texture  is  a  measure  of  the  size  and  proportion  of 
individual  gain  present  in  the  soil.  Texture  can  be  measured  by  the 
behavior  of  a  small  handful  of  soil  when  moistened  and  kneaded  into. a 
ball  and  then  pressed  out  between  thumb  and  forefinger.  The  approxi¬ 
mate  percentage  content  of  clay  (particles  less  than  0.002  mm  in  dia¬ 
meter)  and  silt  (particles  between  0.02  and  0.002  mm  in  diameter)  and 
sand  (range  from  2  mm  to  0.02  mm  diameter  with  an  arbitrary  separation 
of  coarse  and  fine  sand  at  0.2  mm  particle  diameter)  provide  a  guide 
for  classification. 


a.  The  field  method  of  determining  texture 
take  a  sample  of  soil  sufficient  to  comfortably  fit  into 
the  hand.  The  soil  is  moistened  with  water,  a  little  at 


soil  i s  to 
the  palm  of 
a  time,  and 
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kneaded  until  the  ball  of  soil,  so  formed,  just  fails  to  stock  to  the 
fingers.  More  soil  or  water  may  be  added  to  attain  this  condition, 
which  is  known  as  the  sticky  point  and  approximates  field  moisture^ 
capacity  for  that  soil.  Kneading  and  moistening  are  continued  until 
there  is  no  apparent  change  in  the  soil  ball,  usually  a  working  time 
of  one  to  tw6  minutes.  The  soil  ball,  or  bolus,  is  now  ready  for 
shearing  manipulation;  bu  t  the  behavior  of  the  soil  during  bolus 
formation  is  also  indicative  of  its  texture.  The  behavior  of  the 
bolus  and  of  the  ribbon  produced  by  shearing  (pressing  out)  between 
thumb  and  forefinger  characterizes  the  texture. 

b.  A  much  more  reliable  method  of  determining  soil 
texture  is  by  sieving  the  soil  samples. 

Textural  classes 


S  Sand  -  01 

LS  Loamy  sand -  02 

CS  Clayey  sand -  03 

SL  Sandy  loam -  Oh 

FSL  Fine  sandy  loam -  05 

GL  Gravelly  loam -  06 

L  Loam -  07 

SCL  Sandy  clay  loam -  08 

GCL  Gravelly  clay  loam  --  09 

CL  Clay  loam -  10 

SI  CL  Si  lty  clay  loam -  11 

CG  Gravel  1  y  clay -  12 

FSCL  Fine  sandy  clay  loam  -  13 

SC  Sandy  clay -  l*f 

SIC  Si  lty  clay -  15 

LC  Light  clay -  16 

MC  Medium  clay -  17 

HC  Heavy  clay - -  18 


8.  Stones  are  measurable  surface  features.  The  abundance 
is  measured  as  a  percentage  of  cover;  size  consists  of  eight  categories 
ranging  from  gravel  (0.2  -  1  cm)  to  large  boulders  (<  600  cm);  shape 
can  be  angular,  subangular,  rounded,  shaley,  or  tabular;  and  the 
lithology  of  stones  is  given  in  12  petrological  types. 
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(a)  Abundance 

Stoneless 
Slightly  stony 
Stony 

Very  stony 
Extremely  stony 
j  Stone  dominant 

(b)  Size 

Gravel 
Smal  1 
Med i urn 
Large 

Very  large 
Small  boulders 
Medium  boulders 
Large  boulders 


{%  cover) 


<  1  -  0 

1  -  5  - .  1 

5  -  20 -  2 

20  -  50  -  3 

50  -  75  -  * 

>  75 -  5 

(cm) 

0.2-1  - . —  1 

1  -  5 -  2 

5  -  1  -  3 

10-20 -  4 

20  -  60  .  5 

60  -  200  -  6 

200  -  600  -  7 

>  600  . .  8 


(c)  Shape 

Angular 
Subangul ar 
Rounded 

Sha 1 ey  - 

Tabular 


(d)  Lithology 

Flint  - 

Soapstone  - 

Limestone  - 

Mudstone  - 

Siltstone  and  shale 

Sandstone  - 

Ironstone  - 

Quartzite  - 

Granite  - 

Gneiss  or  schist  — 

Ryhol ite  - 

Andesite  - 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 
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9.  Roots  describe  the  presence  of  vegetative  matter  in  the 
soil  measured  in  terms  of  quantity  and  size.  Decay  of  roots  provides 
channels  of  water  flow  within  the  soil. 

(a)  Quantity 


Absent  -  0 

Present  (trace)  -  1 

Very  low -  2 

Low -  3 

Moderate  -  b 

High -  5 

(b)  Size 

Sma  11  -  1 

Sma  1 1  -  Med  ium -  2 

Smal  1  -  Large -  3 

Med i um  -  Sma 11  -  b 

Med ium  -  5 

Med i um  -  Large  -  6 

Large  -  Smal  1  -  7 

Large  -  Medium -  8 

Large -  9 


10.  Structure  is  concerned  with  the  arrangement  of  all  soil 
particles.  Particles  refer  to  all  peds  and  nonpeds  (see  both  Apedal 
and  Peds).  Structure  may  be  described  in  terms  of  the  three  charac¬ 
teristics  of  the  intensity,  geometry,  and  size  of  the  soil  aggregates 
as  foil ows : 


(a)  Intensity  expresses  the  degree  and  strength  of  soil 
aggregation  and  may  vary  with  the  moisture  status  of  the  soil.  For 
purposes  of  standardization,  intensity  should  be  determined  when  the 
soil  is  inthe  moderately  moist  to  moist,  moisture  stages.  Four  inten¬ 
sities  of  structure  are  usually  defined  as  follows: 

(1)  Structureless  in  which  there  is  no  observable 

aggregation  and  the  soil  is  massive,  if  coherent,  and  single  grain,  if 
not  coherent.  (Structureless  is  the  same  as  Apedal.)  -  0 

(2)  Weak  in  which  some  peds  are  discernible  and 

when  disturbed,  less  than  one-third  of  the  soil  material  is  found  to 
consist  of  peds.  -  1 
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(3)  Moderate  in  which  peds  are  clearly  seen  and 

iwhen  disturbed,  one-third  to  two-thirds  of  the  soil  material  is  found 
1  to  consist  of  peds.  - 2 

(4)  Strong  in  which  peds  are  clearly  seen  and 

when  disturbed,  more  than  two-thirds  of  the  soil  material  is  found 
to  consist  of  peds.  - 3 

NOTE  that  both  moderate  and  strong  grades  of 
structure  qualify  a  soil  for  Peds  evident,  whereas  the  structure¬ 
less  state  and  weak  grade  of  structure  qualify  a  soil  for  Peds,  few, 
if  any. 

(b)  Geometry 


Angular  blocky  or  cubic  (ABL)  -  1 

Subangular  blocky  or  cubic  (SABL)  -  2 

Angular  prismatic  (APR)  -  3 

Subangular  prismatic  (SAPR)  -  4 

Polyhedral  >  six  sides  (PH)  -  5 

Paral lelepedal  (PLL)  -  6 

Platy  (PL)  - - —  7 

Spheroidal -sol  id  (SSPH)  - 8 

Spheroidal -vi s i bl e  voids  (VSPH)  -  9 


(c)  Size 

Given  in  inches  or  fraction 

11.  Plasticity  is  a  characteristic  of  wet  soils.  It  is  the 
ability  of  soil  to  change  its  shape  as  a  result  of  stress  and  to  keep 
that  shape  after  the  removal  of  the  stress. 


Nonplastic  - 0 

SI i ght ly  plastic  -  1 

Plast  ic -  2 

Very  plastic -  3 


12.  Friability  is  a  characteristic  of  moist  soils.  It  is  a 
measure  of  resistance  to  crushing. 


Loose -  0 

Very  friable -  1 

Friable  -  2 

Fi  rm -  3 

Very  firm -  4 

Extremely  firm -  5 


13.  Porosity  -  a  porous  soil  material  is  one  characterized 
by  the  presence  of  voids  in  the  soil  moss.  It  is  the  volume  of  the. 
total  soil  bulk  not  occupied  by  solid  particles.  The  first  component 
of  porosity  is  void  quantity,  which  is  the  percentage  of  nonsolid  space 
ranging  from  nonporous  (<  1%)  to  extremely  porous  (>  50%).  The  second 
component  is. pore  size,  measured  in  millimeters.  Porosity  is  largely 
a  function  of  the  arrangement  of  soil  particles  into  aggregates,  the 
type  of  soil  particles,  amount  of  organic  material,  and  degree  of 
compaction. 

(a)  Void  quantity 


Nonporous  <  ~  1  0 

SI  i ghtly  porous  1  -  5 -  1 

Porous  5-10 -  2 

Porous  10-20  -  3 

Very  porous  20  -  50  -  ^ 

Extremely  porous  >  50 -  5 

(b)  Pore  size 

Very  fine  <  ~  0.5 -  1 

Fine  0.5  “  1  2 

Med  ium  1“3 -  3 

Large  3  "  5 - 

Very  large  >  5 -  5 


]b.  Mottling  is  the  presence  of  spots  or  alternate  streaks 
of  color  that  differ  from  the  overall  color  of  the  soil.  Mottling 
occurs  in  zones  of  alternating  periods  of  good  and  poor  aeration. 

The  differences  in  oxidation  states  of  iron  and  magnesium  result  in 
mottled  appearances.  Mottles  are  considered  to  be  masses,  blobs,  or 
blotches  of  subdominant  colors,  such  that  color  places  them  in  a 
different  value/chroma  rating  to  the  matrix  color.  The  glaze  or 
color  skin  which  may  occur  on  the  outside  of  some  peds  ss  not  con¬ 
sidered  as  mottling.  A  soil  is  only  considered  to  be  mottled  when 
10%  or  more  of  mottles  occur  in  the  soil  mass.  For  mottled  soils, 
the  value/chroma  rating  of  the  dominant  matrix  color  is  taken.  Mot¬ 
tling  has  three  components:  abundance,  measured  as  a  percentage; 
size  in  mm;  and  sharpness  of  the  boundary  between  mottled  and  non- 
mottled  soil,  measured  in  mm. 

(a)  Abundance  (%  cover) 

None  <  2 -  0 

Few  2-5  1 
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Common 

5-20 - 

3 

Definite  matrix 

color 

20  -  50  . 

4 

No  matrix  color 

>  20 - 

5 

Definite  matrix 

color 

>  50 . . 

6 

'•(b) 

Size 

(mm) 

Extremely  fine 

<  1  - - 

1 

Very  fine 

1  -  2 - 

2 

Fine 

2-5 . — 

3 

Hed i urn 

5  -  15  — . 

4 

Large 

>  15 - 

5 

(c) 

Sharpness 

(mm) 

Diffuse 

>  2  — . . 

1 

Cl  ear 

c  2 - - 

2 

Sharp 

(knifedged)  - 

3 

15.  Secondary  precipitates  result  from  the  upward  capillary 
movement  and  removal  of  soil  water" through  evaporation  at  the  ground 
surface.  This  process  is  important  in  contaminant  movement  and  con¬ 
centration  in  the  soil. 

(a)  Quantity 


None 

Few 

Common 

Common 

Abundant 

Extra  abundant 


0 


<  2% . .  1 

2  -  5% -  2 

5  -  20% -  3 

20  -  50% .  4 

>  50% -  5 


(b)  Composition 

Lime - 

Magneste  - 

Gypsum  - 

Silica  - 

Fe  oxides  - 

Fe  oxides  +  humus 
Mn  oxides  - 
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(c) 


Type 


Irregular  concretions  ‘  [ 

Regular  concretions  -  2 

Layers  -  8 

16.  *  Bedrock 

(a)  Rock  structure 

(1)  Intensity 

W.  Weak -  ' 

M.  Moderate - -  2 

S.  Strong -  3 

(2)  Geometry 

B.  Banded -  01 

S.  Scoriaceous -  02 

A.  Amygdaloidal  -  03 

P.  Porphyritic -  04 

L.  Laminated -  05 

M.  Massive -  °6 

C.  Cleared -  07 

S.  Sheared - ; -  08 

F.  Contorted -  09 

X.  Cross-bedded -  10 

(3)  Size 

V.  Very  coarse -  1 

C.  Coarse -  2 

F.  Fine -  3 

VF  Very  fine -  ^ 

(b)  Rock  type 

Consolidated  sediments  (0  -  19) 

CG  Conglomerate -  01 

qj  —Quartzite -  02 

SS  Sandstone -  03 

ST  Si  ltstone -  04 

MS  Mudstone -  05 

CS  Claystone -  06 

SH  Shale -  07 

CH  Chert - -  08 
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GW  Greywacke  - 

AK  Arkose  - 

LS  Limestone  - 

IS  Ironstone - - - 

Tl  Til  1  ite . . — 

Igneous  rocks  (20  -  39) 

GR  Granite  - 

GD  Granodiorite  - 

Dl  Diorite  - 

RH  Rhyol i te  - 

BA  Basalt  - 

AG  Agglomerate  - 

TU  Tuff  - 

SC  Scorria  - 

DO  Dolerite  - 

DB  Disabase  - 

PG  Pegmatite  - 

AP  Apl i te  - 

RQ.  Reef  quartz - 

PY  Porphyry  - 

PT  Porphyrite  - 

Metamorphic  rocks  (40  -  59) 

SL  Slate - . - 

PH  Phyl 1 ite  - 

HF  Hornfels  - 

SH  Schist  - 

GN  Gneiss  - 


09 

10 

11 

12 

13 


20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 


40 

41 

42 

43 

44 
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INSTALLATION  RESTORATION 


INSTALLATION  RESTORATION 


INSTALLATION  RESTORATION 


SAMPLING  PLAN  FOR  RMA 


Dynamic  Cluster  Sampling 

An  attractive  attribute  of  dynamic  sampling  is  that  as  informa¬ 
tion  is  gathered,  it  can  be  used  to  direct  further  sampling  without 
biasing  the  probability  properties  of  the  sample. 

Data  on  soils  and  groundwater  will  be  gathered  from  drill  sites 
patterned  in  clusters  of  16  bore  holes  per  cluster,  arranged  in  a 
1)  x  If  square  grid.  Each  sample  will  be  surveyed,  using  the  x  and  y 
coordinates  of  the  State  Plane  Grid  System  which  in  addition,  will 
serve  as  the  identification  code  for  each  of  the  samples.  The  ele¬ 
vation  will  also  be  recorded.  Field  data  on  land  surface  topography, 
physical  properties  of  soils,  groundwater  levels, and  bedrock  will  be 
logged  directly  onto  the  coding  forms.  Soil  samples  will  be  collected 
at  the  appropriate  depth;  sent  to  the  soils  laboratory  where  it  wi 
be  divided  into  two  parts.  One  part  of  the  sample  will  go  to  the  MALD 
(Material  Analysis  Lab  Division)  for  chemical  analysis  of  various  con¬ 
taminants.  The  other  portion  will  be  analyzed  in  the  soi 1  laboratory 
for  soil  moisture,  mechanical,  and  mineral  properties..  While  the 
analyses  are  being  performed,  a  second  cluster  of  16  sites  will  be 
drilled  and  sampled.  The  results  from  the  first  cluster  of  samples 
will  be  mapped  to  determine  the  positioning  of  the  third  cluster  of 
16  sampling  sites.  The  cluster  will  be  located  in  the  direction  of 
the  maximum  concentration  of  the  contaminant  being  mapped.  The  analy¬ 
zed  results  from  the  second  cluster  will  likewise  be  used  to  determine 
the  positioning  of  the  fourth  cluster,  etc.  This  technique  of  mapping 
each  contaminant  and  positioning  the  cluster  of  samples  in  the  d i rec 
tion  of  maximum  concentration  will  permit  tracking  the  pollution  plume 
to  its  origin  with  minimum  exploratory  drilling.  Periodically,  the 
cluster  will  be  combined,  and  analysis  will  be  performed  on  the  total 
combined  samples.  This  will  yield  an  increase  in  the  sample  size  and 
an  increased  relative  homogeneity  of  the  elements  being  sampled.  By. 
combining  the  cluster  into  a  total  sample  one,  in  effect,  is  evaluating 
the  regional  pattern  of  the  pollutants  as  compared  to  more  local,  larger 
scale  behavior  within  each  cluster.  However,  sampling  theory  requires 
that  all  the  elements  of  the  cluster  be  sampled  for  the  results  to  be 
statistical ly  val id. 

It  is  possible  that  several  plumes  will  be  encountered,  each 
originating  from  a  different  source  and  containing  their  own  unique 
grouping  of  compounds.  When  this  occurs,  each  plume  will  be  traced 
to  locate  all  sources. 


On-going  exploratory  drilling  along  the  northern  boundary  of 
RMA  has  revealed  the  existence  of  at  least  two  contaminated  groun 
water  plumes.  One  plume  contains  DIMP,  CP,  and  Na  compounds.  The 
other  is  a  DCPD  and  pesticide  plume.  Because  of  the  high  priority 
placed  upon  containment  of  contaminant  migration,  this  part  of  the 
Arsenal  would  be  a  logical  place  to  start  the  sampling  program.  An 
alternate  site  would  be  at  the  known  contaminant  hot  spots. 
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Samp! ing  Grid  Spacing 

A  test  sampling  grid,  spaced  250  feet  apart,  will  be  constructed 
initially  to  pretest  the  sampling  method  and  data-handl ing  procedure. 
The  distance  between  the  sampling  bore  holes  will  be  set  as  more  data 
becomes  available  on  the  diffusion  rates  of  contaminant  concentration 
with  distance.  The  mathematical  relationship  between  change  in  con¬ 
centration  level  per  unit  change  in  distance  will  be  calculated  for 
each  contaminant,  and  the  grid  spacing  will  be  adjusted  to  the  opti¬ 
mum  distance.  Preliminary  tests  show  that  the  best  fit  regression 
equation,  equating  the  change  in  DIMP  concentration  with  distance 
measured  from  the  center  of  the  DIMP  plume  is  Y  =  A  +  B  log  X.  The 
DIMP  concentration  values  in  the  inclosed  sample  contained  DIMP 
added  from  secondary  sources  carried  by  groundwater.  This  caused 
additive  harmonics,  which  should  be  filtered  out  before  a  usable 
equation  can  be  obtained. 

Sampl i no  Hor i zons 

One-foot  sections  of  the  sampling  core  will  be  collected  to 
serve  as  the  samples  for  both  mechanical  and  chemical  analysis. 

Soil  samples  will  be  taken  at  six  depths,  so  as  to  increase  the 
likelihood  of  intercepting  the  zones  of  contaminant  concentration. 
These  zones  are:  (1)  The  A  soil  horizon,  which  is  the  zone  of  con¬ 
taminant  built  through  (a)  bioaccumulation  in  the  soil  organic 
matter;  (b)  upward  migration  of  contaminants  by  capillary  action 
and  surface  evaporation  of  water;  and  (c)  accumulation  of  pollutants 
on  the  ground  surface  by  both  wind  or  surface  water  transport,  (  ) 

The  B  soil  horizon  which  is  the  zone  of  illuviation  or  accumulation^ 
of  compounds  leached  from  the  upper  horizons;  (3)  A  point  mid-way  in 
the  zone  of  aeration  representing  the  transition  zone;  (A)  Thecapil 
lary  fringe  immediately  above  the  ground  water  table;  (5)  A  point 
below  the  top  of  the  water  table;  and  (6)  At  the  bedrock  surface.  . 
‘The  soil  samples  taken  in  the  zone  of  saturation  will  provide  a  ratio 
of  the  contamination  in  the  soil  compared  with  that  moving  in  the 


i  groundwater.  Additional  samples  will  be  taken  at  the  top  of  clay 
Slayers  or  lenses  whenever  they  are  encountered,  to  determine  the 
1  effect  of  clay  as  a  point  of  contaminant  concentration.  Clay  layers 
in  the  regolith  may  behave  as  selective,  physical,  chemical,  or  osmo 
-tic  filter  aided  by  the  greatly  reduced  rate  of  downward  percolation 
of  groundwater  and  the  possible  presence  of  clays  having  a  high  ionic 
exchange  capacity.  Groundwater  will  also  be  sampled  at  similar  depths. 

The  typical  profile  shown  in  the  diagram  was  constructed  by 
averaging  the  horizon  of  each  of  the  soil  series  occurring  at  RMA. 

The  A  horizon  ranges  from  0-9  inches  in  thickness,  while  the  B  hori¬ 
zon  ranges  from  9  -  26  inches. 

In  several  places  on  the  Arsenal ,  the  B  horizon  does  not  occur; 
and  the  soils  grade  directly  from  the  A  to  the  C  horizons.  This 
usually  occurs  in  younger,  poorly  developed  azonal  soils.  In  order 
to  avoid  missing  data  for  sites  where  horizons  are  absent,  all  sampl¬ 
ing  will  consistently  be  taken  at  uniform  depths  predetermined  to. 
coincide  with  the  soil  profile  horizons.  These  depths  are  given. in 
the  schematic  illustrating  the  sampling  scheme  for  a  typical  boring. 

Two  water  samples  will  be  collected  at  each  drill  one,  one  approxi 
mately  three  feet  below  the  water  table  leve,  the  other  at  bedrock 
surface. 

Numerical  Analysis  and  Mapping 

Each  of  the  contaminants  will  be  mapped  individual!  for  each  of 
the  six  horizons  sampled,  thus  providing  a  three-dimensional  view  of 
the  contaminant  distribution.  The  maximum  concentration  encountered 
at  each  site  will  also  be  mapped  to  show  a  plain  view  of  the  highest 
levels  occurring  on  the  Arsenal. 

Principal  component  factor  analysis  will  be  peformed  to  determine 
the  correlations  and  behavior  of  the  contaminants  as  groups  or  associa 
tions.  The  information  will  be  obtained  from  the  factor  loadings,  which 
will  identify  new  groupings  indicating  that  new  pollution  plumes  were, 
encountered.  The  maps  of  the  factor  scores  for  each  of  the  factors  will 
provide, insight  on  the  spatial  distribution  of  each  of  the  contaminant 
groups  individually  and  will  serve  to  identify  the  location  and  direc 
tion  of  the  plume,  thus  indicating  where  the  next  sampling  cluster. is 
to  be  positioned.  Factor  analysis  of  the  total  set  of  variables  will, 
through  the  correlation  coefficient,  provide  the  linkage  web  between 
the  contaminants  and  the  physical  ecosystem. 


The  relationships  and  the  linkages  will  be  mappable  through  the 
factor  scores.  Multivariate  discriminate  analysis,  using  the  factor 
scores  as  surrogate  weights,  will  serve  as  a  comparison  of  similarity 
and  differences  between  different  regions  representing  either,  dif¬ 
ferent  contaminated  areas  of  different  source  areas.  The  multivariate 
disciminant  function  can  also  serve  as  a  model  for  simulation. 
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PREFACE 


The  production,  testing,  and  disposal  of  various  chemical  and 
biological  substances  at  RMA  has  lead  to  suspected  contamination  of 
land  areas,  water,  ecology,  buildings,  and  equipment.  Little  infor¬ 
mation  exists  as  to  the  concentration  levels,  overall  extent,  or  the 
diffusion  process  of  these  substances  within  the  ecosystem.  The  need 
to  identify  and  assess  the  potential  hazards  created  by  the  presence 
of  toxic  substances  in  the  environment  at  RMA  has  been  recognized, 
and  efforts  are  underway  for  interim  containment  of  off-post  migration. 
A  comprehensive  survey  will  be  initiated  to  gather  and  analyze  all 
necessary  information  which  would  enable  a  rational  assessment  of  the 
degree  of  contamination  on  the  post  and  provide  the  data  base  needed 
for  the  design  of  remedial  measures. 

The  purpose  of  this  study  is  to  provide  the  guidelines  for  soil 
and  water  sampling  procedures  and  field  data  gathering  and  recording 
protocols  to  be  used  in  the  comprehensive  survey  at  RMA.  The 
approach  developed  is  aimed  at  obtaining  the  maximum  likelihood  of 
intercepting  the  contaminated  horizons  in  three  dimensions  and  thus 
gaining  the  optimum  benefit  for  long-run  objectives  from  the  field 
measurements  with  minimum  field  effort  and  thus  reducing  time  and 
cost  by  minimizing  the  possibility  of  gathering  extraneous  data. 


INTRODUCTION 


A  pollutant  is  any  compound  or  substance  that  is  undesirable 
where  it  is.  A  contaminant  or  pollutant  may  become  undesirable  by 
being  introduced  into  an  environment  that  is  foreign  to  it.  The 
pollutant  is  usually  toxic  to  the  host  environment  and  will  disturb 
the  equilibrium  state  of  the  ecosystem.  The  degree  of  concern  of  the 
presence  of  contamination  depends  upon  the  toxicity  level  to  living 
organisms,  especially  man. 

Pollutants  may  be  either  man-made,  released  as  byproducts  of 
industrial  or  military  processes,  or  are  naturally  occurring  and 
become  toxic  when  concentrated,  for  example  salt. 

Transport  processes  often  create  a  pollutant  simply  by  moving  a 
compound  to  a  different  environment.  Sources  of  pollutants  are 
usually  classed  as  "point"  and  "nonpoint. "  Control  measures  are 
simplified  whenever  the  pollution  plume  is  traced  to  a  point  source, 
as  compared  with  diffuse  or  nonpoint  pollutants  which  discharge  into 
waterways  by  dispersal  pathways. 


FACTORS  INFLUENCING  CONTAMINANT  BEHAVIOR 


The  success  in  locating  the  most  probable  areas  of  contamination 
and  in  tracing  the  breadth  (source  and  plume)  of  such  contamination 
within  the  ecosystem  is  dependent  upon  an  accurate  interpretation  of 
the  means  of  contaminant  emplacement.  The  dispersion  and  the  accumu¬ 
lation  of  pollutants  in  the  environment  is  dependent  on  a  number  of 
inter-related  factors.  Some  of  the  principal  factors  are: 

The  Physical  and  Chemical  Properties  of  the  Contaminants  -  The 
important  chemical  properties  of  the  pollutants  are  the  following: 
the  contaminant  phase,  concentration,  chemical  stability,  water 
solubility,  ion  exchange,  adsorption,  equilibrium  with  the  host 
environment,  and  regional  geochemical  properties.  These  properties 
determine  the  chemical  behavior  of  the  contaminants,  whether  they 
are  present  as  a  solution  in  water  or  are  dispersing  as  a  solid. 

The  Form  of  Transport  -  The  form  of  transport  determines  the 
dispersion  process  and  influences  accumulation.  It  is  possible  by 
evaluating  the  manner  that  the  pollutants  were  transported  and 
knowing  their  pathways  to  predict  the  major  contaminated  areas. 
Contaminants  of  greatest  concern  at  RMA  are  transported  either: 

(1)  as  solid  particles  (particulate)  occurring  in  a  crystal ized 
form  or  combined  with  dust  or  soil  particles,  or  (2)  as  a  fluid 
in  solution  with  water.  For  example,  the  movement  of  pesticides 
on  sediment  is  affected  by  the  enrichment  process.  Pesticides  are 
adsorbed  primarily  on  organic  soil  colloids  which  remain  in  suspen¬ 
sion  longer  than  the  coarse  soil  particles.  Most  mobile  herbicides, 
on  the  other  hand,  rarely  contaminate  the  groundwater  because  the. 
adsorption  capacity  in  most  soils  above  the  water  table  is  sufficient 
to  retain  the  chemicals  until  they  decompose.  However,  in  an. area 
with  a  high  percolation  potential,  groundwater  contamination  is 
poss i bl e. 

The  partitioning  of  the  compounds  between  water  and  soil  trans¬ 
port  is  necessary  in  determining  the  mass  movement  of  chemicals  in 
run-off.  The  quantity  of  a  compound  dissolved  in  water  is  obtained 
from  water  solubilities  of  the  chemicals.  The  amount  adsorbed  at 
different  solution  concentrations  can  be  calculated  by  mixing  solu¬ 
tions  of  the  chemicals  with  known  quantities  of  soil.  (This  study 
is  currently  in  progress  for  DiMP  by  the  Aerojet  Ordnance  and 
Manufacturing  Company,  Downey,  California.) 
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The  transporting  medium  for  particulate  contaminants  is  either 
jwind  or  surface  water  movement.  The  contaminants,  therefore,  con¬ 
centrate  in  areas  controlled  by  these  transporting  media  —  on  or 
close  to  the  ground  surface.  The  transporting  media  for  liquid 
contaminant,  on  the  other  hand,  is  either  surface  run-off,  subsurface 
run-off,  or  groundwater  run-off.  These  contaminants  tend  to  concen¬ 
trate  much  deeper  in  the  unconsolidated  mantle. 

The  Characteristics  of  the  Physical  Site  -  The  characteristics 
'of  the  site  include:  CD  Geologic  properties  of  the  aquif ier  such 
as  texture,  porosity,  permeability,  transmissivity,  heterogeneity  of 
the  subsurface  regolith,  thickness  of  the  unsaturated  zone  of  aera¬ 
tion,  thickness  of  the  saturated  groundwater  zone  boundary  conditions, 
and  (2)  The  geomorphology,  the  topography,  and  the  surficial  cover. 

This  includes  slope,  aspect,  surface  relief  (high  relief  is  i nd i cat i ve 
of  high  sediment  delivery  ratio),  channel  density,  the  condition  of  the 
channels  (whether  clogged  or  open,  meandering,  or  straight),  stream 
gradients,  the  relief/length  ratio.  The  geomorphology  of  a  terrain 
provides  potential  channels,  sinks,  or  traps  for  pollutant,  and  (3) 

The  stratigraphy  and  character  of  the  sediments  in  the  aquifer.  This 
would  include  the  ionic  exchange  capacity  of  the  minerals  in  clay  and 
soil,  which  has  a  direct  bearing  on  contaminant  sorption  and  leaching. 
The  type  of  both  primary  and  secondary  minerals  present  and  the  oxida¬ 
tion  reduction  reactions  which  affect  some  of  the  pollutants  and  the 
leach  ability  of  the  sediments. 

Cl imate  -  Humid  environments  are  those  in  which  the  annual  pre¬ 
cipitation  exceeds  evapotranspi rat  ion  rates.  They  have  a  net  water 
surplus  which  transports  the  mobile  compounds  faster  and  to  greater 
depths.  In  arid  climates  such  as  RMA,  much  of  the  water  is  removed  by 
evapotranspi rat  ion ;  and  soil  moisture  is  brought  to  the  surface  by 
capillary  action.  Contaminated  groundwater  may  be  transported  to  the 
surface  and  the  contaminants  are  deposited  on  the  ground  surface.  The 
duration  and  intensity  of  storm  activity  influences  groundwater  recharge 
and  possible  flushing  of  the  contaminant  from  the  soil.  Flow  is  upward 
during  periods  of  high  evapotransp i rat  ion ,  provided  there  is  no  ponding 
and  downward  when  water  is  available  (as  under  the  water  ponds  and  bog 
areas).  The  micro  climate  also  influences  the  viscosity  and  surface 
tension  properties  of  soil  water  solution,  which  is  temperature- 
dependent.  Soils  drain  slowly  during  cold  weather,  and  water  is 
retained  as  storage  when  frozen.  As  a  solution  freezes,  salts  con¬ 
centrate  in  the  nonfrozen  liquid,  and  it  is  not  unusual  to  find  brines 
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in  partially  frozen  soils.  These  brines  become  selectively  mobile 
because  of  their  hydraulic  properties.  During  the  spring  after  the 
frozen  ground  thaws  at  RMA,  much  of  the  groundwater  recharge  results 
from  snowmelt.  If  the  concentrations  of  pollutants  in  soil  water 
reaches  high  enough  levels,  precipitation  (or  salting  out)  may  occur. 

The  Geohydrolic  Conditions  -  Hydrologic  transport  is  the  move¬ 
ment  of  compounds  in  or  by  water  and  may  occur  along  surface  or  sub¬ 
surface  channels  at  RMA.  Hydrologic  transport  is  therefore  dependent 
'upon  the  amount  of  water  and  the  pathways  of:  (1)  Surface  run-off 
traveling  over  the  ground  surface;  (2)  Subsurface  run-off  that  has 
infiltrated  the  surface  soil  and  moves  laterly  through  the  upper  soil 
horizon  within  the  zone  of  aeration;  and  (3)  Groundwater  run-off 
below  the  groundwater  table  in  the  zone  of  saturation.  Each  must  be 
known  before  the  transport  of  pollutants  and  hydrodynamic  dispersion 
can  be  described  in  detail. 

The  hydraulic  conductivity  (K  variable  in  Darcy's  Law)  of  sedi¬ 
ments  is  the  transmissive  properties  of  a  porous  medium.  Generally, 
in  hydrolic  calculations,  this  value  is  averaged  vertically  assuming 
a  homogeneous  aquifer  and  substituting  it  for  transmissivity,  expres¬ 
sed  as  T  =  Kh,  where  h  is  the  saturated  thickness.  The  transmissivity 
can  be  determined  by  pump  tests. 

An  analysis  of  the  dynamics  of  groundwater  flow  permits  a  reliable 
computation  of  the  level  of  contamination  of  the  aquifer  in  space 
through  time.  Changes  in  chemical  concentrations  in  the  groundwater 
system  are  due  to:  (1)  Convective  transport,  expressing  the  amount 
of  dissolved  chemicals  moving  with  the  flowing  groundwater;  (2)  Addi¬ 
tion  or  removal  of  contaminants  through  chemical  and  physical  reactions 
taking  place  in  the  water  or  between  water  and  soil;  (3)  Mixing  with 
recharge  water  of  different  concentrations;  and  (4)  Hydrodynamic 
dispersion  in  which  molecules  or  ionic  diffusion  and  minute  changes 
in  velocity  of  flow  in  the  porous  media  causes  the  paths  to  diverge 
from  the  normal  direction  of  groundwater  movement. 

Grouping  -  Certain  elements  tend  to  group  together  into  compound 
associations  in  a  specific  environment  due  to  similar  migration  capa¬ 
cities  (relative  mobilities).  The  tendency  for  compounds  to  group 
together  simplifies  the  mapping  task  in  that  by  testing  for  either 
indicator  or  pathfinder  elements,  it  is  possible  to  trace  the  spatial 
behavior  of  the  pollutants. 
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Contaminant  Migration  through  Soil  Erosion  -  Several  contami¬ 
nated  areas  exist  on  RMA  which  are  sources  of  ground  surface  pollu¬ 
tion  by  soil  movement.  Field  soil  losses  have  been  estimated  to 
range  from  three  to  eight  tons  per  acre  per  year.  A  ton  of  dry  soil 
is  approximately  one  cubic  yard  in  volume.  Spread  uniformly  over  an 
acre  of  surface,  it  would  have  a  depth  of  less  than  one-hundredth  of 
an  inch.  Erosion  hazards  both  by  wind  and  due  to  rainfall  and  surface 
run-off  are  local  in  nature.  This  can  be  significant  in  areas  which 
are  not  protected  by  a  vegetative  cover  (example  Basins  A,  C,  and 
"most  of  the  Arsenal  during  the  winter  months  when  vegetation  is  dor¬ 
mant).  The  annual  soil  losses  from  sheet  and  rill  erosion  by  rain¬ 
fall  can  be  determined  by  the  Universal  Soil  Loss  Equation  (USLE) . 

The  equation  is: 

A  -  J  ( RKLSCP) 

Where  A  is  estimated  average  annual  soil  loss  measured  in  tons 
per  acre. 

R  is  the  rainfall  and  run-off  erosivity  index. 

K  is  the  soil  erodibility  factor.  it  is  the  average  soil  loss 
per  unit  of  R  and  depends  on  soil  properties  (see  Figure  ). 

LS  is  dimensionless  topographic  factor  that  represents  the  com¬ 
bined  effects  of  slope,  length,  and  steepness. 

C  is  the  cover  and  management  factor. 

P  is  the  support  -  practice  factor  for  land  use. 

'The  USLE  can  provide  a  measure  of  the  contaminated  sediment 
delivered  to  a  stream. 

The  object  of  the  RMA  environmental  sampling  program  is  to  be 
able  to  identify,  measure,  and  characterize  each  of  the  contaminated 
horizons. 

General ly,  immediately  underneath  a  contaminated  source,  each 
layer  in  the  unconsolidated  mantle  collects  the  contaminants  from 
above,  mixes,  absorbs,  reacts,  and  then  drains  the  soil  to  field 
capacity  to  supply  lechate  to  the  layer  below.  Thus,  the  process 
proceeds 
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FIGURE  i 


STRATIFICATION  OF  THE  UNCONSOLIDATED  MANTLE 


When  a  vertical  section  is  made  through  the  unconsolidated  mantle, 
a  series  of  distinct  layers  or  zones  are  observed  which  differ  in  their 
chemical  and  physical  properties.  See  Figure  The  uppermost  A  hori¬ 
zon  varies  from  a  few  inches  to  several  feet  in  thickness  and  repre¬ 
sents  the  eluvial  (washed  out)  layer.  Varying  amounts  of  organic 
material  is  added  to  this.  Materials  removed  from  the  A  horizon  are 
carried  downward  in  solution  and  are  deposited  in  the  B  horizon; 
however,  the  fine  colloidal  clays  are  mechanically  removed  by  descend¬ 
ing  soil  water  (illuvial  layer).  Impervious  hardpan  and  claypan 
layers  may  develop  in  the  B  horizon  by  the  accumulation  of  downward 
moving  material  which  include  iron  and  aluminum  minerals,  including 
clays.  The  A  and  B  horizon  together  form  the  solum.  In  arid  cli¬ 
mates,  calcium  carbonate,  sodium  sulphates,  and  other  salts  may 
accumulate  just  below  the  B  horizon  in  the  C  layer.  The  C  layer  or 
parent  material  does  not  exhibit  a  high  degree  of  alteration;  however, 
it  may  display  weathering,  including  oxidation,  hydration,  and  lime 
accumul at i ons . 

A  mature  soil  profile  exhibits  well  developed  horizons;  however, 
layering  may  be  indistinct  or  lacking  in  young  soils.  Topography 
affects  both  soil  profile  development  and  contaminante  emplacement. 

On  steep  slopes,  the  soil  horizons  are  much  thinner  or  may  never 
develop  because  erosion  removes  weather  products  as  rapidly  as  they 
form.  Topography  controls  the  rate  of  surface  run-off,  the  infiltra¬ 
tion  of  water,  the  position  of  the  water  table,  the  depth  of  penetra¬ 
tion  of  chemical  weathering  processes,  and  contaminant  migration. 

Local  and  regional  relief  affect  the  hydrolic  gradient  which  governs 
the  circulation  of  water  through  the  soil  and  unconsolidated  mantle. 

Knowledge  of  the  soil  profiles  developed  within  a  suspected  con¬ 
taminated  area  is  necessary  to  enable  a  constant  relative  sampling 
procedure  to  be  maintained.  All  samples  must  be  taken  from  the  same 
horizon,  otherwise  they  will  not  be  comparable.  Also,  a  uniform 
sampling  procedure  increases  the  probability  of  locating  the  con¬ 
taminated  zones.  See  Figure 


An  analysis  of  the  vegetative  cover  can  provide  valuable  infor¬ 
mation  about  the  character  of  the  contamination  below  ground  surface. 
For  example,  deep-rooted  plants  may  penetrate  the  pollution  plume  in 
the  overburden  and  be  brought  to  the  surface  and  enriched  in  the  soil 
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THE  SOIL  PROFILE 


Organic,  original  forms  recognized 
Organic,  original  forms  not  recognized 

Mineral,  mixed  with  humus,  dark  colored 

Horizon  of  maximum  eiuviation  of 
silicate  clays,  Fe,  At  oxides  etc. 

Transition  to  B,  more  like  A  than  B 
Transition  to  A,  more  like  B  than  A 

Maximum  illuviation  of  silicate  clays, 

Fe,  Al  oxides,  some  organic  matter 

Transition  to  C,  more  like  B  than  C 

Zone  of  least  weathering,  accumulation 
of  Ca,  Mg  carbonates,  cementation, 
sometimes  high  bulk  density 


Nutrient  recycling  is  an  important  factor  in  determining  the  relation¬ 
ship  between  vegetation  and  the  soil  that  develops. 


Figure  3 


Bioaccumulation  of  DIMP  by  Plant  Parts  in 
20  ppm  Irrigation  (37  days  from  Original  Inoculation) 
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Root 

45.  6 

2.  28 

Stem 

37.  1 

1. 86 

Leaf 
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6.  46 

? 

Carrot  - 

9200 

184 

Root 

12.4 

0.  62 

Stem 

6.6 

0.33 

Leaf 

36.  9 

• 

1.85 

Bean  - 

9200 

184 

Root 

45.4 

2.27 

Stem 

28.  9 

1.45 

Leaf 

150.0  ' 

7.  50 

Wheat  - 

9200 

184 

Root 

31.  5 

1. 58- 

Stem 

14.  2 

0.  71  ’ 

Leaf 

105.  5 

5.  28 

Figure  4 
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Bioaccumulation  of  DIMP  by  plant  parts. 

Soil  culture,  37  days  exposure  to  20  ppm  DIMP  in  irrigation  water 
Figure  5 


,  layer  by  decay  of  the  vegetation  in  the  production  of  humus.  The 
| concent  rat  ion  of  elements  in  plants  is  dependent  on  a  variety  of 
•physical,  chemical,  and  organic  factors  which  control  intake  of 
nutritional  solutions.  Geobotanic  survey  may  be  valuable  to  monitor 
contaminant  movement  in  the  soil,  especially  in  areas  where  soils 
profiles  are  young  or  poorly  developed.  See  Figure  Results 
from  biochemical  sampling  tend  to  be  hetrogeneous ,  and  care  should 
be  taken  to  assure  representative  samples  and  their  interpretation. 
Seasonal  variations  in  concentration  occur  within  plant  tissue  in 
response  to  nutritional  requirements  of  the  plants.  Different  vege¬ 
tative  species  have  selective  absorption  capacities  for  intake  for 
different  elements  which  may  vary  considerably.  Lysimeter  studies 
conducted  by  Aerojet  Ordnance  and  Manufacturing  Company,  Downey, 
California,  have  indicated  that  DIMP  will  bioaccumulate  at  much 
higher  concentrations  in  the  leaves  of  plants  than  either  plant 
roots  or  stems  (see  Figure  ).  The  leaves  account  for  the  greatest 
proportion  of  organic  matter  addition  to  the  soils  which  results  in 
a  concentration  of  pollutants  in  humus  in  the  A1  horizon.  The  infor¬ 
mation  gathered  from  geobotanical  surveys  will  be  used  in  conjunction 
with  soil  and  water  surveys. 


SOILS  OF  THE  ROCKY  MOUNTAIN  ARSENAL 

\ 

There  are  eight  soil  series  and  four  land  types  at  RMA.  A 
soil  series  is  a  subdivision  of  soil  families  --  the  most  generalized 
soil  description.  A  soil  series  is  a  "group  of  soils  developed  from 
the  same  kind  of  parental  material  by  the  same  genetic  combination  of 
processes  and  whose  horizons  are  quite  similar  in  their  arrangement 
and  general  characteristics."  Soils  series  are  divided  into  specific 
mapping  units. 

ASCALON  SERIES  -  Consists  of  well  drained  soils  occurring  on 
level  to  moderately  sloping  land  on  uplands.  They  are  formed  in 
loamy  material  containing  varying  amounts  of  sand  and  gravel.  The 
infiltration  capacity  of  Ascalon  soils  is  moderate  to  high,  and  the 
available  soil  moisture  is  high.  This  soil  is  suitable  for  plant 
root  development. 

Representative  Profile  -  The  surface  layer  is  brown  sandy  loam 
approximately  six  inches  thick  and  noncalcareous.  The  upper  part  of 
the  subsoil  is  brown  sandy  loam  and  sandy  clay  loam  about  15  inches 
thick.  It  is  noncalcareous.  The  lower  part  of  the  subsoil  is  brown 
sandy  loam  about  six  inches  thick,  highly  calcareous,  and  contains 
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much  visible  lime.  The  underlying  material  is  a  highly  calcareous, 
pale-brown  to  very  pale  brown  fine  sandy  loam.  It  is  found  at  a 
depth  of  27  to  60  inches  or  more. 

Mapping  Units 

(AsB)  Ascalon  sandy  loam,  one  to  three  percent  slopes.  This 
unit  is  found  on  upland  ridges  and  benches.  Surface  run-off  is 
slow  to  medium,  and  the  hazard  of  water  erosion  is  slight  to  model 
ate.  Soil  blowing  is  a  severe  hazard  in  unprotected  areas. 

(AsC)  Ascalon  sandy  loam,  three  to  five  percent  slopes.  This 
unit  is  found  in  areas  adjacent  to  major  streams,  but  it  is  only  on 
the  higher  upland  ridges  and  side  slopes.  It  has  a  thinner  surface 
layer  and  subsoil  than  the  representative  profile.  Surface  run-off 
is  medium.  The  hazard  of  water  erosion  is  moderate.  Soil  blowing 
is  severe  in  unprotected  fields. 

(AsD)  Ascalon  sandy  loam,  five  to  nine  percent  slopes.  Located 
on  uplands  near  the  edges  of  the  breaks  adjacent  to  stream  channels. 
The  surface  layer  and  subsoil  are  thinner  than  the  representative 
prof i le. 

(ArC)  Ascalon-vona  sandy  loams,  one  to  five  percent  plagues. 
Ascalon  sandy  loam,  one  to  three  percent  slopes  (AsB)  makes  up  45 
percent  of  this  complex  and  occupies  the  smoother  and  lower  lying 
areas.  Vona  sandy  loam,  three  to  five  percent  slopes  make  up  35 
percent  and  is  normally  in  the  more  sloping  part  of  the  landscape. 

A  few  small  flowouts  are  on  the  higher  ridge  points. 

BLAKELAND  SERIES  -  Consists  of  somewhat  excessively  drained 
soils  occurring  on  gently  sloping  to  moderately  sloping  on  uplands. 
These  soils  formed  in  wind  and  water  worked  sandy  material.  Blake- 
land  soils  absorb  water  rapidly,  and  the  available  water  capacity 
is  low.  The  entire  profile  becomes  unusually  hard  as  it  dries. 

When  dry,  these  soils  hold  a  nearly  vertical  bank.  Permeability 
is  rapid.  The  entire  soil  is  suitable  for  plant  roots. 

Representative  Profile  -  The  surface  layer  is  grayish-brown 
loamy  sand  about  five  inches  thick  and  noncal careous.  The  next 
layer  is  brown  sandy  loam  about  four  inches  thick  and  is  noncal- 
careous.  The  underlying  material  from  a  depth  of  about  3  inches 
to  60  inches  or  more  is  nonca 1 careous ,  brown  loamy  sand  and  pale- 
brown  sand. 
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(BoD)  Blakeland  loamy  sand  three  to  nine  percent  slopes.  Occurs 
on  uplands.  Surface  run-off  is  medium  to  slow.  The  hazard  of  soil 
blowing  is  severe. 

(Bt)  Blakeland  -  truck  ton  association  -  Blakeland  loamy  sand, 
three  to  nine  percent  slopes,  constitutes  at  least  60  percent  of 
this  association,  occurring  in  areas  of  semidune  relief,  though  it 
is  also  on  long,:  narrow  ridges.  Truck  ton  loamy  sand,  three  to 
nine  percent  slopes  makes  up  20  -  40  percent  of  mapped  areas.  It 
is  depressional  or  has  the  more  gentle  slopes.  The  hazard  of  soil 
blowing  is  moderate  but  is  very  severe  where  adequate  cover  is  not 
maintained. 

(Gr)  Gravelly  Land  -  Shale  Outcrop  Complex  -  this  complex  con¬ 
sists  of  steep  and  hilly  land  types.  Areas  are  generally  elongated 
and  roughly  parallel  river  channels.  The  complex  consists  of  very 
shallow,  clayey  soils  and  raw  shale  overlain  by  a  discontinuous 
layer  of  gravel  deposits  four  inches  to  many  feet  thick.  These 
deposits  are  old  enough  that  some  soil  is  forming  in  places  within 
the  gravelly  cap  layer.  Little  if  any  soil  is  forming  in  the  shaly 
exposed  areas,  which  is  eroding  by  moderate  wind  and  water  action. 
Some  areas  contain  steep  colluvium  of  mixed  materials.  All  of  this 
complex  is  subject  to  severe  erosion. 

Loamy  Alluvial-  Only  one  unit  of  loamy  alluvial  is  present  at 

RMA. 


(Lu)  Loamy  alluvial  land  is  found  in  dra i nageways .  Slopes  are 
less  than  three  percent.  It  is  subject  (in  varying  degrees)  to 
flooding  from  adjacent  slopes  and  main  stream  channels. 

The  profile  consists  of  a  surface  layer,  commonly  dark-colored 
loam  or  clay  loam  about  six  to  ten  inches  thick  and  normally  non- 
calcareous.  The  underlying  material  is  stratified  loam,  silt  loam, 
and  clay  loam  containing  varying  amounts  of  fine  sand,  sand,  and 
fine  gravel.  It  is  20  to  rfiore  than  60  inches  thick,  brown  to  dark 
brown  in  color,  and  in  places  it  has  a  weak  structure.  It  is  nor¬ 
mally  calcareous.  Sand  or  sand  and  gravel,  stratified  with  their 
lenses  of  silt  and  loam,  are  in  some  places  below  a  depth  of  36 
i nches. 
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.  Water  is  absorbed  at  a  moderate  rate,  and  their  available  water 

icapacity  is  high.  They  are  normally  well  drained.  Natural  fertility 
'is  high.  The  hazard  of  water  erosion  is  severe  in  all  cultivated 
areas. 

NUNN  SERIES  -  Consists  of  wel  1  drained,  nearly  level  soils  on 
terraces.  These  soils  formed  in  loamy  alluvial  material.  Water  is 
absorbed  at  a  moderate  rate,  and  the  available  water  capacity  is 
high.  Permeability  is  slow  in  soils  with  a  clay  loam  surface  layer 
and  moderately  slow  in  those  with  a  loam  surface  layer.  The  entire 
soil  is  suitable  for  plant  roots. 

Representative  Profile  -  The  surface  layer  is  a  grayish-brown 
loam  about  six  inches  thick  and  noncal careous.  The  subsoil  is  dark 
grayish-brown  and  grayish-brown  clay  about  17  inches  thick.  It  is 
noncal careous.  The  underlying  material  below  a  depth  of  23  inches 
is  light  yellowish-brown  and  grayish-brown  loam  and  silt  loam.  It 
is  highly  calcareous  and  stratified. 

(NIA)  Nunn  loam,  zero  to  one  percent  slopes.  This  soil  is  on 
stream  terraces  along  most  of  the  major  drainageways  in  the  county. 
Surface  run-off  is  slow,  and  the  hazard  of  water  erosion  is  slight. 

The  hazard  of  soil  blowing  is  moderate  in  cultivated  fields  without 
protective  cover. 

The  profile  of  this  unit  is  similar  to  the  reprensentat ive  pro¬ 
file,  except  for  a  thicker  surface  layer  and  subsoil  and  less  sus¬ 
ceptible  to  erosion. 

(NuB)  Nunn  clay  loam,  one  to  three  percent  slopes.  This  soil  is 
found  on  stream  terraces  along  major  streams.  Areas  are  long  and  nar¬ 
row  and  roughly  parallel  to  stream  channels.  Run-off  is  medium,  and 
the  hazard  of  water  erosion  is  moderate.  Soil  blowing  can  be  control¬ 
led  by  adequate  cover.  The  profile  is  similar  to  the  representative 
profile,  but  it  has  a  more  clayey  surface  layer  and  a  thinner  subsoil. 

PLATNER  SERIES  -  Consists  of  well  drained  nearly  level  to  gently 
sloping  soils  on  uplands.  These  soils  are  formed  in  old  alluvium. 

They  absorb  water  slowly,  and  the  available  water  capacity  is  high. 
Permeability  is  slow,  and  the  entire  soil  is  suitable  for  plant  roots. 

Representative  Profile  -  The  surface  layer  is  grayish-brown 
loam  about  nine  inches  thick  and  noncalcareous.  The  upper  part  of 
the  subsoil  is  brown  clay,  about  nine  inches  thick.  It  is  noncal¬ 
careous.  The  lower  part  of  the  subsoil  is  light-gray  clay  loam 
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about  10  inches  thick.  It  is  highly  calcareous,  and  part  of  the  lime 
is  visible  as  splotches.  The  underlying  material  below  a  depth  of 
28  inches  is  very  pale  brown  loam  that  is  highly  calcareous.  At  a 
depth  of  about  kS  inches,  it  is  a  white  and  very  pale  brown  sandy 
loam  that  is  highly  calcareous  and  contains  some  fine  gravel. 

( P 1 B )  Planter  loam,  zero  to  three  percent  slopes.  This  unit  is 
found  on  uplands.  Surface  run-off  is  slow,  and  the  hazard  of  water 
erosion  is  moderate  to  slight.  The  hazard  of  soil  blowing  is  severe 
in  dry  farmed  areas  unless  protected  by  plants  or  stubble. 

(Smj  Sandy  Alluvial  Land  "  Consists  of  an  unstable  accumulation 
of  gravelly  and  sandy  alluvium.  It  is  in  and  adjacent  to  beds  of 
intermittent  streams.  During  periods  of  heavy  rain,  the  streambeds 
are  subject  to  flooding  and  channels  are  relocated  and  sediment  is 
shifted  and  redeposited  at  slightly  different  locations.  During  dry 
spells,  this  land  type  is  extremely  droughty. 

The  material  of  this  soil  was  transported  by  water  from  the 
sand  and  gravel  beds  in  or  adjacent  to  the  area.  It  is  stratified 
because  of  periodic  flooding.  Thin  lenses  or  small  pockets  of  silt, 
clay,  and  sand  are  also  mixed  with  the  gravel. 

This  land  type  differs  from  wet  alluvial  land  primarily  in  that 
it  is  coarser  textured  throughout  and  is  not  affected  by  the  water 
table. 

STONEHAM  SERIES  -  Consists  of  well  drained,  nearly  level  to 
moderately  sloping  soils  on  uplands.  These  soils  formed  in  loamy, 
old  alluvial  material.  Water  is  absorbed  at  a  moderate  rate,  and 
the  available  water  capacity  is  moderate.  Permeability  is  moderate.- 
The  entire  soil  is  suitable  for  roots. 

Representative  Profile  -  The  surface  layer  is  light  grayish- 
brown  loam  about  five  inches  thick  and  is  noncal careous.  The  sub¬ 
soil  is  brown  sandy  clay  loam  about  eight  inches  thich  and  noncal- 
careous.  The  underlying  material  is  light-gray  sandy  loam  and 
gravelly  loam  that  is  highly  calcareous.  At  a  depth  of  30  inches, 
it  is  pale-brown  gravelly  sandy  loam  that  is  highly  calcareous. 

It  extends  to  a  depth  of  about  60  inches. 


(StD)  Stoneham  loam,  three  to  nine  percent  slopes,  is  found  on 
uplands  near  major  stream  channels  in  scattered  areas.  Run-off  is 
generally  medium  (but  rapid  if  the  soil  is  unprotected).  The  hazards 
of  water  erosion  and  soil  blowing  are  severe  in  cultivated  areas. 

TRUCK-TON  SERIES  -  Consists  of  well  drained  to  somewhat  exces¬ 
sively  drained,  nearly  level  to  moderately  sloping  soils  on  uplands. 
These  soils  formed  in  wind-worked  sandy  material.  They  absorb  water 
rapidly,  and  the  available  water  capacity  is  low.  Permeability  is 
rapid  and  the  entire  soil  is  suitable  for  roots. 

Representative  Profile  -  The  surface  layer  is  grayish-brown 
and  dark  grayish-brown  loamy  sand  about  nine  inches  thick.  It  is 
nonca  1  careous .  The  upper  part  of  the  subsoil  is  dark-brown  and 
brown  sandy  loam  about  12  inches  thick.  It  is  nonca 1 careous .  The 
lower  part  of  the  subsoil  is  yellowish-brown  loamy  sand  about  11 
inches  thick;  also  noncalcareous.  The  underlying  material  is  yellow¬ 
ish-brown  coarse  sand  that  extends  to  a  depth  of  60  inches  or  more. 

(TeD)  Truck-ton  loamy  sand,  three  to  nine  percent  slopes. 

Found  mostly  along  major  stream  channels.  Surface  run-off  is  medium, 
and  the  hazards  of  soil  blowing  and  water  erosion  are  severe. 

(TuB)  Truck-ton  sandy  loam,  one  to  three  percent  slopes.  The 
profile  is  similar  to  the  representative  profile  except  for  a  sandy 
surface  layer  and  a  thicker  surface  layer  and  subsoil  surface  run¬ 
off  is  slow  to  medium,  and  the  hazard  of  water  erosion  is  slight. 

The  hazard  of  soil  blowing  is  severe. 

(TuD)  Truck-ton  sandy  loam,  two  to  nine  percent  slopes.  The 
profile  is  similar  to  the  representative  profile  except  for  a  sandy 
loam  surface  layer.  Run-off  is  medium,  and  the  hazard  of  water  ero¬ 
sion  is  moderate.  The  hazard  of  soil  blowing  is  severe. 

VONA  SERIES  -  Consists  of  well  drained,  nearly  level  to  moder¬ 
ately  sloping  soils  on  uplands.  They  were  formed  in  wind-deposited 
material.  The  soils  absorb  water  rapidly,  and  the  available  water 
capacity  is  moderate.  Permeability  is  rapid,  and  the  entire  soil  is 
suitable  for  plant  roots. 

Representative  Profile  -  The  surface  layer  is  light  brownish- 
gray  loamy  sand  about  nine  inches  thick.  It  is  noncalcareous.  The 
subsoil  is  pale-brown  coarse  sandy  loam  about  13  inches  thick;  also 
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is  noncalcareous.  The  underlying  material  is  light  yellowish-brown 
sandy  loam  and  very  pale  brown,  highly  calcareous  loamy  sand  that 
contains  lime  visible  as  splotches. 

(VnD)  Vona  loamy  sand,  three  to  nine  percent  slopes.  Surface 

run-off  is  medium,  and  the  hazard  of  soil  blowing  is  severe. 

(VoC)  Vona  sandy  loam,  three  to  five  percent  slopes.  The  pro¬ 

file  is  similar  to  the  representative  profile  except  for  a  sandy 
loam  surface  layer  and  a  thinner  subsoil.  The  hazards  of  water 
erosion  and  soil  blowing  are  severe. 

WELD  SERIES  -  Consists  of  well  drained,  nearly  level  soil  on 
uplands"!  The  so i  1  s  were  formed  in  wind-worked  loamy  materials.  The 
soils  absorb  water  at  a  moderate  rate,  and  the  available  water  capacity 
is  high.  Permeability  is  slow,  and  the  entire  soil  is  suitable  for 
plant  roots. 

Representative  Profile  -  The  surface  layer  is  brown  loam  about 
six  inches  thick.  It  is  noncalcareous.  The  upper  part  of  the  sub¬ 
soil  is  dark  brown  clay  about  six  inches  thick.  It  is  noncalcareous. 
The  lower  part  of  the  subsoil  is  pale-brown  and  very  fine  sandy  loam 
about  20  inches  thick.  It  is  highly  calcareous,  and  much  of  the  lime 
is  visible  as  splotches.  The  underlying  material  is  highly  calcar¬ 
eous,  pale-brown  silt  loam  and  light  yellowish-brown  fine  sandy  loam. 

(WmB)  Weld  loam,  one  to  three  percent  slopes.  This  soil  is 
found  on  uplands.  Surface  run-off  is  medium,  and  the  hazard  of 
water  erosion  is  moderate  to  severe.  The  hazard  of  soil  blowing  is 
severe  in  clay  farmed  areas,  if  rainfall  is  below  normal. 

(Wt)  Wet  Alluvial  Land  ~  This  is  found  on  the  nearly  level  bot¬ 
tom  lands  of  the  larger  streams  next  to  stream  channels.  The  land 
is  wet  most  of  the  time  and  is  flooded  by  streamflow  once  to  several 
times  a  year  during  periods  of  high  water. 

Materials  are  extremely  variable  in  texture,  consisting  of 
stratified  layers  of  dark-colored  silt,  loam,  and  clay.  The  layers 
are  generally  less  than  six  inches  thick  and  are  underlain  by  sand, 
fine  sand,  and  some  gravel  at  depths  of  one  to  three  feet.  They  are 
wet  at  a  depth  of  two  feet  most  of  the  time  and  are  commonly  wet  to 
the  surface  throughout  the  growing  season.  Water  tolerant  plants 
form  the  vegetative  cover. 
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SAMPLING  METHODS 


One  of  the  objectives  of  the  comprehensive  survey  is  to  appraise 
the  internal  organization  of  contaminant  distribution;  the  location  of 
the  elements  of  the  distribution  with  respect  to  each  other.  The 
"spatial  structure"  represents  the  location  of  each  element  relative 
to  each  of  the  others  and  also  the  location  of  each  element  relative 
to  all  of  the  contaminants  taken  together;  that  is,  the  internal  loca¬ 
tional  organization  of  contaminant  distribution  in  space.  The  "spa¬ 
tial  processes,"  on  the  other  hand,  are  mechanisms  which  produce  the 
spatial  structure  of  the  contaminant  distribution. 

Spatial  structure  is  a  determinant  of  spatial  process  as  much 
as  process  is  a  determinant  of  structure.  For  example,  the  geo¬ 
hydrologics  characteristics  of  the  Arsenal  environment  control 
groundwater  movement,  which  also  affects  contaminant  migration  and 
the  spatial  pattern  of  contamination.  An  understanding  of  the  spa¬ 
tial  process  can  be  achieved  by  the  deductive  approach,  using  the 
system's  concept  which  requires  fewer  samples  then  would  be  required 
to  map  out  the  spatial  pattern. 

Contaminant  emplacement  and  migration  will  be  approached  from 
the  systems  viewpoint.  Systems  modeling  necessitates  a  clear  defi¬ 
nition  of  the  objectives  of  the  study,  to  identify  and  model  the 
components  of  the  system  andsubsystems ,  and  to  elucidate  their  mutual 
relationships,  especially  with  regard  to  the  spatial  patterns. 

System  modeling  has  the  ability  to  handle  a  dynamic  process. 

It  has  the  ability  to  look  at  the  total  ecosystem,  makes  efficient 
use  of  information  gathered,  has  the  ability  of  checking  consistency, 
has  the  ability  of  information  integration,  is  applicable  both  as  a 
planning  tool  and  an  analysis  tool,  and  is  an  excellent  means  of 
organizing  data  for  decision  makers.  The  process  response  system, 
incorporating  both  morphological  and  cascading  subsystems,  can  serve 
as  a  working  model  to  demonstrate  the  linkage  between  form  and  process. 
The  strength  and  direction  of  connectivity  between  the  components  is 
analyzed  by  statistical  methods.  Such  methods  can  range  from  analy¬ 
sis  of  variance,  correlation  and  multiple  regression,  to  cross- 
spectral  analysis. 

The  problem  in  sampling  and  the  selection  of  a  sampling  design 
for  meaningful  data  collection  is  predicated  by  the  research 
obj  ecti ves. 
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The  purpose  of  sampling  is  to  accurately  portray  the  nature  of 
the  environment  from  which  the  samples  were  taken.  These  are  two 
standard  categories  of  sampling  methods,  probability  sampling  and 
nonprobability,  or  judgment  sampling. 

1.  Probab? 1 i ty  Sampl ing  -  Probability  sampling,  based  upon  random 
selection,  i's  used  in  modeling  and  forecasting  the  behavior  of 
naturally  occurring  environmental  systems.  !t  requires  that  each 
element  have  an  equal  chance  of  being  chosen;  that  is,  independent 
of  any  event  in  the  selection  process.  Random  selection  procedures 
eliminate  conscious  or  unconscious  bias  by  the  scientist  and  enables 
the  researcher  to  use  probability  theory  as  a  method  for  estimates 
of  population  parameters  and  estimates  of  error. 

a.  Simple  Random  Sampling  -  After  a  sampling  frame  has  been  esta¬ 
blished,  numbers  are  assigned  to  each  of  the  elements  in  the  list. 
The  elements  to  be  measured  are  then  selected  randomly  from  the  list, 
If  spatial  patterns  are  to  be  sampled,  then  a  grid  overlay  is  used 
and  the  X  and  Y  coordinates  are  selected  randomly  to  identify  either 
the  coordinate  intersection  to  be  measured  as  a  point  sample  or  the 
appropriated  grid  cell  to  be  measured  as  an  area  samples. 

b.  Systematic  Sampling  -  The  first  element  is  selected  randomly, 
then  every  Kth  element  in  the  total  list  is  chosen  for  inclusion  in 
the  sample.  Periodicity  and  cyclical  patterns  that  coincide  with 
the  sampling  interval  may  be  missed  and  may  bias  the  results. 

c.  Strati fied' Sampl ing  -  Stratification  organizes  the  population 
into  homogeneous  subsets  (with  heterogeneity  between  subsets)  and 
an  appropriate  number  of  elements  are  measured  from  each.  Stra¬ 
tification  ensures  proper  representation  of  the  variables  describing 
the  parent  population.  An  appropriate  number  of  elements  are  drawn 
from  homogeneous  subsets  of  the  population.  This  minimizes  sampling 
error. 

d.  Multi-Stage  Cluster  Sampling  -  Cluster  sampling  involves  sampl¬ 
ing  groups  of  elements  (clusters)  followed  by  selection  of  elements 
within  the  clusters.  Multi-stage  cluster  sampling  involves  the 
repetition  of  listing  and  sampling.  As  the  homogeneity  increases 
within  each  cluster,  relatively  fewer  elements  may  be  needed  to 
adequately  represent  a  given  natural  cluster.  However,  a  larger 
number  of  clusters  may  be  needed  as  the  diversity  among  clusters 
increases. 
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e.  Hierarchied  Sampling  (Nested  Sampling)  -  A  region  is  first  sub¬ 
divided  into  major  areas  of  equal  size.  Several  major  areas  are  then 
chosen  at  random  and  partitioned  into  smaller  sub  regions.  The  sub 
regions  are  randomly  selected  and  further  subdivided.  Hierarchied 
sampling  ensures  that  every  part  of  the  target  region  is  represented. 

f.  Disproportionate  Sampling  and  Weighing  -  A  probability  sample 
is  one  in  which  each  population  element  has  a  known  nonzero  prob¬ 
ability  of  selection,  even  though  different  elements  may  have  dif¬ 
ferent  probabilities.  If  controlled,  probability  sampling  procedures 
can  be  used  and  may  be  representative  of  the  population  from  which 

it  was  drawn,  if  each  sample  element  is  assigned  a  weight  equal  to 
the  inverse  of  its  probability  of  selection. 

2.  Nonprobability  Sampl ing;a^  Purpos i ve  or  judgment  sampling  --  a 
researcher  may  sample  with  prior  knowledge  of  the  population,  its 
elements  and  research  objectives.  This  is  useful  to  pretest  a 
sample  plan  or  to  select  a  typical  or  representative  sample. 

b. «  Quota  Sampling  -  Quota  sampling  begins  with  a  matrix  describ¬ 
ing  the  characteristics  of  the  target  population.  The  ration  of 
the  population  occurring  in  various  categories  is  known.  A  pre¬ 
determined  proportion  is  selected  for  sampling  from  each  cell  of 
the  matrix. 

c.  Traverses  -  Traverse  sampling  is  the  technique  used  to  sample 
the  environment  along  a  straight  line  whose  length  and  direction 
have  been  determined.  Traverse  sampling  may  be  used  to  gather  data 
without  disturbing  the  environment.  The  assumption  that  the  sum  of 
the  intercepts  along  a  traverse  line  can  be  used  to  estimate  the 
proportion  of  the  total  area  occupied  by  the  phenomenon  being  inter¬ 
cepted.  The  intercept  lengths  reflect  the  spatial  geometry  and  the 
degree  of  fragmentation  of  the  occurrences  of  the  attribute  being 
sampled  in  the  region. 

Probability  sampling  avoids  conscious  or  unconscious  bias  in 
element  selection  which  closely  represents  the  population  of  all 
elements.  Probability  sampling  permits  estimates  of  sampling  error. 
The  selection  of  a  particular  sampling  design  depends  upon  the  accur¬ 
acy  and  precision  required;  the  efficiency  and  feasibility  of  a  par¬ 
ticular  procedure;  and  the  cost  in  time,  manpower,  and  money.  The 
design  selected  also  depends  upon  the  objectives  of  the  investigation, 
the  structure  or  functional  characteristics  of  the  population. 
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A  combination  of  the  purposive  sampling  approach  (which  allows 
intensive  study  of  structure  and  relationships  and  therefore  hypo¬ 
thesis  formulation)  with  a  probability  sampling  approach  (which 
allows  the  hypotheses  to  be  formally  tested  with  respect  to  a  given 
population)  is  advocated  for  environmental  studies. 

Estimating  the  Sample  Size 


Sampling  involves  a  compromise  between  precision  and  economy. 

A  parameter  can  only  be  defined  with  100  percent  accuracy  by  measur¬ 
ing  or  counting  the  whole  population.  The  smaller  the  sample,  the 
wider  the  margin  of  error  of  any  estimate  made  from  it.  A  homo¬ 
geneous  population  produces  samples  with  smaller  sampling  errors 
than  does  a  heterogeneous  population.  If  the  margin  of  error  that 
can  be  tolerated  for  a  particular  purpose  is  known,  it  is  possible 
to  calculate  the  minimum  sample  size  required  and  thereby  same  time 
in  the  collection  of  an  unnecessarily  large  sample.  Different  methods 
in  environmental  studies  can  be  used  depending  on  whether  the  samples 
are  to  be  measured  or  counted.  Both  involve  prior  collection  of  a 
pilot  sample  of  about  30  observations.  These  samples  must  be  collected 
by  the  same  methods  that  are  to  be  employed  in  the  full  sample;  and 
they  can  be  used  as  part  of  the  full  sample,  provided  they  do  not 
affect  the  selection  of  additional  items. 

The  formular  for  finding  the  required  sample  size  using  measured 
data  is:  n  =  (2i)2,  where  £  is  the  required  size  of  the  sample,  £ 
is  the  standard  deviation  of  the  pilot  sample,  d_  is  the  tolerable 
margin  of  error  at  a  specified  level  of  confidence,  and  z  is  a  con¬ 
stant  taken  from  the  z-table  (table  of  probabilities  associated  with 
values  of  z  in  a  normal  distribution)  corresponding  to  the  same  level 
of  confidence.  Therefore,  if  the  pilot  sample  has  a  standard  devia¬ 
tion  of  10  and  if  the  margin  of  error  that  can  be  tolerated  is  2.5  at 
the  95  percent  level,  then  the  required  sample  size  is  n  =  (2x10/2. 5)^ 
6k. 

The  formula  for  finding  the  sample  size  using  counted  data  is 
n  =  p%  x  q%  x  (z/d)2  percent,  where  £,  z_,  and  d_  have  the  same  mean¬ 
ings  as  in  paragraph  above  and  where  £  and  £  are  the  percentages  of 
the  pilot  sample  belonging  and  not  belonging  to  the  specified  category 
respect i vel y. 
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GEOGRAPHIC  COORDINATE  SYSTEM  USED  AT  RMA 


Several  different  coordinate  systems  have  been  used  in  the  past 
on  RMA.  Each  coordinate  system  has  its  own  merits.  However,  there 
exists  a  need  to  coordinate  all  of  the  data-gather ing  tasks  into  one 
standardized  reference  system.  The  identification  code  of  each  sample 
collected  w i  1.-1  be  located  with  reference  to  the  State  Plane  Coordinate 
System,  and  all  existing  data  will  be  converted  to  this  system. 

1  •  State  Plane  Coordinates  represent  a  system  set  up  for  each 
state  by  the  U.  S.  Department  of  Commerce,  Coast  and  Geodetic  Survey, 
to  be  used  by  the  private  land  surveyors  of  those  states.  The  system 
consists  of  a  rectangular  coordinate  grid  that  covers  the  entire  state, 
which  is  oriented  due  north  and  is  reduced  to  sea  level.  These  coord¬ 
inates  can  be  tied  into  a  more  universal  Geographic  Coordinate,  which 
utilizes  latitude  and  longitude.  Any  point  to  be  referenced  within 
the  state  can  be  identified  on  the  State  Plane  Coordinate  Grid.  This 
has  the  advantage  of  being  more  flexible  and  universal  when  tracing  data 
that  extends  beyond  the  Arsenal  boundaries. 

2.  RMA  Coordinate  System  -  This  system  was  established  in  1 942 
on  the  original  survey  of  the  Arsenal.  The  origin  of  the  system  was 
a  Geological  Survey  monument  named  Sand  1936  (now  destroyed).  The 
grid  is  oriented  to  magnetic  north.  Out  of  about  160  coordinate  points 
established  in  the  past,  only  eight  points  have  been  recorded.  This 
system  is  completely  independent  of  all  other  coordinate  systems  used 
in  the  area.  The  origin  of  this  grid  system  is  located  in  the  southwest 
corner,  such  that  all  X  and  Y  values  on  the  Arsenal  are  positive  numbers. 
Many  of  the  original  Arsenal  design  plans  use  the  RMA  Coordinate  System. 
This  system  is  inadequate  for  IR  needs  because  negative  coordinates 
would  occur  for  locations  beyond  the  Arsenal. 

3-  The  U.  S.  Rectangular  System  is  based  on  a  grid  consisting  of 
six  mile  wide  east-west  tiers  and  north-south  ranges.  Each  grid  cell 

is  a  township  comprised  of  36  one  square  mile  sections  which  are  num¬ 

erically  numbered  in  snake  fashion  beginning  in  the  northeast  corner. 

Each  one  mile  block  of  land  or  section  has  its  own  grid  system  and  is 
therefore  difficult  to  use. 

The  Section  Grid  is  a  coordinate  system  set  up  for  each  one 
mile  block  within  a  section  of  land.  The  southwest  corner  of  the 
section  is  arbitrarily  set  at  1,000  feet  north  and  1,00  feet  east. 

This  small  grid  is  used  with  very  little  control  surveying.  To  set 

up  this  grid,  only  the  boundaries  of  the  section  need  to  be  measured. 

A  control  survey  is  then  completed  to  locate  the  section  corners  on 
another,  more  universal,  system.  The  section  coordinates  can  then 
be  converted  to  the  control  survey  grid.  This  is  a  localized  grid 
system  which  is  easy  to  use  for  specific  tasks. 
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PRINCIPAL  COMPONENT  FACTOR  ANALYSIS 


Factor  analysis  is  a  statistical  technique  which  is  used  to 
reduce  a  square  symmetric  (mxm)  matrix  to  a  related  rectangular 
(mxp)  matrix  such  that  m>p.  The  initial  matrix  is  a  correlation 
coefficient  matrix  consisting  of  m  variables.  The  output  matrix 
presents  restructuring  the  pattern  of  correlations  into  p  clusters 
of  characteristics,  accounting  for  much  of  the  variance  of  the  m 
set  of  initial  variables. 

The  p  set  of  final  vectors  then  constitutes,  a  set  of  surrogate 
variables  (called  factors),  which  represent  the  m  set  of  variables. 

Within  IR,  we  will  deal  exclusively  with  principal  component 
factor  analysis  which  selects  an  orthogenal  rather  than  an  oblique 
(correlated  factor)  solution  and  which  assigns  all  of  the  variance 
for  interpretation  rather  than  only  part  of  it. 

Principal  components  method  is  preferred  because  of  ortho¬ 
gonality  and  the  uniqueness  of  the  solution.  There  is  only  one 
principal  components  solution  to  a  correlation  matrix,  thus  the 
factor  structure  is  not  likely  to  be  an  artifact  of  procedure.  The 
factors  in  the  solution  are  uncorrelated  and  are  statistically 
independent.  Therefore,  surrogate  variables  are  independent  by 
design.  If  there  are  N  observation  drill  sites  in  each  cluster 
and  each  site  provides  information  with  respect  to  m  variables, 
then  at  each  point  in  time  a  data  matrix  (N  x  m)  is  obtained. 
Simulation  of  time  trends  may  result  in  error  because  the  signifi¬ 
cant  correlations  among  the  variables  are  likely  to  be  interdependent. 
However,  the  correlation  structure  among  the  m  variables  can  be  redu¬ 
ced  into  p  orhoganal  components,  with  which  them  variables  are  associ¬ 
ated  by  means  of  linear  equations.  Time  trends  of  the  p  components 
may  now  be  analyzed  without  fear  of  independence  problems.  The' 
results  of  the  simulations  representing  p  components  can  be  recon¬ 
structed  back  to  the  original  variables.  For  example,  the  following 
flow  diagram  illustrates  the  steps  in  factorial  simulation: 
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Principal  component  analysis  can  be  used  in  problems  which 
involve  an  analysis  of  spatial  structure  of  contamination  and  spa¬ 
tial  taxonomy.  For  example,  an  area  of  interest  such  as  the  northern 
boundary  is  subdivided  into  N  subareas  or  clusters.  For  each  subarea, 
information  on_m  descriptive  variables  is  measured.  This  constitutes 
a  data  matrix  X  (N  x  m) .  The  questions  we  now  pose  are: 

1.  How  may!  the  internal  variations  of  the  m  variables  over  the 
entire  area  be  parsimoniously  described? 

2.  How  may  the  subareas  be  arranged  to  reveal  the  similarity 
and  difference  so  as  to  assign  a  meaningful  regional  structure? 

A  data  matrix  consisting  of  variables  measured  at  strategically  . 
located  sites  positioned  to  describe  the  environment  can  be  subjected 
to  correlation  analysis  to  reveal  how  the  descriptive  variables  are 
intercorrel ated  and  therefore  the  degree  to  which  they  report  redun¬ 
dant  information. 

Principal  components  analysis  can  be  used.  The  oringal  variables 
may  be  correlated  and  numerous.  The  patterns  of  variation  may  be  inde¬ 
pendent  and  confusing.  These  can  be  reduced  to  p  factors  which,  on  the 
other  hand,  are  few  and  independent  and  thus  describe  the  region  more 
parsimoniously.  Furthermore,  there  is  no  way  of  knowing  which  of  the 
variables  are  the  most  powerful  descriptive  tools;  whereas,  the  princi¬ 
pal  components  method  ranks  the  factors  in  descending  order  of  ability 
to  interpret  the  correlation  matrix  variance.  If  factor  -fj  interprets 
*42  percent  of  the  variance  and  factor  f£  3**  percent,  we  can  ascertain 
which  is  more  general  and  by  how  much. 

For  mapping  and  regionalizing,  it  is  desirable  to  be  able  to  esti¬ 
mate  the  strength  of  each  factor  within  each  tract.  For  this  purpose, 
the  matrix  A  (m  x  p)  which  gives  the  relationships  between  the  variables 
and  the  factors  is  combined  with  a  matrix  X  Z  (N  x  m)  of  standardized 
values  of  the  initial  variables.  The  result  is  a  matrix  of  dimension 
(N  x  p)  which  estimates  the  strength  of  each  factor  in  each  subarea  of 
observation.  This  is  called  the  factor  score  matrix  F.  The  sequence 
is  as  fol lows : 


Flow  Diagram  of  Principal  Components  Analysis 
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Each  row  of  the  factor  score  matrix  of  a  given  sampling  bore  hole 
(or  subarea)  across  the  independent  dimensions  represented  by  the  columns. 
It  is  therefore  actually  a  point  in  the  p-d imensional  space  of  the  fac¬ 
tor  structure.  One  may  therefore  evaluate  the  similarity  or  difference 
between  subarea  i  and  subarea  j  by  means  of  the  pythagorean  theorem  which 
evaluates  the  distance  between  any  two  points  in  a  p-dimension  space. 

(1)  d?j  =  :(ej]  -  e j ]  )2  +  (ej2  ~  ej2  +  •••  +  (eip  "  ejp  ^ 

Equation  (l)  may  be  evaluated  for  every  pair  of  subareas,  and  the 
set  of  {d?j}  will  then  form  the  matrix.  Analysis  of  this  matrix  -- 
especially  by  such  techniques  as  discriminant  analysis  —  will  cluster 
those  subareas  which  most  resemble  one  another  into  families  linked  into 
a  regional  hierarchy  of  subareas.  The  grouping  of  the  { d ? j }  by  discri¬ 
minant  analysis  will  be  presented  at  a  later  date. 

The  Procedure  of  Principal  Components  Analysis 

If  we  consider  the  matrix  Z  (N  z  m)  of  standardized  values  of  the 
N  subareas  on  the  m  variables,  we  may  regard  each  variable  as  a  vector 
in  a  N  dimensional  space  of  orthogonal  axes.  The  columns  of  the  matrix, 
in  this  case,  form  the  vectors  which  define  each  variable's  positions 
in  this  space.  Since  we  are  deal ing  with  a  standard i zed  value  matrix, 
all  vectors  are  of  equal  (and  unit)  length.  The  similarity  of  any 
variable  to  any  other  is  therefore  measured  by  the  angle  between  them. 

But,  the  cosine  of  the  angle  between  two  standardized  vectors  in  a 
space  of  N  dimensions  is  given  by  the  inner  vector  product  between 
them.  The  cosine  function  varies  between  -1  and  +1  and  is,  in  fact, 

(2)  cos  °ij  =  zj  .  zj 

the  correlation  coefficient  between  the  i£i2  and  jtll  variable.  Thus, 
we  find  R,  the  correlation  coefficient  matrix,  by  expression  (3). 

_  _/  _ 

(3)  R  =  Z  .  Z 

'V 

R  is  a  matrix  which  measures  the  angles  between  variable  vectors. 

The  principal  components  solution  starts  by  discovering  how  orthogonal 
axes  may  be  passed  through  this  pattern  of  intersecting  vectors  so  that 
each  axis  will,  as  far  as  it  is  possible,  coincide  with  the  center  of  a 
"cone"  of  vectors  (i.e.,  be  correlated  with  it).  This  is  analogous 
with  the  physicist's  problem  of  finding  the  resultant  of  a  set  of  force 
vectors.  In  the  F  matrix  case,  the  resultant  is  found  by  solving  the 
characteristic  equation  of  the  matrix. 
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In  ordinary  algebra,  an  expression  of  the  type 
1  C1*)  ay  =  gy 

implies  that  a  =  g.  In  matrix  algebra,  a  common  factor  on  each  side 
of  the  equation  does  not  imply  the  equality. 

'  (5)  AY  =*  Y 

In  expression  (5) ,  A  may  be  a  (m  x  m)  matrix,  and  Y  may  be  a 
(m  x  1 )  vector.  A  may  be  a  scalar,  and  it  will  then  be  true  that  the 
left  hand  term  and  right  hand  term  will  both  be  (m  x  1)  vectors;  yet 
A  ^  A..  Scalars  such  as  A  and  vectors  such  as  Y,  which  have  this 
property  with  respect  to  a  given  matrix,  A,  are  called  eigenvalues 
and  eigenvectors.  When  the  correlation  coefficient  matrix  R  is  fac¬ 
tored  into  eigenvalues  and  eigenvectors  (factor  ana  lysis),  the  eigen¬ 
values  are  orthogonal  and  have  the  properties  which  we  seek.  Further¬ 
more,  the  eigenvalues  prove  to  be  scaling  coefficients  such  that  they 
are  proportional  to  the  amount  of  variance  interpreted  by  each  factor. 
Any  factor  with  a  A>  1.0  interprets  more  variance  than  an  unmodified 
variable  alone,  while  any  factor  with  a  A< 1-0  is  of  relatively  weak 
explanatory  power.  Generally,  the  researcher  rejects  factors  with 
X<  1.0,  and  the  E A i  of  these  rejected  factors  is  proportional  to  the 
unexplained  variance  in  the  model. 

The  factor  loading  matrix  (A)  is  synthesized  from  the  eigenvalues 
and  eigenvectors  so  that 

(5)  R  A  =  A  A  . 

where  A  is  a  diagonal  matrix  whose  elements  on  the  principal  diagonal 
are  the  A;‘s. 

A  matrix  F  can  then  be  derived  whose  columns  are  vectors  fj,  f2, 
corresponding  to  the  scores  of  each  factor  in  a  certain  tract,  by  seek 
ing  some  relationship  based  upon  the  property  that 

(6)  z  =  A  f 

which  is  to  say  that  z  (a  column  vector  of  the  standardized  variables) 
for  the  tract  and  f  (a  column  vector  of  the  factors)  are  linearly  rela 
ted  through  A.  Then,  premultiplying  by  A  gives 

(7)  A'z  =  AAf 

and  finally 

(8)  7  =  (a'  A)'1  A'z 
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Expression  (8)  provides  the  means  for  estimating  the  values  of 
1  the  p  factors  for  all  the  tracts.  These  are  the  values  which  may  be 
1  mapped. 

Rotation  of  the  factor  structure  to  any  one  of  an  infinite 
'number  of  configurations  will  yield  different  resultants.  The  most 
useful  rotation  for  iR  is  that  which  uses  the  varimax  criterion. 

The  purpose  of  this  is  to  place  the  component  axes  in  the  unique 
position  where  as  much  of  each  factor  as  possible  is  interpreted  by 
■  strong  loadings  by  as  few  variables  as  possible.  This  rotation  sim¬ 
plifies  factor  structure  in  that  each  factor  axis  is  placed  as  para¬ 
llel  as  possible  to  a  limited  number  of  variable  vectors. 

Rotation  serves  a  very  useful  purpose  in  clarifying  t,he  inter¬ 
pretation  of  the  model  by  a  rule  of  parsimony.  Furthermore,  the  vari¬ 
max  rotation,  oke  the  principal  components  solution,  yields  a  unique 
solution  whose  factors  are  orthogonal  --  but  with  the  added  advantage 
of  simplicity  of  interpretation. 


FACTOR  ANALYSIS  OF  CONTAMINANTS  ALONG  THE  NORTHERN  BOUNDARY 

The  following  example  will  illustrate  the  use  of  principal  com¬ 
ponent  factor  analysis  for  (l)  data  reduction,  (2)  analysis  of  spa¬ 
tial  structure,  and  (3)  spatial  mapping. 

Twenty-four  bore  hole  samples  were  drilled  250  feet  apart  in  two 
parallel  east  to  west  lines  along  the  northern  boundary  at  RMA.  Chemi¬ 
cal  analysis  was  performed  on  the  groundwater  collected  from  each  sample 
and  tested  for  eight  compounds  (eight  variables).  The  raw  data  matrix 
[2k  x  8)  is  shown  on  Table  1,  in  which  the  rows  represent  the  bore  holes 
(samples,  N) ,  and  the  columns  represent  the  variables  (m) . 

Table  3  is  an  eight  by  eight  correlation  matrix  (M  x  M) ,  which 
represents  the  product  moment  correlation  between  each  of  the  eight 
compounds  analyzed.  When  the  number  of  variables  in  the  set  is  small, 
it  is  easy  to  evaluate  those  variables  which  correlate  highly  with  each 
other.  When  the  number  of  variables  increases,  it  becomes  difficult  to 
identify  those  variables  which  are  intercor related  and  group  together 
as  clusters.  These  groups  can  be  isolated  as  orthogonal  factors.  Table 
6  shows  three  factors  represented  as  columns.  The  variables  in  this 
(m  x  p)  matrix  are  shown  as  the  rows  where  the  factor  loadings  and  the 
factors  represent  their  relationships. 

The  cl  ion,  DIMP  and  the  Na  ion  load  up  on  factor  1;  therefore, 
they  must  also  be  highly  intercorrelated  and  occur  as  a  group  in  this 
data  set. 
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SAMPLES 


RAW  DATA  MATRIX 


VAR  I 


A- 1, 

A- 2 

A- 3 

A-ij 

A- 5 

A- 6 

A- 7 

A- 3 

D-l 

D-2 

D-3 

D-4 

D-5 

D-6 

D-7 

D-8 

UD- 

UD- 

UD- 

UD- 

UD- 

UD- 

UD- 

UD- 


1 

2 

4 

5 

6 

7 

8' 

3 


ABLES  /  /'-S'* 

PPM'  PPM  PPM  rP PM  "PPM, 
V2 


,'44541204550 
i'27630  001950 
.'33530002100 
07507500850 
1230&6304151 
122027306552 
150034307500 
;20037309552 
^28031302590 
!214;20801520 
09.507501300 
24526805503 
(52552505300 
36527002000 
24519501800 
22532508552 
1242214020  00 
124023403250 
44561501550 
32527003950 
S17034510  003 
128529002100 
j20  420  0  0  30  02 
|169j25306153 


.206407 
. G 432 07 
.043857 
.091207 
.404707 
.956507 
.127257 
.557057 
.043507 
.062757 
.061808 
.006407 
.186507 
.  04.4757 
.213407 
.008557 
.144207 
.045007 
.  05740;7 
.044357 
.158507 
.  062957 
.004658 


.6222000031 
.8609400290 
.4807301600 
.90.00920073 
.5904000000 
.80017  00  000 
.4015000000 
.9100820000 
.7712000290 
.9105300240 
. 0100660000 
.6702700000 
.5529000000 
.6507501200 
.7004800088 
.6701500010 
.6085700320 
.6584200010 
.49-19000000 
.7487000093 
.6280630000 
.9809600093 
.2603400012 


.306657.5202400000 
i  I  1  i 


TABLE  1 


VARIABLE 

MEAN 

ST. DEV. 

cl  “ 

1 

273.125011 

115.321591 

Na 

2 

289.125000 

117.819796 

SO  4= 

3 

407.333340 

273.620590 

NO2/NO3 

4 

.906667 

1.241141 

Hardness 

5 

506.249996 

203.855494  . 

pH 

6 

7.722917 

.200531  ‘ 

D1MP 

7  ' 

735.541687 

735.724869 

DCPD 

8 

181.250010 

394.055733 

TABLE  2 
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000005  ,947572 

000005  , 889357 


SAMPLES 


ROTATED  FACTOR  MATRIX 


VARIABLE 

1 

FACTOR 

!  -1 

3 

1  1  0*7  A 

cl'  -1 

Na  ^ 

S°4=  ? 

N02/N03  4 

Hardness  ^ 

pH  6 

DIMP  7 

DCPD  8 

.95432 
.  86626 
.07578 
-.36283  | 

.42726 
-.59616 
'.93632 
.02320 

. 05771 
-.32869 
-.93095 
-.86278 
-.87938 
.32518 
.23634 
.31009 

• liV/o 
-.04775 
-.09502 
-.10680 
-.00878 
-.52109 
-.12260 
.89032 

TABLE  6 


A-1 
A- 2 
A- 3 
A- 4 
A- 5 
A- 6 

A-7 
A- 8 
D-l 
D-2 
D-3 

D-4 

D-5 

d-6 

D-7 

D-8 

UD-1 

UD-2 

UD-4 

UD-5 

UD-6 

UD-7 

UD-8 

UD-3 


FACTOR  scores 


SE 

Factor  1 

1 

1 . 63071 

2 

.04633 

3 

. 13449 

4 

-1.46349 

5 

-.28705 

6 

-.63546 

7 

1.47940 

8 

-.50271 

9 

.27072 

10 

-.55679 

11 

— 1 . 44698 

12 

-.46282 

13 

2.46890 

14 

. 16364 

15 

-.35982 

16 

-.22932 

17 

-.21799 

18 

-.17469 

19 

2.12827 

20  - 

15308 

21 

-.48184 

22 

.00252 

23 

24 


-1.01234 

-.64515 


Var  1  •  cl 
2  =  Na 


|v 


7  =  D1MP 


Factor  2 
.  14408 
.92987 
.57079 
1.25532 
-.26986 
-1.16170 
-.57263 
-1.52047 
.78231 
1.02852 
1.19991 
-1.09225 
.09476 
.59346 
.69941 
-1.60303 
.51007 
.23713 
.10989 
.39788 
-2.18231 
1.04124 
.24606 
-1.44025 


ar  3  =  SO4 

4  =  NO2/NO2 

5  =  Hardness 


Factor  3 _ 

-.63347 
-.24377 
3.43871 
-.37236 
.07464 
-.26449 
.13603 
-.57213 
-.12060 
-.26528 
-.75781 
.01025 
-.85149 
2.25481 
-.03682 
. 04838 
.64720 
-.05952 
-.45420 
-.23061 
.18710 
-.93460 
-1.39069 
.39619 


Var  8  =  DCPD 


TABLE  7 


Factor  2  represents  the  S0]J  ion,  NC^/NO^ ,  and  hardness  variable 
group.  DCPD  behaves  independently  represented  as  Factor  3  in  this 
data  set.  It  is  possible  from  the  factor  loading  to  identify  which 
variable  (model  variable)  has  the  highest  correlation  with  the  respec 
tive  factor.  This  model  variable  can  then  represent  the  group  in 
future  tests  and  thus  permit  data  reduction  with  minimum  information 
loss. 

Table  4  gives  the  eigenvalues  which  can  be  interpreted  as  the 
.proportion  of  the  variance  in  the  data  matrix  which  is  explained  by 
each  of  the  factors.  By  weighting  the  proportion  each  eigenvalue  is 
of  the  total,  we  can  determine  the  relative  explanatory  power  of  each 
of  the  factors  in  terms  of  the  total  common  variance  in  the  data  set. 
Factor  1  accounts  for  42  percent  of  the  total  variance,  while  Factor 
2  explains  33  percent.  All  three  factors  together  explain  88  percent 
of  the  total  variance  in  this  example. 

Table  7  gives  the  factor  scores.  These  numbers  represent  the 
weight  of  each  of  the  factors  at  each  of  the  24  drill  sites.  By  map¬ 
ping  and  contouring  the  factor  scores  for  each  of  the  factors  separa¬ 
tely,  we  can  construct  maps  showing  multivariate  distribution  repre¬ 
senting  the  spatial  behavior  of  groups  of  contaminants  which  behave 
spatially  together.  See  map 

This  information  is  invaluable  in  (l)  tracing  each  of  the  con¬ 
taminant  migration  paths  independently,  (2)  determining  the  source 
of  contamination,  and  (3)  determining  the  position  and  direction  for 
further  sampl ing. 

SAMPLING  SITE  CHARACTERIZATION  &  DRILL  CORE  LOGGING 

OVERVIEW 

Soils  are  organized  natural  bodies  which  vary  markedly,  both 
horizontally  across  the  surface  of  land  and  vertically  as  seen  in 
cross  section.  The  vertical  section  provides  the  best  means  of 
examining  soil.  The  component  layers  that  form  the  basic  unit  for 
the  classification  of  soils,  in  turn,  provides  a  basic  framework 
for  the  analysis  and  mapping  of  the  behavior  of  the  pollutants  in 
the  ground. 

The  profile  form  which  expresses  the  overall  visual  impact  of 
the  physical  soil  properties  in  their  association  with  one  another 
serves  as  the  basis  for  the  following  key.  The  profile  is  regarded 
as  a  physical  system  and  those  physical  properties  capable  of 
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observation  and  record  which  may  incorporate  other  features  or  prop¬ 
erties,  either  physical,  chemical,  or  biological  are  used  as-  the 
critical  ones  to  distinguish  subtle  changes  in  the  environment  and 
contaminant  behavior. 

An  attempt  has  been  made  to  identify  and  code  most  of  the 
relevant  parameters  which  best  describe  the  physical  characteristics 
of  the  soil.  -  The  following  list  of  variables  are  to  be  measured  at 
each  sampling  site  in  the  comprehensive  survey  at  RMA.  This  list 
includes  a  coding  key  which  has  been  successfully  used  in  previous 
studies  and  should  be  used  to  enter  the  data  gathered  into  the  com 
puterized  data  bank.  The  scale  for  each  of  the  variables  describing 
topography,  site  drainage,  physical  soil  properties,  and  bedrock 
lithology  has  been  designed  so  that  the  scale  can  serve  as  input 
directly  into  numerical  analysis  and  computer  graphics  routines. 

The  key  is  flexible  in  that  if  necessary,  it  can  be  expanded  without 
detriment  to  those  portions  of  the  key  already  completed. 

PHYSICAL  VARIABLES 

Each  sampling  site  will  be  identified  by  a  ten  digit  index 
located  on  the  State  Plane  Coordinate  System.  The  first  five 
characters  of  the  index  represent  the  X  axis  of  the  east-west  State 
Plane  Coordinates  measured  to  the  nearest  foot,  and  the  last  five 
characters  represent  the  Y  axis. 

A.  Topographic  Variables 

1.  Elevation  is  the  height  of  a  given  point  above  mean 
sea  level.  Elevation  is  a  critical  variable  which  is  used  in  most 
calculations  to  provide  information  on  relative  position. 

Units  =  feet  (height  above  mean  sea  level) 

2.  Local  relief  is  the  difference  in  elevation  between  the 
highest  and  lowest  points.  It  measures  the  physical  size  of  the 
topographic  features. 

Units  =  feet  (difference  between  highest  and  lowest 

points) 

3.  Ruggedness  is  a  combination  of  gully  depth  measured 
in  feet  and  gul ly  count  per  100  square  yard  site.  Ruggedness  mea¬ 
sures  the  degree  of  direction  of  a  landform  and  provides  information 
about  the  degree  and  stage  of  landform  erosion.  It  also  indirectly 
measures  the  pattern  of  surface  water  flow. 
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(a)  Depth  in  feet 


>  1  foot  -  1 

i  1-6  feet  r -  2 

6-15  feet -  3 

i  ;  >15  feet  -  4 

'(b)!  Count  across  100  yards  square  site 

None -  0 

1-5 - 1 

5-10 - 2 

>  10 -  3 


4.  Slope  describes  the  angular  inclination  of  the  surface 
at  the  sampling  site.  It  is  measured  as  a  ratio  of  the  amount  of 
drop  (h)  in  100  feet  of  length  (C ) . 

Units  =  ratio  h  ;  (recorded  as  a  decimal  number) 

5.  Slope  geometry  represents  the  shape  of  the  sloping  sur¬ 
face  and  is  an  important  variable  in  characterizing  soil  erosion  or 
soil  accumulation  on  the  ground  surface. 

(a)  Shape 


Convex  -  1 

Concave  -  2 

Straight  -  3 

(b)  Relative  location 

Top -  1 

Middle  -  2 

Bottom  -  3 


6.  Aspect  describes  the  directional  orientation  of  a  topo¬ 
graphic  feature  with  reference  to  north.  It  is  most  convenient  to 
divide  the  compass  direction  into  16  divisions. 


N.N.E.  -  01 

N.E.  -  02 

E.N.E.  -  03 

E.  -  Qk 
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E.S.E. - 05 

S.E.  -  06 

5.5. E. . . -  07 

S.  - - —  08 

5.5. W.  — - 09 

S.W.  . 10 

w.s.w. - - -  n 

W.  - 12 

W.N.W.  - 13 

N.W. . . — .  1^ 

N.N.W.  -  15 

N.  - 15 


7 .  Drainage  is  characterized  by  two  different  components. 
Site  drainage  describes  the  movement  or  potential  for  movement  of 
water  and  moisture  —  the  amount  of  water  contained  in  the  soil. 

(a)  Site  drainage 


Shedding  -  01 

Normal  -  02 

Receiving -  03 

Regular  flooding,  Short 

Term -  04 

Regular  flooding.  Long 

Term - 05 

Occasional  flooding,  Short 

Term -  06 

Occasional  flooding,  Long 
Term -  07 


(b) 


Mo isture 

Dry  (becomes  darker  when  wet) 
D,  (air  dry)  - 

o2 - 

d3 - - - 

Moderately  moist 

MM,  - 

mm2 - 

MM3  (wilting  point)  - 

Moi  st 

M1  - 

m2 - 

M3  (field  capacity)  - 


01 

02 

03 


04 

05 

06 


07 

03 

09 
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Wet 


(material  wi  1 1  wet 
and  stick  to  fingers 
when  molded) 


1.  Depth  to  sample  in  the  vertical  distance  measured  from 
the  ground  surface  to  the  top  and  the  bottom  of  the  sample. 

units  =  feet  (measured  to  the  nearest  hundreth  of  foot) 

2.  Horizon  is  a  layer  within  the  soil  profile  haying  mor¬ 
phological  characteristics  and  properties  different  from  those  layers 
which  occur  below  and/or  above  it.  Horizons  are  usually  parallel  to 
the  land  surface,  except  where  tongues  of  material  from  one  horizon 
may  penetrate  vertically  into  the  neighboring  horizon.  The  horizons 
are  recognized  by  the  nature  of  their  organization.  Aq  results  largely 
from  biologic  processes.  The  A  horizon  results  from  biologic  processes 
and  those  physico-chemical  processes  that  occur  at  surfaces.  The 
horizon  results  from  physio-chemical  processes  which  occur  below  sur 
faces  that  are  relatively  deep  seated.  C  horizons  are  recogn i zed  by^ 
their  lack  of  pedologic  development  and  remains  of  geologic  organiza¬ 
tion.  The  D  horizons  result  from  some  earlier  cycle  of  soil-forming 
process;  and  thus,  their  pedologic  organization  contrasts  with  that 
found  in  the  solum  of  the  soil  with  which  they  are  associated. 

Weathered  bedrock  is  an  example. 

Eight  categories  will  be  recognized  in  this  study.  They 
are  given  in  descending  order:  the  AQ  horizon  is  the  surface  soil 
layer  and  is  composed  of  organic  material;  the  A]  horizon  is  the  top¬ 
most  mineral  layer  with  a  high  content  of  humified  organic  matter, 
relatively  dark  in  color  and  has  a  maximum  biological  activity  tor 
any  given  soil  profile;  the  A2  horizon  is  the  zone  of  maximum  leach¬ 
ing  (eluviation) ,  it  may  be  defined  by  color  difference.  It  may  be 
bleached  to  a  white  color  or  sporadic  by  bleaches  which  appear 
irregular  or  as  blocks  less  than  one-fourth  inch  in  thickness  at 
the  interface  with  the  B  layer.  The  B  horizon  is  a  transition  zone 
between  the  A  and  B  horizons  though  accumulation  of  material  (jllu- 
viation)  starts  to  take  place.  The  B  horizon  is  characterized  by  a 
concentration  of  clay  and/or  iron  and/or  aluminum  and/or  translocated 
organic  material  and/or  having  a  structure  or  consistency  different 
from  A  or  C.  B  also  has  stronger  colors.  The  B2  horizon  is  the  zone 


B.  Soil 


W1  — . 

w2  - 

W3  - 

Physical  Properties 


10 

11 

12 
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of  maximum  illuviation  of  mineral  and  organic  matter;  the  B-C  horizon 
lis  a  transition  zone  between  the  illuvial  B  horizons  and  the  uncon 
Undated  material  of  the  C  horizon.  The  C  horizon  ,s  unconsolidated 

material  which  may  or  may  not  be  of  the  same  parent  ma  ®r'a  r  V  , 

the  upper  horizons  are  formed.  It  is  not  in  the  zone  of  major  biological 

activities;  the  D  horizon  is  weathered  bedrock. 

Horizon 


3.  Horizon  thickness  —  represents  the  size  of  the  horizon, 
units  =  feet  measured  to  nearest  .01  of  a  foot. 

4.  Horizon  boundary  is  a  zone  of  change  occurring  between 
soil  horizons  ranging  from  the  sharp  and  clearly  defined  to  diffuse 
The  boundary  defines  the  nature  of  the  change  and  is  specified  by  two 
terms  —  one  a  measure  of  the  thickness  or  width  of  the  transition 
zone  between  horizons,  the  other  a  measure  of  its  shape. 

(a)  Shape 

1  -  Even  =  boundary  is  almost  a  plane  surface  _ 

2  -  Wavy  =  boundary  waves,  depressions  formed  are 

wider  than  their  depths  - - -  .  , 

3  -  Irregular  =  boundary  waves,  depressions  formed 

are  deeper  than  their  widths 

4  ~  Broken  =  boundary  is  discontinuous 

(b)  Sharpness 

1  -  Sharp  =  <  2  cm 

2  -  Clear  =  2  -  5  cm 

3  -  Gradual  =  5  “  10  cm 

4  -  Diffuse  =  >  1 0  cm 

5.  Color  is  defined  in  terms  of  hue,  value,  and  chroma  using 
the  Munsell  soTTTolor  charts.  Hue  is  a  function  of  the  wave  length  o 
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the  light.  Value  is  the  intensity  of  color  and  chroma  the  relative 
purity  or  saturation  of  a  color.  Color  should  be  recorded  for  moist 

soil  only. 

Color  (moist) 

1  Hue 


2,  5  YR -  32 

5  YR - 

7,  5  YR -  8 

10  YR -  4 

2,  5  Y -  2 


2  Value 

1-2-3-4-5-6-7-8 


3  Chroma 

1-2-3-4-5-6-7-8 

6.  pH  measures  the  hydrogen  ion  concentration  of  soil  .  \ 

is  a  fundamental"!  property  of  both  physical  and  chemical  charactei  is 
tics  of  both  the  organic  and  inorganic  constituents  of  the  soil. 

units  =  dimensionless  number  of  nearest  .1 


7.  Soil  texture  is  a  measure  of  the  size  and  proportion  of 
individual  gain  present  in  the  soil.  Texture  can  be  measured  by  the 
behavior  of  a  small  handful  of  soil  when  moistened  and  kneaded  intoa 
ball  and  then  pressed  out  between  thumb  and  forefinger.  The  approxi¬ 
mate  percentage  content  of  clay  (particles  less  than  0.002  mm  in  dia- 
meter)  and  silt  (particles  between  0.02  and  0.002  mm  in  diameter)  and 
sand  (range  from  2  mm  to  0.02  mm  diameter  with  an  arbitrary  separation 
of  coarse  and  fine  sand  at  0.2  mm  particle  diameter)  provide  a  guide 
for  classification. 


a.  The  field  method  of  determining  texture  soil  is  to 
take  a  sample  of  soil  sufficient  to  comfortably  fit  into  thepalm  of 
the  hand.  The  soil  is  moistened  with  water,  a  little  at  a  time,  and 
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kneaded  until  the  ball  of  soil,  so  formed,  just  fails  to  stock  to  the 
fingers.  More  soil  or  water  may  be  added  to  attain  this  condition, 
which  is  known  as  the  sticky  point  and  approximates  field  moisture, 
capacity  for  that  soil.  Kneading  and  moistening  are  continued  until 
there  is  no  apparent  change  in  the  soil  ball,  usually  a  working  time 
of  one  to  tw6  minutes.  The  soil  ball,  or  bolus,  is  now. ready  for 
shearing  manipulation;  bu  t  the  behavior  of  the  soil  during  bolus 
formation  is  also  indicative  of  its  texture.  The  behavior  of  the 
bolus  and  of  the  ribbon  produced  by  shearing  (pressing  out)  between 
thumb  and  forefinger  characterizes  the  texture. 

b.  A  much  more  reliable  method  of  determining  soil 
texture  is  by  sieving  the  soil  samples. 

Textural  classes 


S  Sand  -  01 

LS  Loamy  sand -  02 

CS  Clayey  sand -  03 

SL  Sandy  loam -  04 

FSL  Fine  sandy  loam -  05 

GL  Gravelly  loam -  06 

L  Loam -  07 

SCL  Sandy  clay  loam -  08 

GCL  Gravelly  clay  loam  --  09 

CL  Clay  loam -  10 

SI  CL  Si  lty  clay  loam -  1 1 

CG  Gravel  ly  clay -  12 

FSCL  Fine  sandy  clay  loam  -  13 

SC  Sandy  clay -  14 

SIC  Silty  clay -  1 5 

LC  Light  clay -  16 

MC  Med  i urn  clay -  17 

HC  Heavy  clay - ; -  18 


8.  Stones  are  measurable  surface  features.  The  abundance 
is  measured  as  a  percentage  of  cover;  size  consists  of  eight  categories 
ranging  from  gravel  (0.2  -  1  cm)  to  large  boulders  (<  600  cm);  shape 
can  be  angular,  subangular,  rounded,  shaley,  or  tabular;  and  the 
lithology  of  stones  is  given  in  12  petrological  types. 
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(a)  Abundance 

Stoneless 
SI ightly  stony 
Stony 

Very  stony 
Extremely  stony 
|  Stone  dominant 

(b)  Size 

Gravel 
Sma  1 1 
Hed i urn 
Large 

Very  large 
;  Small  boulders 
Medium  boulders 
Large  boulders 

(c)  Shape 


{%  cover) 


<  1  - .  0 

1  "  5  - .  1 

5-20 -  2 

20  -  50  .  3 

50  -  75  -  * 

>  75  — . -  5 

(cm) 

0.2-1  . .  1 

1  -  5 -  2 

5  .  \  -  3 

10-20 .  4 

20  -  60  -  5 

60  -  200  -  6 

200  -  600  -  7 

>  600  -  8 


Angular  -- 
Subangul ar 
Rounded  -- 

Shaley  - 

Tabular  -- 


(d)  Lithology 

FI  int - 

Soapstone  - 

Limestone  - 

Mudstone  - 

Siltstone  and  shale 

Sandstone  - 

Ironstone  - 

Quartzite  - 

Granite  - 

Gneiss  or  schist  — 

Ryhol ite  - 

Andesite  - 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 
1 1 
12 
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9.  Roots  describe  the  presence  of  vegetative  matter  in  the 
soil  measured  in  terms  of  quantity  and  size.  Decay  of  roots  provides 
channels  of  water  flow  within  the  soil. 

(a)  Quantity 


Absent  -  0 

Present  (trace)  -  1 

Very  low -  2 

Low -  3 

Moderate  -  4 

High -  5 

(b)  Size 

Sma  11  -  1 

Sma  1 1  -  Med  ium -  2 

Smal  1  -  Large -  3 

Med i urn  -  Sma 1 1  -  4 

Med ium  -  5 

Medium  -  Large  -  6 

Large  -  Smal  1  -  7 

Large  -  Medium -  8 

Large -  9 


10.  Structure  is  concerned  with  the  arrangement  of  all  soil 
particles.  Particles  refer  to  all  peds  and  nonpeds  (see  both  Apedal 
and  Peds).  Structure  may  be  described  in  terms  of  the  three  charac¬ 
teristics  of  the  intensity,  geometry,  and  size  of  the  soil  aggregates 
as  foil ows : 

(a)  Intensity  expresses  the  degree  and  strength  of  soil 
aggregation  and  may  vary  with  the  moisture  status  of  the  soil.  For 
purposes  of  standardization,  intensity  should  be  determined  when  the 
soil  is  inthe  moderately  moist  to  moist,  moisture  stages.  Four  inten¬ 
sities  of  structure  are  usually  defined  as  follows: 

(1)  Structureless  in  which  there  is  no  observable 

aggregation  and  the  soil  is  massive,  if  coherent,  and  single  grain,  if 
not  coherent.  (Structureless  is  the  same  as  Apedal.)  -  0 

(2)  Weak  in  which  some  peds  are  discernible  and 

when  disturbed,  less  than  one-third  of  the  soil  material  is  found  to 
consist  of  peds.  -  1 


33 


.  (3)  Moderate  in  which  peds  are  clearly  seen  and 

| when  disturbed,  one-third  to  two-thirds  of  the  soil  material  is  found 
'to  consist  of  peds.  - 2 

(4)  Strong  in  which  peds  are  clearly  seen  and 
when  disturbed,  more  than  two-thirds  of  the  soil  material  is  found 
to  consist  of.  peds.  - 3 

NOTE  that  both  moderate  and  strong  grades  of 
structure  qualify  a  soil  for  Peds  evident,  whereas  the  structure¬ 
less  state  and  weak  grade  of  structure  qualify  a  soil  for  Peds,  few, 
if  any. 

(b)  Geometry 


Angular  blocky  or  cubic  (ABL)  -  1 

Subangular  blocky  or  cubic  (SABL)  -  2 

Angular  prismatic  (APR)  -  3 

Subangular  prismatic  (SAPR)  -  4 

Polyhedral  >  six  sides  (PH)  -  5 

Paral lelepedal  (PLL)  -  6 

Platy  (PL)  - - —  7 

Spheroidal-solid  (SSPH)  - 8 

Spheroidal -vi s ible  voids  (VSPH)  -  9 


(c)  Size 

Given  in  inches  or  fraction 

11.  Plasticity  is  a  characteristic  of  wet  soils.  It  is  the 
ability  of  soil  to  change  its  shape  as  a  result  of  stress  and  to  keep 
that  shape  after  the  removal  of  the  stress. 


Nonplastic  - 0 

SI ightly  plast ic  -  1 

Plastic  -  2 

Very  plastic -  3 


12.  Fr iab? 1 ity  is  a  characteristic  of  moist  soils.  It  is  a 
measure  of  resistance  to  crushing. 


Loose -  0 

Very  friable -  1 

Friable  -  2 

Firm -  3 

Very  firm -  4 

Extremely  firm -  5 


13.  Porosity  -  a  porous  soil  material  is  one  characterized 
by  the  presence  of  voids  in  the  soil  moss.  It  is  the  volume  of  the 
total  soil  bulk  not  occupied  by  solid  particles.  The  first  component 
of  porosity  is  void  quantity,  which  is  the  percentage  of  nonsolid  space 
ranging  from  nonporous  (<  1%)  to  extremely  porous  (>  50%).  The  second 
component  is. pore  size,  measured  in  millimeters.  Porosity  is  largely 
a  function  of  the  arrangement  of  soil  particles  into  aggregates,  the 
type  of  soil  particles,  amount  of  organic  material,  and  degree  of 
compact  ion. 

(a)  Void  quantity 


Nonporous  <  -  1  0 

Slightly  porous  1  ~  5 -  1 

Porous  5  "  10 -  2 

Porous  10  -  20 -  3 

Very  porous  20  -  50  -  4 

Extremely  porous  >  50 -  5 

(b)  Pore  size 

Very  fine  <  -  0.5 -  1 

Fine  0. 5  "  1  -  2 

Med  ium  1  ~  3 -  3 

Large  3  -  5 .  ** 

Very  large  >  5 -  5 


14.  Mott  1 ing  is  the  presence  of  spots  or  alternate  streaks 
of  color  that  differ  from  the  overall  color  of  the  soil.  Mottling 
occurs  in  zones  of  alternating  periods  of  good  and  poor  aeration. 

The  differences  in  oxidation  states  of  iron  and  magnesium  result  in 
mottled  appearances.  Mottles  are  considered  to  be  masses,  blobs,  or 
blotches  of  subdominant  colors,  such  that  color  places  them  in  a 
different  value/chroma  rating  to  the  matrix  color.  The  glaze  or 
color  skin  which  may  occur  on  the  outside  of  some  peds  is  not  con¬ 
sidered  as  mottling.  A  soil  is  only  considered  to  be  mottled  when 
10%  or  more  of  mottles  occur  in  the  soil  mass.  For  mottled  soils, 
the  value/chroma  rating  of  the  dominant  matrix  color  is  taken.  Mot¬ 
tling  has  three  components:  abundance,  measured  as  a  percentage; 
size  in  mm;  and  sharpness  of  the  boundary  between  mottled  and  non- 
mottled  soil,  measured  in  mm. 

(a)  Abundance  {%  cover) 


None  <  2 -  0 

Few  2-5 -  1 
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- 

1 

Common 

5-20 - 

3 

I 

Definite  matrix  color 

20  -  50  - 

4 

No  matrix  color 

>  20  — . 

5 

Definite  matrix  color 

>  50 - 

6 

...  -■ . 

-(b) 

Size 

(mm) 

Extremely  fine 

<  1  - 

-  1 

w 

Very  fine 

1  -  2 - 

2 

Fine 

2  _  5 - — 

-  3 

Med i urn 

5  -  15 - 

-  4 

Large 

>  15 - 

-  5 

(c) 

Sharpness 

(mm) 

Diffuse 

>  2  — . — 

-  1 

Cl  ear 

<  2 . . 

-  2 

Sharp 

(knifedged)  — 

-  3 

15.  Secondary  precipitates  result 

from  the  upward 

capi 1 lary 

movement  and  removal  of  soil  water  through 

evaporation  at 

the  ground 

surface.  This  process  is  important  in  contaminant  movement  and  con 

- 

centration  in  the 

soil. 

i  .  _■ 

(a) 

Quant i ty 

‘ 

None 

-  0 

- 

Few 

<  2% . 

-  1 

Common 

2-5% - 

-  2 

Common 

5  -  20%  - 

-  3 

Abundant 

20  -  50%  - 

-  4 

Extra  abundant 

>  50% . - 

-  5 

(b) 

Compos i t ion 

Lime - 

-  1 

Magneste  - 

--  2 

Gypsum  -  - 

—  3 

/, 

Silica  - 

H 

Fe  oxides  - 

—  5 

Fe  oxides  +  humus  - 

6 

i 

Mn  oxides  - 

--  7 

•  i 
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(c)  Type 


Irregular  concretions - j 

Regular  concretions  -  2 

Layers  -  3 

1 6 .  ' Bedrock 

(a)  Rock  structure 

(1)  Intensity 

W.  Weak -  1 

M.  Moderate - -  2 

S.  Strong -  3 

(2)  Geometry 

B.  Banded -  01 

S.  Scoriaceous -  02 

A.  Amygdaloidal  -  03 

P.  Porphyritic -  04 

L.  Laminated -  05 

M.  Massive -  08 

C.  Cleared -  07 

S.  Sheared - ; -  08 

F.  Contorted -  °S 

X.  Cross-bedded  -  10 

(3)  Size 

V.  Very  coarse -  1 

C.  Coarse -  2 

F.  Fine -  3 

VF  Very  fine -  ^ 

(b)  Rock  type 

Consolidated  sediments  (0  -  19) 

CG  Conglomerate -  01 

QT  "Quartzite -  02 

SS  Sandstone -  03 

ST  Siltstone -  04 

MS  Mudstone -  05 

CS  Claystone -  06 

SH  Shale -  07 

CH  Chert - -  08 
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GW  Greywacke -  09 

AK  Arkose -  1° 

LS  Limestone - 

IS  Ironstone - ■ -  12 

Tl  Tillite  — - 13 

Igneous  rocks  (20  -  33) 

GR  Granite -  20 

GD  Granod  ior  i  te -  21 

D I  Diorite -  22 

RH  Rhyol  i  te -  23 

BA  Basalt  — -  24 

AG  Agglomerate - -  25 

TU  Tuff . 26 

SC  Scorria -  27 

DO  Dolerite -  28 

DB  Disabase -  29 

PG  Pegmatite -  30 

AP  Apl  i  te -  31 

RQ  Reef  quartz -  32 

PY  Porphyry -  33 

PT  Porphyrite -  34 

Metamorphic  rocks  (40  -  59) 

SL  Slate - - - - . 

PH  Phyl  1  ite -  ^1 

HF  Hornfels - ^2 

SH  Schist  -  ^3 

GN  Gneiss -  ^ 


INSTALLATION  RESTORATION 


INSTALLATION  RESTORATION 


INSTALLATION  RESTORATION 


SAMPLING  PLAN  FOR  RMA 


Dynamic  Cluster  Sampling 

An  attractive  attribute  of  dynamic  sampling  is  that  as  informa¬ 
tion  is  gathered,  it  can  be  used  to  direct  further  sampling  without 
biasing  the  probability  properties  of  the  sample. 

Data  on  soils  and  groundwater  will  be  gathered  from  drill  sites 
patterned  in  clusters  of  16  bore  holes  per  cluster,  arranged  in  a 
i,  x  k  square  grid.  Each  sample  will  be  surveyed , .  us  i  ng  the  x  and  y 
coordinates  of  the  State  Plane  Grid  System  which  in  addition,  will 
serve  as  the  identification  code  for  each  of  the  samples.  The  ele¬ 
vation  will  also  be  recorded.  Field  data  on  land  surface  topography, 
physical  properties  of  soils,  groundwater  levels,  and  bedrock  will  be 
logged  directly  onto  the  coding  forms.  Soil  samples  will  be  collected 
at  the  appropriate  depth;  sent  to  the  soils  laboratory  where  it  will 
be  divided  into  two  parts.  One  part  of  the  sample  will  go  to  the  MALD 
(Material  Analysis  Lab  Division)  for  chemical  analysis  of  various  con¬ 
taminants.  The  other  portion  will  be  analyzed  in  the  soi 1  laboratory 
for  soil  moisture,  mechanical,  and  mineral  properties..  While  the 
analyses  are  being  performed,  a  second  cluster  of  16  sites  will  be 
drilled  and  sampled.  The  results  from  the  first  cluster  of  samples 
will  be  mapped  to  determine  the  positioning  of  the  third  cluster  of 
16  sampling  sites.  The  cluster  will  be  located  in  the  direction  of 
the  maximum  concentration  of  the  contaminant  being  mapped.  The  analy¬ 
zed  results  from  the  second  cluster  will  likewise  be  used  to  determine 
the  positioning  of  the  fourth  cluster,  etc.  This  technique  of  mapping 
each  contaminant  and  positioning  the  cluster  of. samples  in  the  direc¬ 
tion  of  maximum  concentration  will  permit  tracking  the  pollution  plume 
to  its  origin  with  minimum  exploratory  drilling.  Periodically,  the 
cluster  will  be  combined,  and  analysis  will  be  performed  on  the  total 
combined  samples.  This  will  yield  an  increase  in  the  sample  size  and 
an  increased  relative  homogeneity  of  the  elements  being  sampled.  By. 
combining  the  cluster  into  a  total  sample  one,  in  effect,  is  evaluating 
the  regional  pattern  of  the  pollutants  as  compared  to  more  local,  larger 
scale  behavior  within  each  cluster.  However,  sampling  theory  requires 
that  all  the  elements  of  the  cluster  be  sampled  for  the  results  to  be 
stat i st ical  1  y  valid. 

It  is  possible  that  several  plumes  will  be  encountered, 
originating  from  a  different  source  and  containing  their  own 
grouping  of  compounds.  When  this  occurs,  each  plume  will  be 
to  locate  all  sources. 


each 
un i que 
traced 


b2 


On-going  exploratory  drilling  along  the  northern  boundary  of 
RMA  has  revealed  the  existence  of  at  least  two  contaminated  ground 
water  plumes.  One  plume  contains  DIMP,  Cl",  and  Na  compounds.  The 
other  is  a  DCPD  and  pesticide  plume.  Because  of  the  high  priority 
placed  upon  containment  of  contaminant  migration,  this  part  of  the 
Arsenal  would  be  a  logical  place  to  start  the  sampling  program.  An 
alternate  site  would  be  at  the  known  contaminant  hot  spots. 

.i 

1 

Sampl ing  Grid  Spacing 

A  test  sampling  grid,  spaced  250  feet  apart,  will  be  constructed 
initially  to  pretest  the  sampling  method  and  data-handl ing  procedure. 
The  distance  between  the  sampling  bore  holes  will  be  set  as  more  data 
becomes  available  on  the  diffusion  rates  of  contaminant  concentration 
with  distance.  The  mathematical  relationship  between  change  in  con¬ 
centration  level  per  unit  change  in  distance  will  be  calculated  for 
each  contaminant,  and  the  grid  spacing  will  be  adjusted  to  the  opti¬ 
mum  distance.  Preliminary  tests  show  that  the  best  fit  regression 
equation,  equating  the  change  in  DIMP  concentration  with  distance 
measured  from  the  center  of  the  DIMP  plume  is  Y  =  A  +  B  log  X.  The 
DIMP  concentration  values  in  the  inclosed  sample  contained  DIMP 
added  from  secondary  sources  carried  by  groundwater.  This  caused 
additive  harmonics,  which  should  be  filtered  out  before  a  usable 
equation  can  be  obtained. 

Sampl i no  Hor i zons 

One-foot  sections  of  the  sampling  core  will  be  collected  to 
serve  as  the  samples  for  both  mechanical  and  chemical  analysis. 

Soil  samples  will  be  taken  at  six  depths,  so  as  to  increase  the 
likelihood  of  intercepting  the  zones  of  contaminant  concentration. 
These  zones  are:  (1)  The  A  soil  horizon,  which  is  the  zone  of  con¬ 
taminant  built  through  (a)  bioaccumulation  in  the  soil  organic 
matter;  (b)  upward  migration  of  contaminants  by  capillary  action 
and  surface  evaporation  of  water;  and  (c)  accumulation  of  pollutants 
on  the  ground  surface  by  both  wind  or  surface  water  transport,  (  ) 

The  B  soil  horizon  which  is  the  zone  of  illuviation  or  accumulation, 
of  compounds  leached  from  the  upper  horizons;  (3)  A  point  mid-way  in 
the  zone  of  aeration  representing  the  transition  zone;  (4)  The.capil 
lary  fringe  immediately  above  the  ground  water  table;  (5)  A  point 
below  the  top  of  the  water  table;  and  (6)  At  the  bedrock  surface.  . 
'The  soil  samples  taken  in  the  zone  of  saturation  will  provide  a  ratio 
of  the  contamination  in  the  soil  compared  with  that  moving  in  the 


i  groundwater.  Additional  samples  will  be  taken  at  the  top  of  clay 
payers  or  lenses  whenever  they  are  encountered,  to  determine  the 
effect  of  clay  as  a  point  of  contaminant  concentration.  Clay  layers 
in  the  regolith  may  behave  as  selective,  physical,  chemical,  or  osmo 
-tic  filter  aided  by  the  greatly  reduced  rate  of  downward  percolation 
of  groundwater  and  the  possible  presence  of  clays  having  a  high  ionic 
exchange  capacity.  Groundwater  will  also  be  sampled  at  similar  depths. 

The  typical  profile  shown  in  the  diagram  was  constructed  by 
'  averaging  the  horizon  of  each  of  the  soil  series  occurring  at  RMA. 

The  A  horizon  ranges  from  0-9  inches  in  thickness,  while  the  B  hori¬ 
zon  ranges  from  S  -  26  inches. 

In  several  places  on  the  Arsenal,  the  B  horizon  does  hot  occur; 
and  the  soils  grade  directly  from  the  A  to  the  C  horizons.  This 
usually  occurs  in  younger,  poorly  developed  azonal  soils.  In  order 
to  avoid  missing  data  for  sites  where  horizons  are  absent,  all  sampl¬ 
ing  will  consistently  be  taken  at  uniform  depths  predetermined  to. 
coincide  with  the  soil  profile  horizons.  These  depths  are  given  in 
the  schematic  illustrating  the  sampling  scheme  for  a  typical  boring. 

Two  water  samples  will  be  collected  at  each  drill  one,  one  approx i 
mately  three  feet  below  the  water  table  leve,  the  other  at  bedrock 
surface. 

Numerical  Analysis  and  Happing 

Each  of  the  contaminants  will  be  mapped  individuall  for  each  of 
the  six  horizons  sampled,  thus  providing  a  three-dimensional  view  of 
the  contaminant  distribution.  The  maximum  concentration  encountered 
at  each  site  will  also  be  mapped  to  show  a  plain  view  of  the  highest 
levels  occurring  on  the  Arsenal. 

Principal  component  factor  analysis  will  be  peformed  to  determine 
the  correlations  and  behavior  of  the  contaminants  as  groups  or  associa¬ 
tions.  The  information  will  be  obtained  from  the  factor  loadings,  which 
will  identify  new  groupings  indicating  that  new  pollution  plumes  were, 
encountered.  The  maps  of  the  factor  scores  for  each  of  the  factors  will 
provide  ins i ght  on  the  spatial  distribution  of  each  of  the  contaminant 
groups  individually  and  will  serve  to  identify  the  location  and  direc¬ 
tion  of  the  plume,  thus  indicating  where  the  next  sampling  cluster. is 
to  be  positioned.  Factor  analysis  of  the  total  set  of  variables  will, 
through  the  correlation  coefficient,  provide  the  linkage  web  between 
the  contaminants  and  the  physical  ecosystem. 


The  relationships  and  the  linkages  will  be  mappable  through  the 
factor  scores.  Multivariate  discriminate  analysis,  using  the  factor 
scores  as  surrogate  weights,  will  serve  as  a  comparison  of  similarity 
and  differences  between  different  regions  representing  either,  dif¬ 
ferent  contaminated  areas  of  different  source  areas.  The  multivariate 
disciminant  function  can  also  serve  as  a  model  for  simulation. 
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Figure  10 


